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ABSTRACT 


OLD  BARHYARD:  A  CROSS  SECTION  CODE  FOR  SCHOOL  USE 


Charles  <i *  Bridgman 
Daniel  P.  Cannon 
Ernest  Park  Sims 
Robert  H.  Hansen 


v  OLD  BARNYARD  is  a  digital  computer  code  which  calculates  few  group 
fast  and  thermal  neutron  cross  sections  and  constants.  The  first  of 
two  sequential  chains  calculates  the  neutron  flux  needed  to  determine 
the  group  cross  sections.  The  code  uses  the  moments  method  to  find 
the  fast  flux  and  solves  the  Wilkins  equation  numerically  for  the 
thermal  flux.  The  code  is  written  in  Fortran  for  a  small  computer, 
the  IBM  1620,  and  is  intended  for  Use  by  students  and  professors  in 
support  of  classroom  assignments. 
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problems  generally  require,  as  input,  group  cross  sections  and  cross 
section  related  constants  such  as  the  diffusion  coefficient  or  neutron 
age.  However,  these  problems,  unlike  their  counterparts  in  the  real 
world,  usually  involve  only  one  or  a  few  energy  groups.  This  few 
group  distinction  has  always  made  the  acquisition  of  group  cross  sec¬ 
tions  a  particularly  frustrating  experience.  One  borrows  on  experi¬ 
mental  value  from  here,  a  calculated  value  from  there,  and  combines 
them  with  some  fast  cross  sections  collapsed  by  hand  from  a  many  group 
set  (which  inevitably  doesn't  contain  all  the  elements  of  interest.) 
Although  sophisticated  methods  are  always  available,  data  for  one 
problem  of  a  10  problem  assignment  rarely  warrants  the  effort  to 
obtain  them;  nor  does  the  accuracy  required  in  a  school  problem  war¬ 
rant  their  sophistication. 

The  code  described  in  this  report  is  an  attempt  to  provide 
reasonably  good  cross  sections  from  a  single,  easily  used,  and  un¬ 
sophisticated  (therefore  inexpensive)  source.  The  code  is  written  in 
a  version  of  the  Fortran  language,  the  Kingston  Version  of  Fortran  II 
(Kindstran.)  Kingstran  has  not  received  wide  publicity  and,  I  find, 
is  in  limited  use  at  other  schools.  However  the  language  :.'.ay  be  trans¬ 
lated  to  Fortran  IV  with  only  very  minor  modifications.  The  £FIT  ver¬ 
sion  has  been  compiled  for  the  IBM  1620,  a  digital  computer,  installed 
ub  146  (55%)  of  the  268  schools  listed  in  the  1964  Oak  Ridge  Institute 
of  Nuclear  Studies  report,  EDUCATIONAL  PROGRAMS  AND  FACILITIES  IN 


NUCLEAR  SCIENCE  AND  ENGINEERING. 


The  AFIT  1620  has  a  40K  memory,  card  input -output,  floating  di¬ 
vide,  indirect  address  and  the  usual  additional  instructions  (TNS,  TNF, 
MF.)  Our  object  decks  should  execute  without  change  on  any  similar 
machine.  It  is  our  present  plan  to  honor  requests  for  copies  of  our 
source  or  object  decks  directly.  Please  eend  such  requests  to  me  at 
the  address  below. 

Robert  H.  Hansen  carried  out  the  moments  reduction  study  as  his 
M.S.  thesis.  Ernest  Park  Sims  added  the  thermal  portion  of  the  code 
in  conjunction  with  his  M.S.  thesis  research,  and  put  the  code  in 
its  preliminary  form.  Daniel  P.  Cannon  made  the  final  revisions  and 
thoroughly  checked  the  code  out.  In  addition  he  provided  the  bulk  of 
this  revised  document.  This  effort  was  his  M.S.  thesis.  In  addition 
to  project  direction,  I  provided  the  theory  chapter. 

This  report  is  a  revision  of  an  interim  report  which  was  publi¬ 
shed  as  AFIT  TR  66-6  (May  1966.1  This  revision  contains  extensive 
code  changes  as  well  as  increased  documentation. 

Finally  I  feel  compelled  to  make  some  comment  about  the  name  OLD 
BARNYARD.  It  was  originally  coined  by  Capt.  Sims  as  a  joke;  however, 

I  find  it  is  a  name  which  the  student  remembers  and  more  importantly 
associates  with  cross  sections  and  age.  1  have  yet  to  have  a  student 
ask  "What  was  the  name  of  that  code  which  generates  cross  sections?" 


Charles  J.  3ridgman 

Associate  Professor  of  Nuclear  Engineering 
Department  of  Physics 
Air  Force  Institute  of  Technology 
Wright -Patterson  AFB,  Ohio  -  45433 
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OLD  BARNYARD 

A  CROSS  SECTION  CODE  FOR  SCHOOL  USE 

I.  Introduction 

Many  nuclear  engineering  instructors  regard  the  assignment  and 
completion  of  realistic  problems  as  vital  to  student  learning.  How¬ 
ever,  instructors  who  assign  their  students  such  problems  are  continu¬ 
ally  faced  with  the  problem  of  providing  group  neutron  cross  sections 
and  cross  sec  .on  related  constants  such  as  the  diffusion  coefficient 
and  age.  The  alternative  to  providing  such  data  is  to  require  each 
student  to  find  or  calculate  his  own  values.  This  alternative  is  un¬ 
satisfactory  for  bo  h  the  student  and  the  instructor.  For  the  stu¬ 
dent,  except  for  the  first  few  times,  finding  or  calculating  group 
cross  section  is  repetitious s  time  consuming,  and  of  little  educa¬ 
tional  value.  For  the  instructor  this  practice  is  guaranteed  to 
produce  a  considerable  variation  in  final  answers  to  evaluate. 

An  instructor  usually  can  assume  a  Maxwell-Boltzmann  distribu¬ 
tion  to  determine  thermal  group  constants  or  he  can  use  experimental 
values,  as  available.  Epithermal  and  fast  group  values  are  more 
difficult  to  obtain  however,  especially  since  only  a  few  broad  groups 
are  desired  for  learning  purposes.  Computer  codes  are  currently 
available,  such  as  GAM  (Ref  13),  to  provide  fast  cross  sections,  but 
they  usually  require  large,  fast  and  expensive  computers  not  always 
available.  Even  when  such  codes  are  locally  available,  they  usually 
employ  sophisticated  methods  which  are  sketchily  documented.  Such 
a  situation  is  acceptable  when  the  codes  are  used  by  professionals 
but  quite  unacceptable  when  the  codes  are  used  by  students. 
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As  a  result  the  instructor  ends  up  borrowing  some  epithermal  and 
fast  group  values  from  some  few  group  published  sets  such  as  those  in 
ANL  5800  (Ref  1.)  However,  these  may  not  include  all  the  elements  of 
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interest,  and  the  group  energy  boundaries  are  not  flexible. 


Statement  of  the  Problem 

A  need  exists  for  a  few-group  neutron  cross  section  code  for 
school  use  in  suppoi't  of  reactor  criticality  and  neutron  transport 
problems  which  meets  the  following  criteria:  ca’.culates  cross  sections 
and  related  constants  to  a  fair  degree  of  accuracy  (say  10%),  uses 
methods  with  which  a  student  can  be  expected  to  be  familiar,  is 
throughly  documented,  and  can  be  used  on  small  (as  well  as  large) 
digital  computers. 


Purpose  and  Method 

The  purpose  of  this  study  is  to  produce  such  a  code.  The  com¬ 
puter  code  described  in  this  report  calculates  both  fast  and  thermal 
cross  sections  and  related  constants  such  as  age  and  diffusion  length. 
It  is  written  in  Fortran  and  has  been  executed  on  an  IBM  1620  computer 
where  it  requires  about  15  minutes  for  completion.  The  cross  sections 
are  found  by  calculating  the  energy  dependent  flux  in  an  infinite 
homogeneous  mixture  of  the  isotopes  or  compounds  .specified.  The  in¬ 
put  cross  section  data  is  obtained  from  eleven  group  libraries  of  fast 
cross  sections  plus  2200  meter  per  second  thermal  values.  The  energy 
dependent  flux  is  calculated  by  the  neutron  transport,  moments  method, 
for  above  thermal  energies  and  by  a  Wilkins  equation  calculation  for 
thermal  energies. 
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The  eleven  group  input  was  selected  as  a  compromise  between 
accuracy  and  amount  of  calculation  after  a  systematic  reduction  in 
the  number  of  fine  groups  (Ref  8.)  The  above  thermal  calculations 
were  further  simplified  by  computing  resonance  absorption  with  a 
resonance  integral  taken  from  the  formulas  given  by  Murray  (Ref  16:63) 
and  Isbin  (Ref  11:459.)  The  Wilkins  spectra  is  normalized  to  the 
moments  spectra  at  1.125  ev. 

The  code  itself  is  described  in  Chapter  II,  The  Code.  This  is 
followed  by  a  chapter  on  operating  instructions  and  information  on 
preparing  input.  Thirteen  sample  problems  are  given  in  Chapter  IV. 

The  final  chapter  attempts  to  trace  the  theory  of  both  epithermal  and 
thermal  calculations  in  textbook  fashion. 
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II.  The  Code 
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consists  of  two  chains  which  must  be  run  sequentially.  The  first 
chain  computes  the  energy  dependent  flux,  both  fast  and  thermal.  This 
first  chain  also  computes  neutron  age,  resonance  parameters,  and  the 
thermal  constants.  The  output  of  this  chain  consists  of  the  group 
flux,  expressed  as  flux  per  unit  lethargy,  for  11  epithermal  groups 
plus  12  thermal  values  at  half  lethargy  intervals.  The  flux  is  nor¬ 
malized  to  a  maximum  value  of  100.  The  neutron  age  is  also  computed 
Doth  to  the  Indium  resonance  and  to  1.125  ev  (which  is  taken  here  as 
the  arbitrary  thermal  region  upper  boundary).  Should  a  resonance 
calculation  be  performed  the  scattering  per  absorber  atom,  and  effec¬ 
tive  resonance  integral  are  also  calculated. 

The  thermal  flux  values  are  computed  only  if  there  is  moderator  in 
the  system.  If  a  thermal  calculation  is  performed,  in  addition  to  the 
output  above,  the  code  computes  the  thermal  diffusion  length  squared, 
the  migration  area,  the  most  probable  and  average  thermal  neutron  ve¬ 
locities,  and  the  absorption  parameter  of  the  system.  Chain  one  also 
includes,  as  an  option,  a  plot  of  the  flux  as  a  function  of  lethargy 
over  the  entire  range  of  interest. 

The  second  chain  collapses  the  cross  sections  to  the  desired 
group  structure  using  the  fluxes  calculated  in  chain  one.  The 
operator  may  choose  between  microscopic  and  macroscopic  cross  sections 
for  chain  two  output.  A  listing  of  the  source  decks  is  contained  in 
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Appendices  B  and  D.  The  source  decks  shown  there  have  been  executed 
on  an  IBM  1620  with  40  K  memory,  card  input-output,  floating  divide, 
indirect  addressing  and  the  usual  additional  instructions  (TNS,  TNF, 
MF). 

As  a  convenience  to  the  reader,  the  more  important  quanities  cal¬ 
culated  by  the  code  are  listed  in  figure  1  along  with  the  appropriate 
equation  numbers  from  Chapter  V  of  this  report.  Figure  1  also  lists 
page  number  references  in  three  popular  texts  for  the  same  quantities 
The  texts  are  by  Glasstone  and  Edlund  (Ref  7),  Isbin  (Ref  11),  and 
Murray  (Ref  16). 


Quanity  to  be  •  Equation  Humber  in  Textbook  Reference 

Calculated  Chapter  V _  Murray  Isbin  G  £  E 


±  . 

Fast  flux  and  age 

66,  67,  68 

60 

181 

2. 

Resonance  Integral 

71  thru  74 

459 

253 

3. 

Resonance  escape 
probability 

75 

63 

454 

253 

4. 

Thermal  utilization 

78,  79 

87 

507 

264 

5. 

Absorption  para¬ 
meter 

82 

6. 

Thermal  flux 

101,  103,  105 

7. 

Most  probable  thermal 
neutron  velocity 

109 

47 

38 

3. 

Thermal  diffusion 
length  squared 

118 

66 

236 

116 

9. 

Migration  area 

122 

255 

216 

10. 

Collapsed  broad 
group  cross  sections 

125 

Figure  1 

Summary  of  Important  Quanities  Calculated 


The  input  to  chain  I  consists  of  a  cross  section  library  ("fine" 
groups)  of  11  fast  groups  and  the  2200  m/sec  thermal  cross  section. 
Related  thermal  constants  are  included  such  as  the  average  logar-imethic 
energy  decrement  C*  the  absorption  parameter  at  T  =  293  deg  Kelvin,  and 
the  quanity  (1-p)  where  \i  H  2/3A.  There  are  two  cross  section  li¬ 
braries  provided  with  the  code.  The  thermal  values  in  these  libraries 
were  obtained  from  ANL  5800  (Ref  1)  and  BNL  325  and  the  supplements  to 
it  (Ref  9).  The  eleven  fast  group  cross  sections  were  generated  by  the 
General  Atomics  68  group  moments  code,  GAM-1  (Ref  13).  Therefore,  the 
eleven  group  cross  sections ,  used  as  input ,  are  already  flux  averaged 
cross  sections.  As  such,  they  are  a  function  of  the  material  composi¬ 
tion  used  to  generate  them.  Thus  the  two  cross  section  libraries  pro¬ 
vided  represent  flux  weighting  over  both  a  thermal  and  a  fast  reactor 
core  composition.  The  Air  Force  Nuclear  Engineering  Test  Facility 
(NETF)  core,  an  MTR  type  reactor ,  was  selected  as  the  thermal  standard; 
and  the  Godiva  core,  a  bare  unmoderated  system >  was  selected  as  the 
fast  reactor  standard  for  input  to  GAM-1.  The  "thermal"  library  of  11 
fast  groups  was  generated  by  GAM  working  with  the  NETF  core  plus  a 
trace  of  the  nuclide  to  be  considered.  The  calculation  was  repeated 
using  each  nuclide  in  the  library  as  the  trace  nuclide.  The  "fast" 
library  was  generated  in  the  same  way  except  that  the  Godiva  core  was 
used  in  place  of  the  NETF  core.  The  "thermal”  library  is  intended  for 
use  with  thermal  reactor  problems;  the  fast  library  is  used  with  fast 
reactor  problems.  The  libraries  include  inelastic  (n,n*)  and  (n-2n) 
cress  sections  and  these  effects  are  included  in  the  moments  methods 
calculations.  The  libraries  also  include  fission  cross  sections  for 
use  in  chain  two.  The  user  may  elect  to  generate  his  own  library 
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using  some  other  material  composition.  The  only  restriction  is  that 
the  11  fine  group  boundaries  are  the  same  and  the  data  order  is  that 
shown  in  Figure  2. 

The  libraries  provided  with  the  code  currently  include  the 
following  nuclides : 


1. 

Hydrogen 

13. 

Nickel 

2. 

Deuterium 

14. 

Copper 

3. 

Beryllium 

15. 

Zirconium 

4. 

Boron 

16. 

Molybdenum 

5. 

Carbon 

17. 

Thorium-232 

6. 

Oxygen 

18. 

Uranium-235 

7. 

Nitrogen 

19. 

Uranium-238 

8. 

Aluminum 

20. 

Plutonium- 2 39 

9. 

Magnesium 

21. 

Fission  Products 
(of  Uranium-235) 

10. 

Chromium 

22. 

Xenon-135 

11. 

Manvaneco 

23. 

Samarium- 149 

12. 

Iron 

The  fission  products  entries  includes  all  product  nuclides  of 
thermal  Uranium-235  fission  except  Xenon-135  and  Samarium-149.  These 
two  nuclides  are  included  directly  in  the  libraries. 
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Figure  2t  Uranium-c38  Cross  oections,  Thermal  Library 


Figure  2  shows  the  U-238  cross  sections  from  the  "thermal"  library. 
The  entries  are  explained  below  as  codec  in  the  Figure. 

a)  Library  and  nuclide  identification  (Thermal  library,  U-238). 


Columns  1  thru  20 

b)  Number  of  fast  cross  sections  below  which  begins  on  the 


second  line.  Does  not  include  the  vc^.  entries  for  fission¬ 
able  nuclides. 

c)  Epithermal-fast  absorption  code,  o  =  Ho  fast  absorption 
cross  sections  provided  (assumed  negligible);  1  =  Fast 
absorption  cross  sections  provided 

d)  Type  of  inelastic  scatter:  1  =  inelastic;  2  =  (n,2n); 

3  =  both;  4  =  neither 


. - v- 


e)  Thermal  average  logarithmic  energy  decrement,  £ 

f)  Thermal  microscopic  absorption  cross  section,  og 

g)  Thermal  absorption  parameter,  A 

h)  Thermal  microscopic  scattering  cross  section,  5n  barns 

i)  Thermal  (1-p)  where  p  =  2/3A 

j )  Thermal  value  of  the  product  of  the  number  of  neutrons  per 
fission  and  the  microscopic  fission  cross  section,  vo^ 

k.)  Epithermal-fast  fission  code:  0  =  no;  1  =  yes 

l)  Microscopic,  fast,  absorption  cross  sections:  o1,  a2,*  ••a11, 
barns 

m)  Number  of  groups  from  which  neutrons  inelastically  scatter 

n)  Number  of  groups  to  which  neutrons  inelastically  scatter 

o)  Microscopic,  fast,  inelastic  cross  sections  in  the  following 
order: 


in 

in 

in 

a 

) 

a 

9 

«  «  «  0 

1-1 

1-2 

1-9 

in 

in 

in 

a 

9 

a 

9 

•  . «  o 

2-2 

• 

CO 

1 

2-10 

’in 

in 

in 

a 

) 

a 

9 

...  0 

7-7 

7-8 

7-th 

in 

in 

in 

o 

9 

a 

9 

•  •  •  0 

CD 

1 

CO 

8-9 

8-th 

p)  Number  of  groups  in  which  (n,2n)  occurs 

q)  Number  of  groups  to  which  (n,2n)  scatter  occurs 

r)  Microscopic,  fast,  (n-2n)  cross  sections  in  the  same  ox'der 
as  inelastic 

s)  Number  or  groups  from  which  neutrons  elastically  scatter 

t)  Number  of  groups  to  which  neutrons  elastically  scatter 


u)  Microscopic  Po  and  PI  elastic  scatter  cross  sections  in  the 


following  order  (See  the  note  below): 


Po 

Pi 

Po 

Pi 

;  o 

;  o 

;  o 

i-i 

1-1 

1-2 

1-2 

Po 

PI 

Po 

PI 

;  a 

;  o 

;  o 

2-2 

2-2 

2-3 

2-3 

Po  PI  Po  PI 

0  i  o  ;  o  .  ;  a 

11-11  11-11  11-th  l l -th 

NOTE:  The  microscopic  Po  elastic  scatter  cross  section  is 

the  all  angle  elastic  scatter  cross  section  o .  The  micro¬ 
scopic  PI  elastic  scatter  cross  section  is  the  product  of  g 
times  os  where  p  is  the  average  cosine  of  the  scattering 
angle.  GAM-1  outputs  the  value  3oPl  and  this  is  the  value 
used  in  the  cross  section  library.  Therefore,  in  the  code 
all  oPl  are  ai vided  by  3. 

v)  Microscopic,  fast,  va^.  in  the  order  va*,  vo|,  ...  vo*.1 

v)  Grid  background  card  for  the  spectrum  printout. 

NC".:  The  grid  card  is  considered  part  of  the  library  for 
each  nuclide. 

Since  flux  weighting  cross  sections  involves  division  by  the 
group  flux,  the  problem  of  division  by  near  zero  quantities  in  the 
lcwer  energy  groups  of  fast  reactor  systems  might  occur.  To  prevent 
any  problem  the  code  scans  the  fast  fluxes  for  a  minimum  value  of  10~12, 
A  message  is  typed  giving  the  group  number  of  the  first  fine  group  in 
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- - - - 


which  the  neutron  flux  ic  less  than  the  minimum  value.  If  the  tot:  1 


flux  in  a  broad  group  (the  sum  of  all  fine  group  fluxes  in  a  given 
broad  group)  is  less  than  the  minimum,  there  is  no  output  for  that 


—  - - o  r 

Should  the  user  desire  cross  sections  of  "thermal"  type  materials, 
such  as  carbon  or  hydrogen ,  in  the  presence  of  a  fast  flux  spectrum,  he 
can  obtain  these  simply  by  substituting  the  cross  section  library  of 
the  "thermal"  nuclide  when  executing  chain  two. 

The  special  set  of  Kingstran  subroutines*  provided  with  the  code 
must  be  used  to  permit  the  chain  operation. 


*At  AFIT,  if  new  object  decks  are  produced  (from  the  source  decks) 
using  the  Kingstran  compiler,  seven  (7)  Kingstran  trailer  cards  must  be 
added  to  the  back  of  the  chain  one  and  chain  two  object  decks  produced. 
This  is  necessary  to  use  the  special  set  01  subroutines. 


p*  .vvV(i? 


III.  Operating  Instructions:  IBM  1620 

Input  data  for  chain  one  and  chain  two  are  prepared  according  to 
instructions  contained  in  the  next  two  sections  respectively.  To 
execute  the  code  follow  the  operating  instructions  below: 

1.  Set  all  console  switches  to  PROGRAM  except  the  Parity  switch  is 
set  to  STOP. 

2.  Place  the  special  Old  Barnyard  short  subroutines  in  the  read 
hopper.  Press  RESET  and  LOAD. 

3.  When  the  last  card  of  the  subroutines  is  reached,  the  console 
READER  NO  FEED  light  will  illuminate  steadily.  Press  READER  START 
and  read  in  last  card. 

4.  Remove  the  subroutines  from  the  out  hopper.  Place  the  Old  Barn¬ 
yard  chain  one  object  deck  in  the  reader  hopper.  Place  the  data  cards 
for  chain  one  on  top  of  the  chain  one  object  deck.  Press  LOAD  and 
PUNCH  START. 

5.  Chain  one  program  will  run  automatically.  When  the  last  data  card 
enters  the  reader,  the  console  READER  NO  FEED  light  will  illuminate 
steadily.  Press  READER  START  to  read  in  last  card. 

6.  The  console  typewriter  will  indicate  the  beginning  of  each  section 
of  the  code.  After  all  calculations  to  be  executed  are  complete,  the 
typewriter  will  type  "SET  SWITCH  1  ON  FOR  SPECTRUM  PRINTOUT.  PRESS 
START."  If  a  plot  of  the  neutron  spectrum  is  desired,  set  switch  1  ON 
and  press  START.  If  no  spectrum  plot  is  desired,  set  switch  1  OFF  and 
press  START. 
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7.  At  the  end  of  chain  one  the  typewriter  will  type  "END  OF  CHAIN  1. 


TO  COLLAPSE  CROSS  SECTIONS  LOAD  CHAIN  2."  If  no  cross  sections  are 
desired,  the  program  is  finished.  Remove  chain  one  ard  data  cards  and, 
if  no  cross  sections  are  desired  the  chain  one  output  cards,  from  the 
out  hoppers . 

8.  If  cross  sections  are  desired,  place  Old  Barnyard  Chain  2  object 
deck  in  the  read  hopper.  Place  data  cards  for  chain  two  on  top  of  it. 
Press  RESET  and  LOAD.  Chain  two  will  run  automatically.  After  the 
first  data  card  is  read  in  a  typewritten  message  summarizes  the  group 
structure  requested.  When  the  last  card  enters  the  reader  the  console 
READER  NO  FEED  light  will  illuminate  steadily.  Press  READER  START  to 
reaa  in  the  last  card.  Remove  chain  two  and  data  cards  from  the  out 
hopper. 

9 .  When  chain  two  is  complete  the  typewriter  wili  type  "END  OF  PRO¬ 
GRAM".  Remove  output  cards  from  the  output  hopper.  It  is  not 
necessary  to  press  NON-PROCESS  RUN  OUT. 

10.  List  output  answer  cards  on  the  407  lister.  Put  switch  number  4 
on  the  right  side  of  the  407  UP.  This  causes  automatic  advance  of  a 
new  sheet  of  paper  for  each  section  of  output  answers.  If  the  paper 
advances  at  unwanted  times  (as  it  sometimes  does)  put  switch  4  DOWN. 

The  best  plan  is  to  put  switch  4  UP  for  chain  one  output,  DOWN  for 
chain  two  outputs. 

11.  To  run  a  second  problem  go  back  to  step  (1)  and  repeat.  NOTE: 
if  anything  should  go  wrong  during  a  run  (such  as  an  erroneous  data 
card,  or  an  I/O  error,  or  a  SKIP  CHECK),  the  particular  chain  in  use 
may  be  restarted  by  the  following  procedure : 


1)  Press  STOP. 


2)  Remove  remaining  data  cards,  if  any,  from  the  read  hopper. 

3)  Clear  the  reader  of  any  internal  date  cards  by  pressing 
NON-PROCESS  RUN  OUT. 

4)  Replace  corrected  data  cards  in  the  read  hopper. 

5)  Press  RESET  and  INSERT. 

6)  Type  on  the  console  typewriter  the  numbers  4900936. 

7)  ^ress  RESET  and  START  or  punch  the  R-S  key  on  the  typewriter. 

8)  Press  READER  START  to  read  in  data  cards. 

9)  Program  will  run. 


Preparation  of  Input  Data,  Chain  One 

The  following  data  cards  are  required  to  operate  chain  one.  Input 
data  may  be  punched  in  any  format.  If  the  E  format  is  used  for  data, 
no  blank  spaces  are  permitted  between  the  last  digit  of  the  mantissa 
and  the  E  of  the  exponent,  but  blanks  are  permitted  between  the  E  and 
the  exponent.  If  the  exponent  is  positive  the  plus  sign  is  optional. 
Individual  numbers  on  a  card  must  be  separated  by  one  or  more  blank 
spaces;  commas  are  not  allowed.  Examples  of  correct  and  incorrect 
methods  are  shown  below. 


Correct 

1  1.5 

1.67E+01 
1.67E  01 


Incorrect 

1,1.5  or  1,  1.5 

1.67  E+01 


Card  No .  1 : 

This  card  contains  up  to  65  columns  of  alphameric  data  of  the 
User's  choosing.  Identification,  problem  title,  date,  etc.,  are 


possible  information  to  be  placed  on  this  card.  The  data  on  this  card 
will  be  typed  out  on  the  console  typewriter  and  will  be  punched  at  the 
head  of  the  output  data. 

EXAMPLE:  AFIT  Subcritical  Core.  5  April  67  Mackrill 

Card  No.  2: 

This  card  has  six  (5)  numbers  as  follows; 

a)  A  number  between  two  (2)  and  twenty-one  (21)  indicating  the 
source  of  neutrons  the  user  wants  to  use.  The  meaning  of  the  numbers 
is  shown  below: 


Number 

2 

3 

a 

5 

5 

7 

8 
9 

10 

11-21 


Source 

U235  fission  (Cranberg  Spectrum) 

U233  fission 
Pu239 

Pu24l 

cf2S2 

Pu-Be  (Whitmore-Baker) 

Pu-Be  (Cochrane--Henry ) 

Ra-Be  (Hill) 

The  user's  own  source  (see  Card  No.  4) 

A  unit  source  in  group  1-11  respectively 


EXAMPLE:  15  =  A  unit  source  in  group  5. 

b)  A  number  giving  the  number  of  nuclides  to  be  used  in  the 


problem 


EXAMPLE : 


2  for  Water  (H  and  0) 


2  for  BeO  (Be  and  0) 

2  for  C 2 ■  i g  (0  and  H) 

3  ror  an  mixture  (Ai,  0,  and  H) 

1  for  pure  carbon. 

c)  A  number  giving  the  total  nuclei  density  for  the  problem  in 
atoms  times  10"24  per  cm3,  (atoms  per  barn-cm) 

EXAMPLE:  For  Carton  0.08025 

For  Water  (HO)  0.1005 

d)  A  number  giving  the  temperature  of  the  system  in  degrees 
Kelvin. 

e)  A  number  indicating  whether  or  not  a  resonance  escape  calcu¬ 
lation  should  be  performed;  as  below: 

0  -  NO 
1  =  YES 

f)  A  number  indicating  whether  or  not  the  mixture  is  considered 
homogeneous  or  heterogeneous;  as  below: 

0  =  homogeneous 
1  =  heterogeneous 


Card  No.  3: 

This  "card"  is  a  small  deck  of  cards  containing  the  neutron 
source  data  for  the  program.  It  is  included  as  part  of  the  input 
library  and  normally  requires  no  preparation  on  the  part  of  the  user. 
The  source  deck  must  always  be  included  even  when  the  user  desires  to 
read  in  his  own  source  spectrum.  A  listing  of  the  sources  is  shown  in 
Appendix  E. 
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Card  No.  4: 


This  card  is  required  only  if  the  first  number  on  card  no.  2  is 
10;  i.e.,  only  if  the  user  wants  to  read  in  and  use  his  own  source. 
This  card  contains  eleven  numbers  representing  the  fraction  of  the 
source  in  each  of  the  eleven  epithermal-fast  groups. 

Card  No.  5: 

This  card  is  included  only  if  the  fifth  number  on  card  no.  2  is 
1;  i.e.,  if  a  resonance  escape  calculation  is  to  be  performed.  This 
card  contains  three  (3)  numbers  as  follows: 

a)  A  number  indicating  whether  the  resonance  nuclide  is  U238, 
Th232  or  W  (Tungsten)  as  shown  below: 


1  =  U  or  W 

2  =  Th 


b)  A  number  giving  the  nuclei  density  of  the  resonance  nuclide(s) 
in  atoms  per  barn-cm. 

c)  A  number  indicating  whether  the  effective  resonance  integral 
is  being  supplied  by  the  user,  or  whether  it  is  to  be  calculated  by  the 
program : 


0  =  to  be  calculated 


non-zero  =  value  of  RI  «  to  be  used 

eff 


Card  No.  6: 

This  card  is  included  only  if  the  sixth  number  on  card  no.  2  is 
1;  i.e.,  only  if  the  mixture  is  considered  heterogeneous.  This  card 
contains  two  (2)  numbers  giving  the  "F"  and  "£"  factors,  in  that  order, 
for  calculating  the  effect  of  heterogeneous  cell  structure  on  thermal 


17 


utilization.  (See  the  section  on  thermal  utilization  for  the  equations 
used  to  calculate  the  "F"  and  "E"  factors). 

Cards  number  7  and  8  are  required  for  each  nuclide  in  the  sequence 

7,  8,  7,  8,  ...  . 

Card  No.  7: 

This  card  contains  six  numbers  as  follows: 

a)  Twenty  (20)  columns  of  alphameric  identification  to  identify 
the  nuclide  of  the  problem.  This  is  optional,  for  the  user's  use 
only,  and  may  be  left  blank  (the  next  number  must  always  start  after 
column  20  however). 

b)  A  number  giving  the  nuclei  density  of  the  ntfl  nuclide  in  the 
problem  in  atoms  per  barn-cm. 

c)  A  number  indicating  whether  the  nuclei  is  a  moderator,  fuel, 
or  other  as  below: 

1  =  moderator 

2  =  fuel 

3  =  other 

d)  A  number  indicating  whether  microscopic  or  macroscopic  cross 
sections  are  to  be  output  for  this  n  nuclide.  A  number  is  always 
required  here  even  if  cross  section  averaging  is  not  to  be  performed 
(chain  two  of  the  program).  The  number  should  be: 

1  =  microscopic 

2  =  macroscopic 

e)  A  number  indicating  if  this  is  a  resonant  nuclide  (U238, 

Th232  or  W).  The  number  should  be: 

0  =  NO 

1  =  YES 
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Card  No.  8: 


This  "card"  is  a  deck  of  cards  containing  the  eleven  (11)  group 
plus  thermal  cross  sections  of  the  n  ‘  nuclide.  This  deck  is  selected 
from  one  of  the  two  program  libraries  provided,  and  requires  no  prepa¬ 
ration  on  the  part  of  the  user.  The  thermal  library  is  for  use  with 
thermal  reactor  problems ,  and  the  fast  library  is  for  use  with  fast 
reactor  problems.  The  Uranium-238  cross  section  deck  from  the  thermal 
library  is  shown  in  Figure  2.  The  grid  background  card  (the  last  card 
of  each  cross  section  deck)  must  be  included  for  each  nuclide. 

Figure  3  shows  the  arrangement  of  the  entire  input  deck  for  chain 
one.  The  following  are  some  notes  on  chain  one  data: 

a)  The  total  nuclei  density  on  the  first  card  should  be  the  sum 
of  the  individual  nuclide  densities  on  the  separate  card  no.  7's. 

b)  The  nuclei  density  of  the  resonance  nuclide  on  card  no.  3 
should  agree  with  that  on  card  no  7  for  the  resonance  absorber. 

Preparation  of  Input  Data,  Chain  Two 

Chain  two  of  the  program  requires  the  preparation  of  only  one 
additional  card.  This  card  has  up  to  twelve  (12)  numbers.  The  first 
is  the  number  of  epithermal- fast  broad  groups  for  which  cross  sections 
are  to  be  calculated.  The  remaining  numbers  are  the  numbers  of  the 
lowest  fine  group  in  each  of  the  selected  broad  groups.  The  last 
number  on  this  card  is  eleven  (11)  indicating  that  the  last  epithermal- 
fast  broad  group  extends  down  to  .4\4  ev.  Thus  if  the  first  number  is 
five  there  will  be  five  more  numbers  on  the  card.  If  the  first  number 
is  eleven,  there  will  be  eleven  more  numbers  on  the  card. 
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EXAMPLES: 


1  11 
2  6  11 

11  123456789  10  11 

This  additional  card  is  the  first  input  data  card  for  chain  two. 
It  is  followed  by  the  collection  of  cards  7  and  8  from  chair,  one  as 
shown  in  Figure  4. 


dame  for  nuc^iae  #3 

dame  for  nuclide  #2 

Nuc.den.  and  cross 
sect. ,nuclide^l 

Heterogeneous  cell 
factors  (optional) 

xiesonance  data  (optional) 

Users  own  source  (optional) 

“  dource  deck(supplied) 

'dource  number,  no.  of  nuclides, total  nuclei  density, 

T,  type,  media 

-Problem  li),  title,  users  name,  date,  etc. 


Figure  3;  Chain  One  Input 


Figure  4:  Ohain  Two  Input 


IV.  Sample  Problems 


In  the  proceeding  chapter,  instructions  were  given  showing  how  to 
prepare  input  for  the  code  anu  how  to  run  the  code  using  the  IBM  ]G?0 
digital  computer.  This  chapter  is  devoted  to  exhibiting  "OLD  BARNYARD" 
output  for  selected  problems.  Further  explanation  of  some  of  the  prob¬ 
lems  is  given  below.  The  following  is  a  list  of  the  sample  problems 
included  in  this  chapter: 


i-.y 


I  • 


(1)  Pure  Carbon 

a)  1  fast  group,  T  =  293  deg.  Kelvin 

b)  2  fast  groups,  T  =  293  deg  Kelvin 

c)  1  fast  group,  T  =  1273  deg.  Kelvin 


(2)  Water 


‘V 


a)  1  fast  group 

b)  2  fast  groups 

(3)  Pure  Uranium-235 

a)  1  fast  group 

b)  3  fast  groups 

(4)  Air  Force  Nuclear  Engineering  Test  Facility  (NETF)  core 

a)  1  fast  group 

b)  2  fast  groups 

(5)  Air  Force  Institute  of  Technology  (AFIT)  Subcritical 


< 


a)  1  fast  group  - 

b)  2  ^art  groups 

( 6 )  Romaska  Core 


?  ■ 
r. 


a)  1  fast  group 

b)  3  fast  groups. 

The  first  two  problems  listed  above  are  seif  explanatory.  Prob¬ 
lem  three,  pure  Uranium-235  is  shown  in  figures  10  and  11  for  one  and 
three  fast  groups  respectively.  This  is,  of  course,  a  fast  system 
with  .10  moderator  present.  The  reader  will  note  that  in  this  case 
thermal  calculations  are  not  performed,  as  was  explained  earlier  in 
Chapter  II.  Therefore,  only  ten  epithermal  group  fluxes  are  output 
plus  the  age  to  the  Indium  resonance  in  chain  one.  Also,  the  flux 
plot  is  terminated  at  lethargy  16  in  this  case. 

The  Air  force  Nuclear  Engineering  lest  Facility  (NETF)  is  a  ten 
megawatt,  approximately  90%  enriched,  unpressurized,  thermal  research 
reactor.  The  core  has  a  rectangular  parallelepiped  arrangement  of  MTR 
flat-plate  type  fuel  elements,  cooled  and  moderated  by  light  water. 

The  Air  Force  Institute  of  Technology  (AFIT)  subcritical  reactor 
is  a  natural  uranium  reactor  with  cylindrical,  aluminum  encased,  fuel 
elements  arranged  in  a  cylindrical  pattern.  The  core  is  moderated 
with  light  water. 

Finally,  the  Romaska  reactor  is  a  Russian  direct  energy  converter 
reactor  with  a  UC2 ,  C  core.  The  Romaska  core  has  a  Carbon-Uranium 
atom  ratio  of  approximately  6  tc  1.  Thus,  the  Romaska  core  is  a  fast 
reactor  system.  Output  for  this  core  is  shown  in  figures  16  and  17 
for  one  and  three  fast  groups  respectively.  Since  there  is  some  mod¬ 
erator  in  the  system,  thermal  calculations  are  performed  even  though 
the  thermal  flux  is  very  low. 
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JOB  OLD  BARNYARD.  CHAIN  1 


CARBON  1  FAST  GROUP  21  FEB,  1967 


THE  SOURCE  USED  IS  U235  FISSION  (CRANBERG  SPECTRUM). 
SYSTEM  TEMPERATURE  IS  293.0  KELVIN 
ATOM  FRACTIONS  ARE 

CARBON  1.00000000 


RELATIVE  GROUP 


GROUP 

LETHARGY 

ENERGY,  E V 

SOURCE 

FLUX,  /UNIT 

1 

.50 

6.0653E+06 

2.3023E-02 

5.3003E-02 

2 

1.00 

3.6788E+0& 

1.0824E-01 

2  «  8939E-0 1 

3 

1.50 

2.2313E+06 

2, 1044E-01 

7.6130E-01 

4 

2.00 

1.3534E+06 

2. 3139E-01 

1 . 3565E+00 

5 

2.50 

8.2085C+05 

1 . 8048E-01 

1 • 3975E+00 

6 

3.00 

4.9787E+05 

1. 1483F-0’ 

1 .4450E+00 

7 

5.00 

6. 7379E+04 

1.2439E-01 

1 . 3667E+00 

8 

8.00 

3.3546E+03 

1.2094E-03 

1.1381E+00 

9 

12.00 

6.  1442E-*-01 

• OOOOE+OO 

1 .0195E+00 

10 

16.00 

1.  1254E+00 

.0000E+00 

1.045^E+00 

11 

17.00 

4. 1399E-01 

. 0000E+00 

1.1 5 16E+00 

12 

17.50 

2.51 10E-0 1 

• OOOOE+OO 

2. 1834E+00 

13 

18.00 

1 . 5230E-01 

• OOOOE+OO 

1 . 7762E+0 1 

14 

18.50 

9.2374E-02 

. OOOOE+OO 

6.5395E+01 

15 

19.00 

5.6028E-02 

. 0000E+00 

1 . 0000E-+  02 

16 

19.50 

3  •  3983E-02 

.OOOOE+OO 

8.7568E+C1 

17 

20.00 

2.0612E-02 

.OOOOE+OO 

5.4531E+01 

18 

20.50 

1  •  2502E-02 

.OOOOE+OO 

2.7603E+01 

19 

21.00 

7  •  5826E-03 

.OOOOE+OO 

1 .232 1E+01 

20 

21.50 

4.5991E-03 

*  OOOOE+OO 

5.0958E+00 

21 

22.00 

2.7895E-03 

.OOOOE+OO 

2.0123E+00 

22 

22.50 

1.6919E-03 

.OOOOE+OO 

7.7271E-01 

23 

23.00 

1 . 0262E-03 

.OOOOE+OO 

2.9172E-01 

AGE  TO  INDIUM  RESONANCE  (1.46EV)  IS 
AGE  TO  ARBITRARY  THERMAL  (1.12EV)  IS 
THERMAL  DIFFUSION  LENGTH  SQUARED  IS 
TOTAL  MIGRATION  AREA  IS 

MOST  PROBABLE  THERMAL  NEUTRON  VELOCITY  IS 
AVERAGE  THERMAL  NEUTRON  VELOCITY  IS 
ABSORPTION  PARAMETER  IS 


3.2282E+02  CM2 
3.2706E+02  CM2 
3.0762E+03  CM2 
3.4032E+03  CM2 
2 . 5493E+03  M/SEC 
2.5200E+03  M/SEC 
2.4099E-02 


figure  5a:  Carbon,  1  Fast  Group,  Chain  1 


PLOT  OF  SPECTRUM,  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY,  FOLLOWS 


JOB  <UD  BARNYARD  CHAIN  2Z 


THE  TOTAL  NUMBER  OF  BROAD  GROUPS  !S  2 
OUTPUT  WILt  BE  FOR  2  BROAD  GROUPS 
THE  BOUNDARY  FINE  GROUPS  ARE  11  12 


NUCLIDE  IS  CARBON 

ITS  NI'MBER  DENS!  tv  ;c  S.0230E-02  PtR  BARN-CM 


MICROSCOPiC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 


SCATTER 

elastic 

ELASTIC 

total 

FROM  TO 

INELASTIC 

N-2N 

t  PO) 

(PI) 

transfer 

1  1 

r 00068 

.00000 

4.21159 

.28014 

4.21227 

1  2 

.00000 

.00000 

.04126 

-.01204 

.04126 

2  2 

.00000 

.00000 

4.80000 

.26688 

4.80000 

GROUP 

3IGTR 

SIGR 

SIGA  NUSI6F 

SOURCE 

1 

3  •  9  f  5  A  3 

04  1 26 

00000  .00000 

1.00000 

2 

4.53661 

00349 

00549  .00000 

MAXWELL- 

BOLTZMAN  FACTOR 

=  1.128, 

AVERAGE  X  =  1 • 1468881 F  +00 

Figure 

5c:  Carbon,  1 

Fast  Group, 

Chain  2 

JOB  OLD  BARNYARD. 

CHAIN  1 

CARBON 

1 

Fast  GROUP 

S  21  FFB,  1967 

THF 

c OIJRCF  USED  IS  IJ235 

F  I  sc,  I  ON  'CRANBERG  SPECTRUM!. 

SYSTfM  TEMPERATURE  IS 

293.0  LEI  VIN 

ATOM 

FRACTION  ARE 

CARBON 

1  .00000000 

RELATIVE  GROUP 

GROUP 

LETHARGY 

ENERGY,  EV 

SOURCE 

FLUX,  /UNIT  U 

1 

.SC 

6 . 065  3E+06 

2. 3023E-02 

5.3003E-02 

2 

1.00 

3.6788E+06 

1.0824E-01 

2.8939t:-0l 

3 

1.50 

2.23 1 3E+06 

2. 1044E-01 

7.6130  01 

A 

2.00 

1.353AE+06 

2. 3139E-01 

1 . 3565E+00 

5 

2.50 

8. 20856+05 

1.804BE-01 

1.39  /5E+00 

6 

3.00 

A.9787E+05 

1.  1A83E-01 

1 . 44506+00 

7 

5.00 

6.  7379E  +  0A 

1.2A39E-01 

1 . 3667E+00 

8 

8.00 

3 • 35466+03 

7.2094E-03 

1  • 1 3  8 i E  +00 

9 

12.00 

6. 1A42E+Q1 

.oooot+oo 

1.0195E+00 

10 

16.00 

1 . 1254E+00 

. OOOOE  +00 

1 . 04546+00 

11 

17.00 

A. 13996-01 

.OOOOE+OO 

1 . 1516E+00 

12 

17.50 

2.51 10E-01 

. 0000E+00 

2. 1834E+0C 

13 

3  8.00 

1.5230E-01 

.OOOOE  +  00 

1 .7762E+01 

1 A 

18.50 

9.2374E-G2 

.OOOOE+OO 

6.5395E+01 

15 

19.00 

5.6028E-02 

.OOOOE+OO 

l.OOCOE+02 

16 

19.50 

3.3983E-02 

.OOOOE+OO 

8 . 7568E+0 1 

17 

20.00 

2.0612E-02 

.OOOOE+OO 

5.4531E+01 

18 

20.50 

1.2502E-02 

.OOOOE  i-OO 

2.7603E+01 

19 

21.00 

7.5826E-03 

.OOOOE+OO 

1.2321E+01 

20 

21.50 

a. 5991E-03 

.OOOOE+OO 

5.0958E+00 

21 

22.00 

2.7895E-03 

.OOOOE+OO 

2.0I23E+00 

22 

22.50 

1 .69196-03 

.OOOOE+OO 

7.7271E-01 

23 

23.00 

1.0262E-03 

.00006+00 

2.9172E-01 

AGE  TO 

INDIUM  RESONANCE  (1.A6EV) 

IS  3.2232E+02  CM2 

AGE  TO 

ARBITRARY  THERMAL  (1.12EV) 

IS  3.2706E 

+02  CM2  ] 

THERMAL 

DIFFUSION  LENGTH  SQUARED 

IS  3.0762E+03  CM2 

TOTAL  MIGRATION  AREA  IS 

3.4Q32E+03  CM2 

MOST  PRC9ABLE  THERMAL  NEUTRON  VELOCITY  IS  2.5A93E+03  M/SEC  I 

AVERAGE 

THERMAL  NEUTRON  VEL0CITV 

IS  2.5200E+03  M/SEC 

ABSORRT 

ION  PARAMETER  IS 

2. 4Q99E 

-02 

1 

figure  6a; 

Carbon,  2 

fast  Groups,  Chain  1 

PLOT  OP  SPECTRUM,  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY,  FOLLOWS 


JOB  OLi>  BARNYARD  CHAIN  2Z 

T Hr  TOTAL  N'JMFER  Of-'  BROAD  GROUPS  IS  3 
miTDiir  urILL  BE  rOP  3  BROAD  GROUPS 
THf  ROUND AR v  -INF  GROUPS  ARE  6  11  12 

NUCLIDE  IS  CARBON 

ITS  NUMBER  DENSITY  IS  8.0230E-02  PER  8ARN-CM 

MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS i  FOLLOW. 


SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N 

I  PO) 

(  PI  ) 

transfer 

1 

1 

. 00403 

.00000 

2.14892 

.32287 

2.15295 

1 

2 

.00061 

.00000 

. 24404 

-.06270 

.24465 

1 

3 

,00000 

.00000 

.00000 

.00000 

.00000 

2 

? 

.00000 

.00000 

4.52366 

.28357 

4.52366 

2 

3 

.00000 

.00000 

.04834 

-.0  1 A  10 

.04834 

3 

3 

.  00000 

.00000 

4,80000 

.26688 

4.80000 

GROUP 

SIG7R  SIGR 

SIGA  NUSIGF 

SOURCE 

l 

2 . 1 3  7  A  3  .24465 

• 

0C000  .00000 

.86840 

2 

4.30254  .04834 

• 

00000  .00000 

.13160 

3 

4.53661  .00340 

« 

00349  .00000 

MAXWEL 

L- 

80LTZV3N  FACTOR  =  1. 

128, 

AVERAGE  X  =  1. 1468881E+00 

figure  6c:  Carbon,  2  fast  Groups,  Chain  2 


r 

jor>  oio  rarnvarp 

.  CHAIN  1 

CARBON  AT  ]OO0  r'fG 

i 

C  f  N  T  I  CwADr  1  t 

AST 

GROUP 

THF 

i 

AOt  iRCc  USED  ! a  "23 

5  F  T  c  A  I  0 N  ('RANPERG 

SPECTRUM  )  . 

A  Y  S  T  r  M  TEMPERATURE  !S 

127  <.0  < L L V 1 N 

A  T  •'(V 

FRACTIONS  ARE 

CARBON 

1.00000000 

CROIJ 

RELATIVE  GROUP 

o  LETHARGY 

ENERGY,  EV 

SOURCE 

flux,  /unit  u 

1 

.50 

6 . 065  3  E  +  06 

2.3023E-02 

6.641 OE -02 

2 

1  .oc 

3.6788E+06 

1.0824E-01 

3.6259E-01 

3 

1.50 

2.2313E+06 

2  « 1 OA  4  E-0 1 

9. 5386E-0 1 

A 

2.00 

1 . 3534E+Q6 

2.3139E-01 

1 .6996E+00 

5 

2.50 

8. 2035E+05 

1.80-8E-6 

1  .  7 509E  +00 

6 

3.00 

4.978  7E  +  05 

1. 1483E-0, 

1 . 3 1 0  5  E  +00 

I  7 

5.00 

6. 73 7 /E +04 

1 .2439E-01 

1.7 1 24E+00 

8 

8.00 

3.3546E+03 

7. 2094F-03 

1 . 4260E  +00 

9 

12.00 

6. 1442E---01 

. OOOOE+OO 

1.2  773E  +  00 

10 

16.00 

1 . 1 254E+00 

. OOOOE+OC 

1 .3098E+00 

11 

17.00 

4.1 399E-01 

.OOOOE+OO 

6. 1758E-»0i 

12 

17.50 

2.51 10E-01 

.OOOOE+OO 

9.954CE+01 

13 

18.00 

1.5230E-C1 

. OOOOE+OO 

8.9865E-01 

14 

18.50 

9.2374E-0 2 

.OOOOE'-OO 

5 , 6999E-*-0  1 

15 

19.00 

5.6C28E-02 

.OOOOE+OC 

2 • 9 1 77E+0 1 

16 

19.50 

3.3983E-02 

.OOOOE+OO 

1.3114E+01 

17 

20.00 

2.0612E-0 2 

.OOOOE+OO 

5  •  4468E  +  00 

18 

20.50 

1 .25C2E-02 

.OOOOE+OO 

2.1  567E+0C 

19 

21.00 

7.5S26E-03 

.OOOOE+OO 

8.2955E-01 

20 

21.50 

4. 599 1 f -  03 

.OOOOE+OO 

3. 1352E-01 

21 

22 <  00 

2.7895E-03 

.OOOOE+OO 

1.1724E-01 

22 

22.50 

1.6919E-03 

.OOOOE+OO 

4.3556E-02 

23 

23.00 

1 .0262E-C3 

.OOOOE-t  00 

1.6119E-02 

AGE  TO  INDIUM  RESONANCE  (1.46EV) 

I  s 

3.2A76E+02  CM2 

AGE  TO  ARBITRARY 

THERMAL  (1.12EV) 

IS 

3.2902E+02  CM2 

thermal  diffusion 

LENGTH  SQUARED 

IS 

6  . 3847E  +  03  CM2 

TOTAL  MIGRATION  AREA  IS 

( .7137E+03  CM2 

MOST  PROBABLE  THERMAL  NEUTRON  VELOCITY  IS  5.3968E+03  M/SEC 

AVERAGE  THERMAL  NEUTRON  VELOCITY 

I  S 

5.1951E+03  M/SEC 

i 

ABSORPTION  PARAMETER  IS 

1.1 561E-02 

.Figure  7a:  Carbon, 

1  .Fast  Group,  1000 

Deg,  KeiVin, 

Chain  1 

JOB  OLD  BARNYARD  CHAIN  2Z 

THE  TOTAL  NUMBER  OF  8RC\D  GROUPS  IS  2 
OUTPUT  WILL  BE  FOR  2  3ROAD  GROUPS 
THF  BOUNDARY  FINE  GROUPS  ARE  11  12 


NUCLIDE  IS  CARBON 

ITS  NUMBER  DENSITY  IS  7.9990E-02  PER  BARN-CM 


MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 


SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N  (PO) 

(Pi  > 

transfer 

1 

1 

.00068 

00000  4.21159 

.28014 

4.21227 

1 

2 

.00000  . 

00000  .04126 

-.01204 

.04126 

2 

2 

.00000  . 

00000  4.80000 

. 26688 

4.80000 

GROUP 

SIGTR  SIGR 

£  IGA  NUSIGF 

SOURCE 

1 

3.98543  .04126 

.00000  .00000 

1.00000 

2 

4.53665  .00353 

.00353  .00000 

MAXWELL-BOLTZMAN  FACTOR  *  1.128*  AVERAGE  X  ■  1. 1343391E+00 


Figure  7ot  Carbon,  1  Fast  Croup,  1000  Deg.  Kelvin,  Chain  2 


JOB  OLD  BARNYARD  *  CHAIN  i  | 

WA  ER  1  »'AST  CROUP  21  FEB,  19 6  7 


THF  SOURCE  used  IS 

U235  r  I  SSI  ON  (CRANBERG  SPECTR'JM). 

N  v  S  T  r  m  t  r  MPgR  a  TL'RP 

IS  29  3.0  K t  L  V  !  N 

ATOM  FRACTIONS  ARP 

HvDROG£X 

,  66660666 

OXYGEN 

.  33333333 

RELATIVE  GRQ. 

GROUP 

Lfc'  THARGY 

ENERGY,  £V 

SOURCE 

Flux,  /unit 

1 

.SO 

6.0653E+C9 

2, 3023E-02 

5 .6299c -0 1 

2 

1.00 

3.6788E+06 

l.O62^E-0i 

2. 1 18  ’E-00 

3 

1.30 

2.2313E+06 

2 . 1044E~Gl 

4 . 1 6O0E  +00 

A 

2.00 

1 . 353A-E  +  06 

2.3139E-01 

4. 1480E  +00 

u 

2.30 

8.2085E  05 

1 .8048E-01 

3 . 5  71 2E+00 

6 

3.00 

4.9787E+05 

1. 1483E-01 

3 , 00B6E+00 

7 

5.00 

fa.  73  79E  +04 

1.2439E-01 

1 .6096E+C0 

8 

8.00 

3.3546E-*-03 

7. 2094 C -03 

S.5745E-C1 

9 

12.00 

6. 1442E+G1 

. OOOOE+OO 

7 • 6440E-Q 1 

10 

16.00 

1  .  1254E-.00 

•OOOOE+OO 

7.5715E-01 

1 1 

17.00 

4.13  99  E  —  0 1 

.0000E+00 

8 . 3546E-0 1 

12 

17.50 

2. 5 1 10E-01 

•OOOOE+OO 

1.8262E+00 

13 

18.00 

1.5230E-G1 

.00002+00 

1.7363E+01 

14 

18.50 

9.2374F-P2 

. OOOCE+OO 

6.5127E+C1 

IS 

19.00 

5. 6023E-02 

.OOOOE+OO 

1.0000E+02 

16 

19.50 

3 • 3983E-02 

.OCOOE+OO 

8.7712E+01 

17 

20.00 

2.0612E-02 

.OOOOE+OO 

5 . 46? 1 E+0 1 

18 

20.50 

1 • 2502E-C2 

.OOOCE+OO 

2.7  690E+0 1 

19 

21.00 

7. 5826E-03 

.OOOOE+OO 

1 . 2365E+0 1 

20 

2  1. 50 

-..5991E-03 

.OOOOE+OO 

5*11 54E+00 

21 

22.00 

2 . 78°5E-03 

.OOOOE+OO 

2.0205E+00 

22 

22.50 

1.69 192-03 

.OOOOE+OO 

7 . 7599E-0 1 

2  3 

23.00 

1.0262E-03 

.OOOOE+OO 

2o9300E-0l 

AGE  TO 

INDIUM  RESONANCE  (1.46EV) 

IS 

2 . 5224P+0 1  CM2 

AGE  TO 

ARBITRARY  THERMAL  (1.12EV) 

IS 

2.5338E 

+01  CM2 

THERMAL 

DIRFUS! 

ON  LENGTH  SQUARED 

IS 

6.2752E+00  CM2 

TOTAL  MIGRATION 

AREA 

IS 

3.1614E+01  CM2 

MOST  PROBABLE  THERMAL 

NEUTRON  VELOCITY 

IS  2.5493E 

+03  M/SEC 

AVERAGE 

thermal 

NEUTRON  VELOCI Tv 

I  S 

2.51 11E+03  M/SEC 

ABSORPTION  PARAMETER 

IS 

1.7454E 

-02 

Figure  8a;  Water,  1  Fast  Group,  Grain  l 


33 


PLOT  OF  SPECTRUM,  iN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY,  FOLLuwS 


JOB  OLD  BARNYARD  CHAIN  27 


THF  TDTAI  MUMP ; R  nc  R  R  fl  A  H  r,PnijP$  IF  ? 
OUTPUT  WILL  SC  FOR  2  BROAD  GROUPS 
THE  BOUNDARY  FINE  GROUPS  ARE  11  12 


NUCLIDE  IS  hydrogen 


ITS  ►■UMBEP  DENSITY  IS  6.680CE-02  PER  BARN-CM 


MICROSCOPIC 

CROSS  SEC 

7  IONS, 

IN  BARNS,  FOLLOW. 

SCATTER 

ELASTIC 

ELASTIC 

FROM 

TO 

INELAST 

IC 

N-2N 

(  PO) 

(PI  > 

1 

1 

.oooor 

.00000 

10.63433  7 

.21949 

l 

2 

.00000 

.00000 

.70545 

.31223 

2 

2 

.00000 

.00000 

38.00000  25 

.13320 

GROUP 

SIGTR 

SIGR 

SIGA  NUSIGF 

SOURCE 

1 

3. 

81308 

.71048 

00502  .00000 

1.00000 

2 

13. 

15730 

.29050 

• 

29050  .00000 

NUCLIDE  IS  OXYGEN 

ITS  NUMBER  DENSITY  IS  3.3400E-02  PER  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS,  FOLLOW. 


SCATTER 

FROM  TO 

INELASTIC 

N-2N 

ELASTIC 

(PO) 

ELASTIC 
(PI  ) 

1 

1 

.00118 

.00000 

3.34723 

.25328 

1 

2 

,00000 

.00000 

.01400 

-.00420 

2 

2 

.00000 

.00000 

4.20000 

.17514 

GROUP 

SIGTR 

SIGR 

SIGA 

NUS I Gr 

SOURCE 

1 

3. 11941 

.02009 

.00609 

»00000 

1.00000 

2 

4.02661 

.00175 

.00175 

« GOOOO 

MAXWELL-BOL  T7.MAN  FACTOR  ^  1,123,  AVERAGE  X  =  1  .  1428491E+00 


Figure  8e:  Water,  1  Fast  Group,  Chain  2 
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total 

TRANSFER 

10.63433 

.70545 

38.00000 


TOTAL 

TRANSFER 

3.34840 

.01400 

4.20000 


s ss<2&*iS9 


UTtvr 

•  r  * 


JOB  OLD  barnyard.  CHAIN  1 

WATER  2  FAST  GROUPS  21  FEB.  1967 

1 

THE 

SOURCE  USED  IS 

U235  FISSION  (CRANBERG 

SPEC  I KUM i . 

! 

SYSTEM  TEMPERATURE 

IS  293.0  <EL V i N 

ATOM 

FRACTIONS  are 

hydrogen 

.66666666 

OXYGEN 

.33333333 

RELATIVE  GROUP 

GROUP  letharc 

v  ENERGY.  EV 

SOURCE 

FlUX.  /UNI T  u 

1 

.50 

6 . C65  3E  *06 

2.3023E-02 

5 • 6299E-0 ! 

2 

1.00 

3.6788E+06 

1.0824E-01 

2.1 187E+00 

a 

1.5C 

2.231 3Ea06 

2. 1044E-01 

4. 1606E+00 

4 

2.00 

1.3:>34£-i06 

2. 3139E-01 

4. 1480E+00 

5 

2.50 

8.2085E+C5 

1.8048E-01 

3.5712E+00 

6 

3.00 

4.9787E+05 

1. 1483E-01 

3 . 0086E+00 

7 

5.00 

6. 7379E+0A 

1.2439E-01 

1.6096E+00 

8 

8.00 

3. 3546E+03 

7.2094E-03 

8.574SC-01 

9 

12.00 

6. 1442E-*-Ql 

. OOOOE  +00 

7.6C40E-01 

10 

16.00 

1 . I254E+00 

. 0000E+00 

7.5715E-01 

11 

17.00 

4, 1399E-01 

.OOOOE+OO 

8 . 3546E-0 1 

12 

17.50 

2.51 i0£-0i 

.0000E+00 

1.8262E+00 

13 

18. OC 

1 • 5230E-Q1 

.OOOOE+OO 

1 • 7363E+0 i 

14 

11.50 

9.2374E-02 

.0000E->-00 

6 . 5 127E+0 1 

15 

19.00 

5.6028E-02 

.OOOOE+OO 

1.0000E+02 

16 

19.50 

3.3983E-02 

.OOOOE+OO 

8.7712E+01 

1? 

20.00 

2.0612E-02 

.OOOOE+OO 

5.4671E+C1 

ie 

20.50 

1.2502E-02 

.OOOOE+OO 

2.7  690E  +  0 1 

19 

21.00 

7.5826E-03 

.OOOOE+OO 

1.2365E+01 

20 

21.50 

4.5991E-03 

»0000t+00 

5.1 154EJ-00 

21 

22.00 

2.7895E-03 

.OOOOE+OO 

2.0205E+00 

22 

22.50 

1.6919E-03 

.OOOOE+OO 

7.7  599E-0 1 

23 

23.00 

1 .0262E-03 

.OOOOE+OO 

2.9300E-01 

AGE  TO  INDIUM 

RESONANCE  U.46EV) 

I  S 

2.5224E+01  CM2 

AGE  TO  ARBITRA 

SY  THERMAL  C 1.126V) 

IS 

2.5335E 

+01  CM2 

thermal  diffus 

ION  length  SGUAREO 

I  S 

6.2752E+00  CM2 

total  migration  area  is 

3.1 6 14E+0 1  CM2 

MOST  PROBABLE 

THERMAL  NEUTRON  VELOCITY  IS  2.5493E+03  M/5E0 

AVERAGE  THERMAL  NEUTRON  VELOCITY 

IS 

2.5111E+03  M/SEC  i 

1 

absorption  parameter  IS 

1.7454E 

-02 

Figure  9a:  Water, 

2  Fast  Groups,  Chain  i 
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PLOT  OF  SPECTRUM,  5  N  ci_UX  PER  UNIT  LETHARGY  VERSUS  LETHARGY*  FOLLOWS 


PLOT  OF  SPECTRUM,  IN  clUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY*  FOLLOWS 
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JOB  OLD  BARNYARD  CHAIN  21 

THE  TOTAL  NUMBER  OF  SRC  AD  GRGurS  IS  3 
OUTPUT  WILL  BE  FOR  3  3ROAD  GROUPS 
I  THF  BOUNDARY  F1NF  GROUPS  ARE  6  11  12 


vV 


NUCLIDE  IS  HYDROGFN 

ITS  NUMBE T  DENSITY  IS  6.6800E-02  PER  BARN-CM 
MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS,  FOLLOW. 


SCATTER 


ELASTIC 


ELASTIC 


TOTAL 


GROUP 

1 

2 

3 


TO 

INELASTIC 

N-2N 

(  PO> 

(PI  ) 

TRANSFER 

1 

.00000 

.00000 

1.82518 

1.46149 

1.82518 

2 

.00000 

.00000 

1.45849 

.72167 

1.45849 

3 

.00000 

.00000 

.00000 

.00000 

.00000 

2 

.00000 

.00000 

15.76453 

10.73447 

15.76453 

3 

.00000 

.00000 

1.19781 

.53015 

1.19781 

•3 

.00000 

.00000 

38.00000 

25.13320 

38.00000 

SIGTR 

SIGR 

SIGA  NUSIGF 

SOURCE 

1.10052 
5.70625 
13.  15730 


1.45849 

1.20634 

.29050 


i  OuOOO 
,00853 
,29050 


,00000 

>00000 

>00000 


.86840 

.13160 


NUCLIDE  IS  OXYGEN 

ITS  NUMBER  DENSITY  IS  3.3400E-02  PER  8ARN-CM 
MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS,  FOLLOW. 


SCATTER 
FROM  TO 


INELASTIC 

N-2N 

ELASTIC 

(PO) 

ELASTIC 
(PI  ) 

TOTAL 

“RANSFE 

.00244 

.00000 

2.44289 

.39193 

2.44532 

.00043 

.00000 

.11221 

.00256 

.11264 

.00000 

.00000 

.00000 

.00000 

.00000 

.COOJO 

.00000 

3.90007 

.15472 

3.90007 

.00000 

.00000 

.02377 

-.00713 

.02377 

.00000 

.00000 

4.2000C 

.17514 

4.20000 

GROUP 

1 

2 

3 


SIGTR 

2.17829 

3.77624 

4.02661 


SIGR 

,12746 

,02377 

.00175 


SIGA 

.01481 

.00000 

.00175 


NUSIGF 

.00000 

.00000 

.00000 


SOURCE 
.86840 
. 13160 


MAXWELL-BOLTZMAN  factor  =  1.128, 


AVERAGE  X  =  1 . 1428491 E  +  00 


.Figure  9c:  Water,  2  Pa3t  Groups,  Chain  2 
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JOB  OLD  BARNYARD. 

CHAIN  1 

PURE  URANIUM-235  1 

1 

FAST  GROUP  21 

FEB 

,  1967 

THE 

SOURCE  USED  IS  U235 

FISSION  I CRAN3ERG 

SPECTRUM) . 

SYSTEM  TEMPERATURE  IS 

293.0  KELVIN 

ATOM 

FRACTIONS  ARE 

URANIUM  235 

1.00000000 

RELATIVE  GROUP 

6R0UP  LETHARGY 

ENERGY.  EV 

SOURCE 

FLUX,  /UNIT  U 

1 

.50 

6.0653E+06 

2.3023E-02 

6.3595E+00 

2 

1.00 

3.6788E+06 

1.0824E-01 

3.3570E+01 

3 

1.50 

2.2313E+06 

2.1044E-01 

6.6725E+01 

4 

2,00 

1.3534E+06 

2.3139E-01 

8.9339E+01 

5 

2  <.50 

8.2085E+05 

1.8048E-01 

1.0000E+02 

6 

3 "  00 

4.9787E+05 

1.1483E-01 

9.9793E+01 

7 

5. CO 

6.7379E+04 

1.2439E-01 

6.6472E+01 

8 

8  .  CO 

3.3546E+03 

7.2094E-03 

4.7531E+00 

9 

12.00 

6.  14'/2E+01 

.oonoE+oo 

3.3288E-03 

10 

16.00 

1 . 1254E+00 

•OOOOE+OO 

8.0036E-08 

AGE  TO  INDIUM  RESONANCE  (1.46EV) 

IS 

1.9578E+01  CM2 

figure  lOat  Uranium-235 »  1  fast  Group, 

Chain  1 
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PLOT  OF  SPECTRUM,  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY,  FOLLOWS 


xxxxxxxx 


xxxxxxxxxxxx 


— i _ i _ i  — i  — _ j _ i_i _ t _ i_j _ i _ i _ r_i _ i _ i  — i  — i  — i  — : _ i _ j  — f 

l  l  t  l  l  l  i  i  i  i  *  l  i  >  i  i  i  i  i  i  i  t  i  l  i  i  <  i  i  i  t  !  i  I  l  i  l  t  l  l  i  l  l  l  i  :  1  l  )  l  l  t. 

.-i  o  in  op'®  r~-  «o  it  >a- 
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figure  10b %  Uranium-235,  1  Past  Group,  Flux  Plot 


JOB  OLD  BARNYARD  IHAIN  2Z 

THE  TOT  AL  NUMBER  OF  BROAD  GROUPS  IS  i 
OUT,JU'  ILL  BE  FOR  1  BROAD  GROUPS 
THE  8(  4uARY  FINE  GROUT  A  ARE  11  12 


NUCLIDE  iS  URANIUM  235 

ITS  HUMBER  DENSITY  IS  4.7900E-0 2  PER  BARN-CM 


MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS*  FOLLOW. 


SCATTER 
FROM  TO 

INELASTIC 

J 

ELASTIC 

N-2N  (PO) 

Elastic 

(Pi) 

TOTAL 

TRANSFER 

1  1 

1  2 

1.23713 

.00000 

• 

• 

00000  5.64214 

00000  .00000 

1.68886 

.00000 

6.87927 

.OCOOO 

GROUP 

l 

SIGTR 

6.63000  1. 

SIGR 

43959 

SIGA  NUSIGF 

1.43959  3.61426 

SOURCE 

1.00000 

figure  10c i  Uranium-235,  1  fast  group,  Chain  2 
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JOB  OLC  BARNYARD,  CHAIN  1 


l 


i 


t 


PURE  URANIUM-235  3  FAST  GROUPS  2\  FEB. 1967 


THE 

SOURCE  USED  IS  U235 

FISSION  <  CRANBERG 

SPECTRUM). 

SYSTEM  TEMPERATURE  IS 

293.0  KELVIN 

ATOM 

FRACTIONS  ARE 

URANIUM  235 

1.00000000 

RELATIVE  GROUP 

GROUP  LETHARGY 

ENERGY.  EV 

SOURCE 

FLUX.  /UNIT  U 

1 

.50 

6.0653E+06 

2.3023E-02 

6. 3595E+0C 

2 

1.0G 

3.6788E+06 

1.0824E-01 

3.3570E+0  *• 

3 

1.50 

2.2313E+06 

2.1044E-01 

6.6729E+0 1 

4 

2.00 

1 • 3534E+06 

2.3139E-01 

8.9339E+0 1 

5 

2.50 

8.2085E+05 

1.8048E-01 

1.0000E+02 

6 

3.00 

4.9787E+05 

1.1483E-01 

9.9793E+01 

7 

5.00 

6.7379E+04 

1.2439E-01 

6.6472E+0 1 

8 

8.00 

3.3546E+03 

7.2094E-03 

4. 7531E+00 

9 

12.00 

6.1442E+01 

•0000E+00 

3.8288E-03 

10 

16.00 

1.1254E+00 

•0000E+00 

8.0036E-08 

AGE  TO  INDIUM  RESONANCE  I1.46EV)  IS 

1.9578E+01  CM2 

ligure  Hat  Uranium-235 »  3  J?ast  Groups ,  Chain  1 
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PLOT  OF  SPECTRUM*  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY*  FOLLOWS 


JOB  OLD  BARNYARD  CHAIN  2Z 

THE  TOTAL  NUMBER  OF  BROAD  GROUPS  IS  3 
OUTPUT  WILL  BE  FOR  3  BROAD  GROUPS 
THE  BOUNDARY  FINE  GROUPS  ARE  3  6  11  12 


NUCLIDE  IS  URaNIUM  235 


ITS  NUMBER  DENSITY  IS  4.7900E-02  PER  BARN-CM 


MICRO 

SCOP  I  c 

CROSS  SECTIONS. 

IN  BARNS. 

FOLLOW. 

SCATTER 

ELASTIC 

ELASTIC 

FROM 

TO 

INELASTIC 

N-2N 

(PO) 

(PI) 

1 

1 

.05410 

.00000 

4.76586 

4.03780 

1 

2 

-1.17413 

.00000 

.02081 

-.61377 

1 

3 

.70041 

.00000 

.00000 

.00000 

1 

4 

.00000 

.00000 

.00000 

.00000 

2 

2 

.44240 

.00000 

4.07620 

1.57196 

2 

3 

.97542 

.00000 

.01892 

-.06592 

2 

4 

.00000 

.00000 

.00000 

.00000 

3 

3 

.80920 

.00000 

7.47117 

1.23983 

3 

4 

.00000 

.00000 

.00000 

.00000 

GROUP 

SIGTR 

SIGR 

SIGA 

NUSIGF 

SOURCE 

1 

4.57927 

3.18334 

1.28798 

3.91628 

.34170 

2 

5.26345 

2.25088 

1.25654 

3.26017 

.52670 

3 

8.71479 

1.67426 

1.67426 

3.85258 

.13160 

TOTAL 

TRANSFER 

4.81996 

1.19494 

.70041 

.00000 

4.51861 

.99434 

.00000 

8.28037 

.00000 


Figure  lies  Uranium-255*  3  Fast  Groups,  Chain  2 
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JOB  OLD  BARNYARD*  CHAIN  I 


NETF  CORE  (AN  MTR  TYPE  REACTOR)  1  FAST  GROUP  21  FEB.  1967 

THE  SOURCE  USED  IS  U235  FISSION  (CRANBERG  SPECTRUM). 

SYSTEM  TEMPERATURE  IS  293.0  KELVIN 
ATOM  FRACTIONS  ARE 


HYDROGEN 
OXYGEN 
ALUMINUM 
URANIUM  235 
URANIUM  238 


,48160826 

.24079822 

.27604967 

.00139125 

.00015263 


GROUP 

LETHARGY 

ENERGY,  EV 

SOURCE 

RELATIVE  GROUP 
FLUX.  /UNIT  U 

1 

.50 

610653E+06 

2.3023E-02 

1.3110E+01 

2 

1.00 

3. 6738E+06 

1.0824E-01 

5.0607E+01 

3 

1.50 

2.2313E+06 

2.1044E-01 

1.0000E+02 

4 

2.00 

1.3534E+06 

2.3139E-01 

9.9896E+01 

5 

2.50 

8.2085E+05 

1.8048E-01 

8.6102E+01 

6 

3.00 

4.9787E+05 

1.1483E-01 

7 . 2588E+0 1 

7 

5.00 

6.7379E+04 

1.2439E-01 

3.8833E+01 

8 

8.00 

3.3546E+03 

7.2094E-03 

2.0604E+01 

9 

12.00 

6.1442E+01 

•0000E+00 

1.8118E+01 

10 

16.00 

1.1254E+00 

•OOOOE+OO 

1.7 140E+0 1 

11 

17.00 

4.1399E-01 

•OOOOE+OO 

1 • 7765E+0 1 

12 

17.50 

2.51 10E-01 

.OOOOE+OO 

1.9334E+01 

13 

18.00 

1 • 5230E-01 

.OOOOE+OO 

2.5553E+01 

14 

IE. 50 

9.2374E-02 

.OOOOE+OO 

3.5858E+01 

15 

19.00 

5.6028E-02 

.OOOOE+OO 

3.6817E+01 

16 

19.50 

3.3983E-02 

•OOOOE+OO 

2.6404E+01 

17 

20.00 

2.0612E-02 

.OOOOE+OO 

1.4587E+01 

18 

20.50 

1.2502E-02 

.OOOOE+OO 

6.8143E+00 

19 

21.00 

7.5826E-03 

•OOOOE+OO 

2.8722E+00 

20 

21.50 

4. 5991E-03 

.OOOOE+OO 

1.1384E+00 

21 

22.00 

2 • 7895E-03 

.OOOOE+OO 

4.3532E-01 

22 

22.50 

1.6919E-03 

.OOOOE+OO 

1.6308E-01 

23 

23.00 

1.0262E-03 

.OOOOE+OO 

6.0405E-02 

AGE  TO  INDIUM  RESONANCE  (1.46EV)  IS 
AGE  TO  ARBITRARY  THERMAL  (1.12EV)  IS 
THERMAL  DIFFUSION  LENGTH  SQUARED  IS 
TOTAL  MIGRATION  AREA  IS 

MOST  PROBABLE  THERMAL  NEUTRON  VELOCITY  IS 
AVERAGE  THERMAL  NEUTRON  VELOCITY  IS 
ABSORPTION  PARAMETER  IS 
SCATTERING  PER  RESONANCE  ATOM  IS 
EFFECTIVE  RESONANCE  INTEGRAL  IS 

Figure  I2as  NETF,  i  Feat  Group,  Chain  1 


6.3752E+01  CM2 
6.3939E+01  CM2 
3, 1354E+00  CM2 
6.7075E+01  CM2 
2.8887E+03  M/SEC 
3.6? 31E+03  M/SEC 
1.0226E+00 
6.9315E+C4  BARNS 
2.5056E+02  8ARNS 


PLOT  OF  SPECTRUM.  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY.  FOLLOWS 


_  j.  _  _  •  t  , 


X 

X.  x 

X  X 


xxxxxxxxxxxxxxxxxxxx 


_  W  ‘  1  1  ‘  '  Jl,  ’  1  1  o  1  '  1  ’  1  *  ’  *  1  '  1  «  1  •  '  I  I  1  I  1  1  «  I  »  )  •  »  I  I  I 

.  ^  ^  °  O'  cd  in  *r\ 


figure  12bi  NJSTP,  1  Paat  Group,  Plux  Plot 


>-r-  -r 


'V  V>,»^ 


-’  *■«'-*  (^l--  w//  ■ 


X  X  X  X 


X  X 


X  X 


X  X 


xxxxxxxx 


xxxxxxxxxxxxxx*x  xxxxxx 


xxxxxxxxxxxx 


XXX 


5  I  I  *  I  l  I  ?  I  I  I  I  I  I  I  I  I  i  I  I  I  I  I  I  I  l  1  I  I  I  I  I  I  I  I  !  I  I  i  8  I  I  I  I  I  I  I  I  I  I  I  I  I  I 
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b 


100  1000 


J03  01.0  BARNYARD  CHAIN  2Z 

THE  TOTAL  NUMBER  OF  BROAD  GROUPS  IS  2 
OUTPUT  WILL  BE  FOR  2  BROAD  GROUPS 
THE  BOUNDARY  FINE  GROUPS  ARE  11  12 


NUCLIDE  IS  HYDROGEN 

ITS  NUMBER  OENSiTY  15  4.0713E-02  PER  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS,  FOLLOW. 


SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N 

1  PO) 

(PI) 

TRANSFER 

1 

1 

.00000 

.00000 

10.55896 

7 . 16066 

10.55896 

1 

2 

.000  '0 

.00000 

.66517 

.29354 

.66517 

2 

2 

.00000 

.00000 

38.00000 

25.13320 

38.00000 

GROUP 

SIGTR 

SIGR 

SIGA  NUSIGF 

SOURCE 

1 

3.77468 

66992 

.00475  .00000 

1.00000 

2 

13.06540 

19860 

.19860  *00000 

NUCLIDE 

IS  OXYGEN 

ITS 

NUMBER  DENSITY 

IS  2.0336E- 

02 

PER  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS 

,  FOLLOW. 

SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N 

(PO) 

(PI) 

TRANSFER 

1 

1 

.00116 

.00000 

3.34652 

.25410 

3.34767 

1 

2 

.00000 

.00000 

,01311 

-.00393 

.01311 

2 

2 

.00000 

.00000 

4.20000 

.17514 

4.20000 

GROUP 

SIGTR 

SIGR 

SIGA  NUSIGF 

SOURCE 

3.11666 

0191b 

.00605  .00000 

1.00000 

p 

4.02606 

00120 

.00120  .00000 

figure  12c i  NOT?,  1  fast  Group,  Chain  2 


'  -  T*»„  ..  , 


NUCLIDE  IS  ALI  MINUM 


ITS  NUMBER  DENSITY 

IS  2. 

3336E-02 

PER 

BARN-CM 

MICROSCOPIC  CROSS  SECTIONS* 

IN  BARNS 

,  FOLLOW. 

SCATTER 

ELASTIC 

elastic 

FROM  TO  INELASTIC 

N-2N 

(PO) 

(PI ) 

1  1  .05433 

.00000 

.00010 

.00093 

1  2  .00000 

.00000 

.OCOOO 

.00000 

2  2  .00000 

i 

.00000 

1 

.40000 

.03444 

GROUP  SIGTR 

SIGR 

SIGA 

NUSIGF 

SOURCE 

1  .06197 

00847 

.00847 

.00000 

1.00000 

2  1.50972 

14416 

.14416 

.00000 

NUCLIDE  IS  URANIUM  235 

ITS  NUMBER  DENSITY 

IS  1. 

1761E-04 

PER 

BARN-CM 

MICROSCOPIC  CROSS  SECTIONS. 

IN  BARNS 

,  FOLLOW. 

SCATTER 

elastic 

ELASTIC 

FROM  TO  INELASMC 

N-2N 

(PO) 

(PI  1 

1  1  .80758 

.00000 

7 

.81200 

1.1975V 

1  2  .00000 

.00000 

.00363 

-.00108 

2  2  .00000 

.00000 

14 

.80000 

.04144 

GROUP 

1 

2 


SIGTR 

24.27*99 

420.44«.80 


SIGR 

16.85293 

405.68824 


SIGA 

16.84930 

405.68824 


NUSIGF 

27.34165 

829.68092 


SOURCE 

1.00000 


figure  12d»  Hi£TPv  1  Past  Group,  Chain  2,  Continued 


TOTAL 

TRANSFER 

•05443 

.00000 

1.40000 


TOTAL 

transfer 

8.61957 

.00363 

14.80000 
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NUCLIDE  IS  URANIUM  238 


ITS  NUMBER  DENSITY  IS  1.2903E-05  PER  BARN-CM 
MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 


SCATTER 

FROM  TO 

INELASTIC 

N-2N 

ELASTIC 
(P  0) 

ELASTIC 
(PI ) 

1 

1 

.95131 

.00806 

12.89255 

1.33711 

1 

2 

.00000 

.00000 

.00208 

-.00063 

2 

2 

.  ocooo 

.00000 

13.80000 

.03864 

GROUP 

SIGTR 

SIGR 

StGA 

NUSlGF 

SOURCE 

I 

21.26361 

8.74012 

8.74610 

•  A  C 1 6 1 

1.00000 

2 

15.39440 

1.63304 

1.63304 

.00000 

MAXWELL-BOLTZMAN  FACTOR  =  1.128.  AVERAGE  X  =  1 .6717270F.+00 


figure  12e:  NiSTF,  1  fast  Croup,  Chain  2,  Continued 


total 

TPANSFER 

13.85997 

.00208 

13.80000 


r 

JOB 

OLD  BARNYARD.  CHAIN  1 

NETF 

CORE  (AN  MTR 

TY°E  REACTOR)  2  FAST  GROUPS  28 

FEB,  1967 

THE 

SOURCE  USED  IS  U235  FISSION  (CRANBERG 

SPECTRUM). 

SYSTEM  TEMPERATURE  IS  293.0  KELVIN 

ATOM 

FRACTIONS  ARE 

HYDROGEN 

.48160826 

OXYGEN 

.24079822 

ALUMINUM 

.27604967 

URANIUM  235 

.00139125 

URANIUM  238 

.00015263 

RELATIVE  GROUP 

GROUP 

LETHARGY 

ENERGY,  EV 

SOURCE 

FLUX,  /UNIT  U 

1 

.50 

6.0653E+06 

2.3023E-02 

1.3110E+01 

2 

1.00 

3.6788E+06 

1.0824E-01 

5.0607E+01 

3 

1.^0 

2.2313E+06 

2  •  1044E -01 

1.0000E+02 

4 

2.00 

1.3534E+06 

2.3139E-01 

9.9896E+01 

5 

2.50 

8  •  2085E+05 

1.8048E-01 

8.6102E+01 

6 

3.00 

4.9787E+05 

1.1483E-01 

7.2588E+01 

7 

5.00 

6.7379E+04 

1.2439E-01 

3.8833E+01 

3 

3.00 

3.3546E+03 

7.2094E-03 

2.0604E+01 

9 

12.00 

6.1442E+01 

•0000E+00 

1.6118E+01 

10 

16.00 

1.1254E+00 

•OOOOE+OO 

1.7140E+01 

11 

17.00 

4.1399E-01 

•OOOOE+OO 

1.7765E+01 

12 

17,50 

2.51 10E-01 

.0000E+00 

1.9334E+01 

13 

18.00 

1.5230E-01 

•OOOOE+OO 

2.5553E+01 

14 

13.50 

9.2374E-02 

•OOOOE+OO 

3.5858E+01 

15 

19.00 

5.6028E-0? 

•OOOOE+OO 

3.6817E+01 

16 

19.50 

3.3983E-02 

•OOOOE+OO 

2.6404E+0 1 

17 

20.00 

2.0612E-0 2 

.OOOOE+OO 

1.4587E+01 

18 

20.50 

1.2502E-02 

.OOOOE+OO 

6.8143E+00 

19 

21.00 

7 • 5826E-03 

•OOOOE+OO 

2.8722E+0C 

2C 

21.50 

4.5991E-03 

•OOOOE+OO 

1 . 1 384E+0C 

21 

22.00 

2 • 7895E-03 

.OOOOE+OO 

4.3532E-01 

22 

22.50 

1.6919E-03 

•OOOOE+OO 

1.6308E-01 

23 

23.00 

1.0262E-03 

.OOOOE+OO 

6.0405E-02  j 

AGE 

TO  INDIUM  RESONANCE  (1.46EV)  IS 

6.3752E+01  CM2 

AGE 

TO  ARBITRARY 

THERMAL  (1.12EV)  IS 

6.3939E+01  CM2 

THERMAL  DIFFUSION  LENGTH  SQUARED  IS 

3.1354E+00  CM2 

TOTAL  MIGRATION 

AREA  IS 

6.7075E+01  CM2 

MOST 

PROBABLE  THERMAL  NEUTRON  VELOCITY  IS  2.8887E+03  M/SEC 

AVERAGE  THERMAL 

NEUTRON  VELOCITY  IS 

3.6731E+03  M/SEC 

ABSORPTION  PARAMETER  IS 

1.0226E+00 

SCATTERING  PER  RESONANCE  ATOM  IS 

6.9315E+04  BARNS 

EFFECTIVE  RESONANCE  INTEGRAL  IS 

2.5056E+C2  BARNS 

Figure 

13a:  HfflW, 

2  Fast  Groups,  Chain 

1 
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PLOT  OF  SPECTRUM.  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY*  FOLLOWS 


o 

o 

o 


o 

o 


X 

X 


X 


X 

X 


X  X 

X 

X  X 

xxxxxxxxxxxx xxxxxxxxxx 

X 


X 


X 


X 


X 


X 


X 
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Figure  13b;  NjSIF,  2  Fast  Groups,  Flux  Plot 


i 

c\j 


JOB  OLD  RARNfARD  CHAIN  22 

THE  TOTAL  NUMBER  OF  BROAD  GROUPS  IS  3 
OUTPUT  WILL  BE  FOR  3  BROAD  GROUPS 
THE  BOUNDARY  r INF  GROUPS  ARE  6  11  12 


NUCLIDE  IS  HYDROGEN 

ITS  NUMBER  DENSITY  IS  4.C713E-02  PER  BARN -CM 


MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 


SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

fro;. 

TO 

INELASTIC 

N-2N  (PO) 

(PI ) 

TRANSFER 

1 

1 

.00000 

.00000  1.82604 

1.46253 

1.82604 

1 

2 

.00000 

.00000  1.46266 

.72397 

1.46266 

1 

3 

.00000 

.00000  .00000 

.00000 

.00000 

2 

2 

.00000 

.00000  15.74018 

10.70556 

15.74018 

2 

3 

.00000 

.00000  1.13920 

.50273 

1.13920 

3 

3 

.00000 

.00000  38.00000 

25.13320 

38.00000 

GROUP 

SIGTR  SI6R 

SIGA  NUSIGF 

SOURCE 

1 

1.10220  1.46266 

.00000  .00000 

.86840 

2 

5.67925  1.14734 

.00814  .00000 

.13160 

3  13.C6540  .19660  .19860  .00000 

NUCLIDE  IS  OXYGEN 

ITS  NUMBER  DENSITY  IS  2.0356E-02  PER  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 


SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELAST I c 

N-2N 

(PO) 

(PI) 

TRANSFER 

1 

1 

.00236 

00000 

2.44589 

.39171 

2.44825 

1 

2 

.00042 

00000 

.11264 

.00257 

.11306 

1 

3 

.PG000 

00000 

.00000 

.00000 

.00000 

2 

2 

.00000 

00000 

3.90808 

.15421 

3.90808 

2 

3 

.OOOOP 

00000 

.02245 

-.00673 

.02245 

3 

3 

.00000 

00000 

4.20000 

.17514 

4.20000 

GROUP 

SIGTR  SIGR 

SIGA  NUSIGF 

SOURCE 

1 

2.18157  .12760 

• 

0i'+54  .00000 

.86840 

n 

3.78305  .02245 

• 

00000  .00000 

.13160 

3 

4.02606  .00120 

• 

00120  .00000 

figure  13c;  NiSTf,  2  fast  Groups,  Chain  2 
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NUCLIDE  15  ALUMINUM 


ITS  NUMBER  DENSITY 

IS  2.3336E-02 

PER  6ARN-CM 

MICROSCOPIC  CROSS  SECTIONS,  IN  8ARNS 

,  FOLLOW. 

SCATTER 

ELASTIC 

ELASTIC 

total 

FROM 

TO  INELASTIC 

ft-  2N 

<PO> 

(PI  ) 

TRANSFER 

1 

1  .13057 

.00000 

.00024 

.00224 

.13081 

1 

2  .00000 

.00000 

.00000 

.00000 

.00000 

1 

3  .00000 

.00000 

.00000 

.00000 

.00000 

2 

2  .C0000 

.00000 

.00000 

.00000 

.00000 

2 

3  .00000 

.00000 

.00000 

.00000 

.00000 

3 

3  .00000 

.00000 

1.40000 

.03444 

1.40000 

GROUP 

SIGTR 

SIGR 

SIGA  NUSIGF 

SOURCE 

1 

.13;:89 

00432  .00432  .00000 

.86840 

2 

.01143 

01143  .01143  .00000 

,13160 

3 

1.50972 

14416  .14416  .00000 

NUCLIDE  IS  URANIUM  235 

ITS  NUMBER  DENSITY 

IS  1.1761E-04 

PER  BARN-CM 

MICROS 

COPIC  CROSS  SECTIONS,  IN  BARNS 

,  FOLLOW. 

SCATTER 

ELASTIC 

ELASTIC 

total 

FROM 

TO  INELASTIC 

N-2N 

(PO) 

(PI  > 

TRANSFER 

1 

1  .77234 

.00000 

4.31206 

2.30494 

5.08440 

1 

2  .85655 

.00000 

.00884 

-.03074 

.86540 

1 

3  .00000 

.00000 

.00000 

.00000 

.00000 

2 

2  .22226 

.00000 

10.29996 

.43032 

10.52222 

2 

3  .00000 

.00000 

.00621 

-.00184 

.00621 

3 

3  .00000 

.00000 

14.80000 

.04144 

14.80000 

GROUP 

SIGTR 

SIGR 

SIGA  NUSIGF 

SOURCE 

1 

4.95087  2. 

14068  1.27528  3.55094 

.86840 

2 

38.04826  27. 

95452  27.94830  44.29636 

.13160 

3 

420.44680  405. 

68824  405.68824  829.68092 

1 

figure  13d *  KBTF,  2  Fast  Groups,  Chain  2  Continued 


NUCLIDE  IS  URANIUM  2  38 

ITS  NUMBER  DENSITY  IS  1.2903E-05  PER  BARN-CM 

MICROSCOPIC  Cl  OSS  SECTIONS,  IN  RARNS*  FOLLOW. 


SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N  IPO) 

(  PI ) 

TRANSFER 

1 

X 

1.03487 

.01530  4.95898 

2.44602 

6.02445 

1 

2 

1.0608S 

.00406  .01103 

.01275 

1.08000 

1 

3 

.00000 

.00000  .00000 

.00000 

.00000 

2 

2 

.  13573 

.00000  18.53863 

.53775 

18.67436 

2 

3 

.00000 

.00000  .00356 

-.00107 

.00356 

3 

3 

.00000 

.00000  13.80000 

.03864 

13.80000 

GROUP 

SIGTR  SIGR 

SI6A  NUSIGF 

SOURCE 

1 

5.03425  1.46857 

.40793  .96514 

.86840 

2 

32.82964  14.69195 

14.68839  .00000 

.13160 

3 

15.39440  1.63304 

1.63304  .00000 

MAXWFLL-BOLT2MAN  FACTOR  =  1.128,  AVERAGE  X  =  1.6717270E+00 


figure  13e*  HETF,  2  Fast  Groups,  Chain  2  Continued 


JOB  OLD  BARNYARD.  CHAIN  1 


AFIT  SUBCR  I  T  I  CAL  REACTOR  CORE  1  FAST  GROUP 


THE  SOURCE  USED  IS  U235  FISSION  (CRANBERG  SPECTRUM). 
SYSTEM  TEMPERATURE  IS  293.0  KELVIN 
ATOM  FRACTIONS  ARE 


HYDROGEN 
OXYGEN 
ALUMINUM 
URANIUM  235 
URANIUM  238 


.46238076 

-.23119038 

.08513929 

.00161649 

.22934297 


RELATIVE  GROUP 


GROUP 

LETHARGY 

ENERGY,  EV 

SOURCE 

FLUX,  /UNIT 

1 

.50 

6.0653E+06 

2.3023E-02 

9.4603E+00 

2 

1.00 

3.6788E+06 

1.0824E-01 

4.2270E+01 

3 

1.50 

2.2313E+06 

2. 1044E-01 

8.4042E+01 

4 

2.00 

1.3534E+06 

2.3139C-01 

9.3733E+01 

5 

2.50 

8.2085E+05 

1 • 8048E-0 1 

1.0000E+02 

6 

3.00 

4.9787E+05 

1.1483E-01 

9.9286E+01 

7 

5.00 

6.7379E+04 

1.2439E-01 

6.0211 £+01 

8 

8.00 

3.3546E+03 

7.2094E-03 

3.1922E+01 

9 

12.00 

6.1442E+01 

• OOOOE+OO 

2.6493E+01 

10 

16.00 

1 . 1254E+00 

•OOOOE+OO 

1 • 7623E+0 1 

11 

17.00 

4, 1399E-01 

• 000CE+00 

1.8736E+01 

12 

17.50 

2.5110E-01 

•OOOOE+OO 

2.1139E+01 

13 

18.00 

1.5230E-01 

•0000E+00 

3.3059E+01 

14 

18.50 

9.2374E-02 

.OOOOE+OO 

5.8141E+01 

15 

IS. 00 

5.6028E-02 

.OOOOE+OO 

6.8442E+01 

16 

19.50 

3.3983E-02 

.OOOOE+OO 

5 . 3053E+0 1 

17 

20.00 

2.0612E-02 

•OOOOE+OO 

3.0782E+01 

18 

20.50 

1.2502E-02 

.OOOOE+OO 

1.4872E+01 

19 

21.00 

7.5826E-03 

.OOOOE+OO 

6.4232E+00 

20 

21.50 

4.5991E-03 

.OOOOE+OO 

2 • 5928E+00 

21 

22.00 

2.7895E-03 

.OOOOE+OO 

1.0053E+00 

22 

22.50 

1.6919E-0? 

.OOOOE+OO 

3.8067E-01 

23 

23.00 

1.G262E-03 

.OOOOE+OO 

1.4218E-01 

AGE  TO  INDIUM  RESONANCE  (1.46EV)  IS 
AGE  TO  ARBITRARY  THERMAL  (1.12EV)  IS 
THERMAL  DIFFUSION  LENGTH  SQUARED  IS 
TOTAL  MIGRATION  AREA  IS 
MOST  PROBABLE  THERMAL  NEUTRON  VELOCITY 
AVERAGE  THERMAL  NEUTRON  VELOCITY  IS 
ABSORPT’CN  PARAMETER  IS 
SCATTERING  PER  RESONANCE  ATOM  IS 
EFFECTIVE  RESONANCE  INTEGRAL  IS 


2.9947E+01  CM2 
3.0004E+01  CM2 
1.8123E+00  CM2 
3.1 8 17E+0 1  CM2 
2 • 5493E+03  M/SEC 
3.2063E+03  M/SEC 
5  9 8382E-0 1 
7 • 6640E+0 1  BARNS 
2.0765E+01  BARNS 


figure  14a i  Afll  Subcrit.,  1  fast  Group,  Chain  1 
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PLOT  OF  SPECTRUM*  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY.  FOLLOWS 


o 

c 

c 


4 


o 

o 


X 

X  x 


XXXXXXXXX.  XXXXXXKXXXX 


o  ••••••••••• 

*-«  X 


X 


X 


X 


X 


X 

X 

X 


1  1  X 

1  r-t  Cf\ 

1  1  1  1 
CNJ 

1  1  1  1  1  1 
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•* 
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■# 
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Figure  14b t 

JLFIT 

Subcrit., 

1  Fast  Group,  Flux  Plot 

JOB  OLD  BARNYARD  CHAIN  2Z 


THE  TOTAL  NUMBER  OF  BROAD  GROUPS  IS  2 
OUTPUT  WILL  BE  FOR  2  BROAD  GROUPS 
THE  BOUNDARY  FINE  GROUPS  ARE  11  12 


NUCLIDE  IS  HYDROGEN 

ITS  NUMBER  DENSITY  IS  3.5040E-02  PER  BARN-CM 


MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 

SCATTER  ELASTIC 


SCATTER 
FROM  TO 


ELASTIC 


TOTAL 


GROUP 

1 

2 


INELASTIC 

N-2N 

(PO) 

!  PI  > 

TRANSFER 

.00000 

.00000 

11.44637 

7.72347 

11.44837 

.00000 

.00000 

.55438 

.24351 

.55438 

.00000 

.00000 

38.00000 

25.13320 

38.00000 

R 

SIGR 

SIGA  NUSIGF 

SOURCE 

4.03*'72 

13.09431 


.55834 

.22751 


.00396 

.22751 


.00000 

.00000 


1.00000 


NUCLIDE  IS  OXYGEN 

ITS  NUMBER  DENSITY  IS  1.7520E-02  PER  BARN-CM 


MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 


SCATTER 
FROM  TO 


INELASTIC 

N-2N 

ELASTIC 

(PO) 

ELASTIC 
(PI ) 

TOTAL 

TRANSFER 

. C0068 

.00000 

3.52317 

.24283 

3.52385 

.00000 

.00000 

.01084 

-.00325 

.01084 

.00000 

.00000 

4.20000 

.17514 

4.20000 

GROUP 

1 

2 


SIGTR 

3.29896 

4.02623 


SIGR 

.01469 

.00137 


SIGA 

.00385 

.00137 


NUSIGF 

.00000 

.00000 


SOURCE 

1.00000 


figure  14c:  Afll  Subcrit. ,  1  fast  Broup,  Chain  2 
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nuclide  is  aluminum 


ITS  NUMBER  DENSITY  • S  6.4520E-03  PER  6ARN-CM 
MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS,  FOLLOW. 


SCATTER 

elastic 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N 

(PO) 

!  PI  I 

TRANSFER 

1 

1 

.04296 

.00000 

.00006 

.00055 

.04303 

1 

2 

.COOOC 

.OOOOC 

.00000 

.00000 

.00000 

2 

2 

.00000 

.00000 

1.40000 

.03444 

1.40000 

GROU° 

SIGTR 

SIGR 

S!GA 

NUSIGF 

SOURCE 

1 

.05049 

.00800 

.00300 

.coooo 

1.00000 

2 

1.53071 

.16515 

.16515 

.00000 

NUCLIDE  IS  URANIUM  235 

ITS  NUMBER  DENSITY  IS  1.2250E-04  PER  EARN-CM 
MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS.  FOLLOW. 


SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC  N-2N 

(PO) 

vPl  ) 

TRANSFER 

1 

1 

.70799  .00000 

6.06126 

1.04389 

8.76924 

1 

2 

.00000  .OOCOO 

.00300 

-.00089 

.00300 

2 

2 

.00000  ,00000 

14.80000 

.04144 

14.30000 

GROUP 

SIGTR  SIGR 

SiGA  NUSIGF 

SOURCE 

1  2 3.28539  15.55915  15.55615  25,34965  1.00000 

2  479.51266  464 .75410  464.7S41Q  950,47764 


figure  i4d:  AFII  Subcnt.,  1  fast  Group,  Chain  2,  Continued 
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NUCLIDE  IS  URANIUM  238 


ITS  NUMBER  DENSITY  IS  1.7380E-02  PER  BARN-CM 


MICROSCOPIC 

CROSS  SECTIONS, 

IN  BARNS* 

FOLLOW. 

SCATTER 

FROM  TO 

INELASTIC 

N-2N 

ELASTIC 
(  PO) 

ELASTIC 
t  PI  ) 

TOTAL 

TRANSFER 

1  1 

1  2 

2  2 

.77921 

.00000 

.00000 

.00474 

.00000 

.00000 

13.38083 

.00172 

13.80000 

1.19954 

-.0005? 

.03864 

14.16952 

.00172 

13.80000 

GROUP 

SIGTR 

SIGR 

1 

13.77623 

.80574 

2 

15.63216 

1.87080 

SIGA 

.80875 

1.87080 


NUSIGF 

.28073 

.00000 


SOURCE 

1.00000 


MAXWELL-BOLTZMAN  FACTOR  =  1.128, 


figure  Hex  API!  Subcrit.,  1 


AVERAGE  X  =  1.4592663E+00 


fast  Group,  Oiiain  2,  Continued 
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JOB  OLD  BARNYARD  .  CHAIN  1 


AFIT  SU8CR  I  T  I CAi.  REACTOR  CORE  2  FAST  GROUPS 


THE  SOURCE  USED  IS  U235  FISSION  (CRANtJtRG  SPtCTRuM). 
SYSTEM  TEMPERATURE  IS  293.0  KELVIN 
ATOM  FRACTIONS  ARE 


HYDROGEN 
OXYGEN 
ALUMINUM 
URANIUM  23? 
URANIUM  238 


.46238076 

,23119038 

.08S13929 

.00161649 

.22934297 


RELATIVE  GROUP 


GROUP 

LETHARGY 

ENERGY.  EV 

SOURCE 

FLUX.  /UNIT 

1 

.50 

6.0653E+06 

2.3023E-02 

9.4603E+00 

2 

1.00 

3.6788E+06 

1.0824E-01 

4-.2270E+01 

3 

1.50 

2.2313E+06 

2.1044t-01 

0.4042E+C1 

4 

2.00 

1 ♦ 3534E+06 

2.3139E-01 

9.3733E+01 

5 

2.50 

8.2085E+05 

1.8048E-01 

1.0000E+02 

6 

3.00 

4.9787E+05 

1.1483E-01 

9.9286E+01 

7 

5.00 

6.7379E+0A 

1.243yE-0i 

6.0211E+01 

8 

8.00 

3.3546E+03 

7.2094E-03 

3.1922E+01 

9 

12.00 

6.1442E+01 

• OOOOE+OO 

2.6493E+01 

1C 

16.00 

1 ■ 1254E+00 

, OOOOE+OO 

1.7622S+01 

11 

17.00 

4.1399E-01 

•OOOGE+OO 

1 .8736E+01 

12 

17.50 

2.5110E-01 

.OOOOE+OO 

2.1139E+01 

13 

18.00 

1.5230E-0? 

•OOOOE+OO 

3.3059E+01 

14 

18.50 

9.2374E-02 

.OOOOE+OO 

5.8141F-r0l 

15 

19.00 

5.6O28E-02 

•OOOOE+OO 

6.8442E+01 

16 

19.50 

3.3983E-02 

.OOOOE+OO 

5 • 3053E+0 1 

17 

20.00 

2 . 0612E-02 

•OOOOF+OO 

3.0/82E+01 

18 

20.50 

1.2502E-02 

.OOCGE+OG 

1.4872E+01 

19 

21.00 

7.582&E-03 

.OOOOE+OO 

6.4232E+00 

20 

21.50 

4.599.-E-03 

•OOOCE+OO 

2 . 5928E+00 

21 

22.00 

2 ^  7895E-03 

.OOOOt+OO 

1.0C53E+00 

22 

22.50 

1.6919E-03 

.OOOOE+OO 

3»8067E-0l 

23 

23.00 

1  < 0262E-03 

.OOOOE+OO 

1.4216E-01 

AGE  TO  INDIUM  RESONANCE  (1.46EV)  IS  2.9947E+01  CM2 

AGE  TO  ARBITRARY  THERMAL  (1.12EVJ  IS  3.00Q4E+01  CM2 

THERMAL  DIFFUSION  LENGTH  SQUARED  IS  1.8123E+00  CM2 

TOTAL  MIGRATION  AREA  IS  3.1817E+01  CM2 

MOST  PR08ABLE  THERMAL  NEUTRON  VELOCITY  IS  2.5493E+03  M/SEC 

AVERAGE  THERMAL  NEUTRON  VELOCITY  IS  3.2063E+03  M/SEC 

ABSORPTION  PARAMETER  IS  5.8382E-01 

SCATTERING  PER  RESONANCE  ATOM  IS  7.664QE+01  BARNS 

EFFECTIVE  RESONANCE  INTEGRAL  iS  2.0765E+G1  BARNS 


■figure  15at  AFIT  Subcrit.,  2  Fast  Groups,  Chain  X 
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LOT  OF  SPECTRUM,  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY,  FOLLOWS 


o 

o 

o 


o 

o 


X 

X 


X 


X 


X 


X 


X 


X 


X 

X 

X  X 

xxxxxxxxxxxxxxxxxxxx 


X  X  X  > 


X 


X 


X 


X 


X 


X 


X 


X 


O  1  X 

»h  ff» 

till 

eg 

i  i  i  i  i  i 

H 

1 

O 

i  i  i  i  i  i  i  t 
O'  00 

1  1  1  1  1  1 

ti  •  r\i 

rsj 

f\J 

rv 

r-«  ,-H 

fH 

figure  15 b» 

AFIT 

Suberit., 

2 

Fast  Groups, 

Flux  Plot 

A 


:  i  f  i  i  i  t  i  s  i  (  I 

^  cn  CNJ 


JOB  OLD  BARNYARD  CHAIN  2Z 


THE  TOTAL  NUMBER  OF  BROAD  GROUPS  IS  3 
OUTPUT  WILL  BE  FOR  3  BROAD  GROUPS 
THE  BOUNDARY  FINE  GROUPS  ARE  6  11  12 


NUCLIDE  IS  HYDROGEN 


ITS  NUMBER  DENSITY  IS  3.5040E-02  PER  BARN-CM 


MICRO 

SCOP  I C 

CROSS  SECTIONS. 

IN  BARNS. 

FOLLOW. 

SCATT 

ER 

ELASTIC 

ELASTIC 

FROM 

TO 

INELASTIC 

N-2N 

(PO) 

(PI) 

1 

1 

.00000 

.00000 

1.79736 

1.45938 

1 

2 

.00000 

.00000 

1.73146 

.88732 

1 

3 

.00000 

.00000 

.00000 

.00000 

2 

2 

•OOOOC 

.00000 

15.60306 

10.54419 

2 

3 

.00000 

.00000 

.84521 

.37127 

3 

3 

.00000 

.00000 

38.00000 

25.13320 

GROUP 

SIGTR 

SIGR 

SIGA 

NUSIGF 

SOURCE 

1 

1.18213 

1.73147 

.00000 

.00000 

.86840 

2 

5.53885 

.35125 

•00604 

.00000 

.13160 

3 

13.09431 

.22751 

.22751 

.00000 

N-  .LIDE  IS  OXYGEN 

ITS  NUMREP  DENSITY  IS  1.7520E-02  PER  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS*  IN  BARNS*  FOLLOW. 


SCATT 

FROM 

ER 

TO 

INELASTIC 

N-2N 

ELASTIC 

(PO) 

ELASTIC 

(PI) 

1 

1 

.00168 

.00000 

2.56149 

.41282 

1 

2 

.00030 

.00000 

.15174 

.00346 

1 

3 

.00000 

.00000 

.00000 

.00000 

2 

2 

.00000 

.00000 

3.94807 

.15183 

2 

3 

.00000 

.00000 

.01653 

-.00495 

3 

3 

.00000 

.00000 

4.20000 

.17514 

GROUP 

SIGTR 

SIGR 

SIGA 

NUSIGF 

SOURCE 

1 

2.31011 

. 16321 

.01118 

.00000 

. 86840 

2 

3.81772 

.01653 

.00000 

.00000 

.13160 

3 

4.C2623 

.00137 

.00137 

.00000 

figure  15c*  API! 

Subcrit., 

2  Past  Groups* 

Chain  2 

TOTAL 

TRANSFER 

1.79736 

1.73146 

.00000 

15.60306 

.84521 

38.00000 


TOTAL 

TRANSFER 

2.56317 

.15204 

.00000 

3.94807 

.01653 

4.20000 


NUCLIDE  IS  ALUMINUM 


ITS  NUMBER  DENSITY  IS  6.4520E-03  PCR  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 


SCATTER 

FTOM  TO 

inelastic 

N-2N 

ELASTIC 
(  P0> 

ELASTIC 
(PI  ) 

1 

1 

.12489 

.00000 

.00017 

.00159 

1 

2 

.00000 

.00000 

.00000 

.00000 

1 

3 

.00000 

.00000 

.00000 

.00000 

2 

2 

.00000 

.00000 

.00000 

.00000 

2 

3 

.oocco 

.00000 

.00000 

.00000 

3 

3 

.00000 

.00000 

1.40000 

.03444 

GROUP 

SIGT'R 

SIGR 

SIGA 

NUSIGF 

SOURCE 

1 

.12679 

.00331 

.00331 

.00000 

.36840 

2 

.01046 

.01046 

.01046 

.00000 

.13160 

3 

1.53071 

.16515 

.16515 

.00000 

NUCLIDE  IS  URANIUM  235 

ITS  NUMBER  DENSITY  IS  1.2250E-04  PER  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS*  FOLLOW. 


SCATT 

FROM 

ER 

TO 

INELASTIC 

N-2N 

ELASTIC 

(P0> 

ELASTIC 

(PI) 

1 

1 

•69454 

.00000 

4.28541 

2.16289 

1 

2 

.88675 

.00000 

.01191 

-.04141 

1 

3 

.00000 

.00000 

.00000 

.00000 

2 

2 

.24984 

.00000 

10.03586 

.47857 

2 

3 

.00000 

.ooooc 

•  00457 

-.00136 

3 

O 

.OOCOC 

.00000 

14.80000 

*04 144 

GROUP 

SIGTR 

SIGR 

SIGA 

NUSIGF 

SOURCE 

1 

5.02510 

2.  i6663 

1.26798 

3.47598 

•  86840 

2 

32.86496 

23.05647 

23.05189 

36.82483 

•i3160 

3 

479.51266 

464.75410 

464.75410 

950.47764 

ifigure  15dt  AFIT  Subcrit.,  2  Fast  Groups,  Ciaain  2,  Continued 


total 

transfer 

.  1250/' 
.00000 
.00000 
.oocoo 
.00000 
1.4Q000 


TOTAL 

transfer 

4.97995 

.89866 

.00000 

10.26571 

.00457 

14.80000 


NUCLIDE  IS  URANIUM  238 

ITS  NUMBER  DENSITY  IS  1.7380E-02  PER  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS,  FOLLOW. 


SCATTER 
FROM  TO 


INELASTIC 

N-2N 

ELASTIC 

IPO) 

ELASTIC 
(PI  ) 

TOTAL 

TRANSFE 

.92311 

.01087 

5.04366 

2.33414 

5.98852 

1.05061 

.00289 

.01486 

.01718 

1.07124 

.00000 

.00000 

.00000 

.00000 

.00000 

.15256 

.00000 

17.74682 

.59531 

17.89938 

.00000 

.00000 

.00262 

-.00079 

.00262 

.00000 

.00000 

13.80000 

.03864 

13.80000 

GROUP 

1 

2 
3 


SIGTR 

5.06187 

18.34789 

15.63216 


SIGR 

1.42467 

1.04302 

1.87080 


SIGA 

.36719 

1.04040 

1.87080 


NUSIGF 

.81584 

.00000 

.00000 


SOURCE 

.86840 

.13160 


MAXWTJLL-BOLTZMAN  FACTOR  =  1.128,  AVERAGE  X  =  1 .4592663E+00 


figure  15e:  AFIT  Subcrit.,  2  fast  Groups,  Chain  2,  Continued 


JOB  OLD  BARNYARD.  CHAIN  1 


ROMASKA  CORF  (A  DIRECT  ENERGY  CONVERTER  REACTOR)  1  FAST  GROUP 

THE  SOURCE  US  ED  IS  U235  FISSION  (CRANBERG  SPECTRUM). 

SYSTEM  TEMPERATURE  IS  1273.0  KELVIN 
ATOM  FRACTIONS  ARE 


CARBON 

.86548694 

uranium 

235 

.12213834 

URANIUM 

238 

.01361218 

GROUP 

LETHARGY 

ENERGY,  EV 

SOURCE 

RELATIVE  GROUP 
FLUX,  /UNIT  U 

1 

.50 

6.0653E+06 

2.3023E-02 

4.4849E+00 

2 

1.00 

3.6788E+06 

1.0824E-01 

2 . 3895E+0 1 

3 

1.50 

2.2313E+06 

2. 104^E-01 

5.8781E+01 

4 

2.00 

1 . 3534E  +-06 

2.3139E-01 

8.9394E+01 

5 

2.50 

8. 2035E+05 

1.8043F-01 

9.1312E+01 

6 

3.00 

4.9787E+05 

1.1483E-01 

9 • 0867E+0 1 

7 

5.00 

6. 7379E+04 

1.2439E-01 

1.0000E+02 

8 

8.00 

3.3546E+03 

7.2094E-03 

3.4544E+01 

9 

12.00 

6.1442E+01 

. OOOOE+OO 

6.9175E-02 

10 

16.00 

1  •  1254E+00 

•0000E+00 

9.7373E-07 

11 

17.00 

4.1399E-01 

. OOOOE+OO 

1.0213E-11 

12 

17.50 

2.51 10E-01 

•0000E+00 

4.7077E-14 

13 

18.00 

1.5230E-01 

. 0000E+00 

3.0930E-16 

14 

18.50 

9.2374E-02 

• OOOOE+OO 

2.9839E-18 

15 

19.00 

5.6028E-02 

.OOOOE+OO 

4.2305E-20 

16 

19.50 

3.3983E-02 

•00C0E+00 

8.6734E-22 

17 

20.00 

2.0612E-02 

.OOOOE+OO 

2.5074E-23 

18 

20.50 

1 . 2502E-02 

.OOOOE+OO 

9.9264E -25 

19 

21.00 

7.5826E-  03 

.0000E+00 

5.2209E-26 

20 

21.50 

4.5991E-03 

.OOOOE+OO 

3 . 5424E-27 

21 

22.00 

2.7895E-03 

.OOOOE+OO 

3.0158E-28 

22 

22.50 

3.6919E-03 

.OOOOE+OO 

3.1401E-29 

23 

23.00 

1.0262E-03 

.OOOOE+OO 

3.9066E-30 

AGE  TO  UDIUM  RESONANCE  U.46EV)  IS 
AGE  TO  ARBITRARY  THERMAL  (1.12EV)  IS 
THERMAL  DIFFUSION  LENGTH  SQUARED  IS 
TOTAL  MIGRATION  AREA  IS 

MOST  PROP ABLF  THERMAL  NEUTRON  VELOCITY  IS 
AVERAGE  THERMAL  NEUTRON  VELOCITY  IS 
ABSORPTION  PARAMETER  IS 
SCATTERING  PER  RESONANCE  ATOM  IS 
EFFECTIVE  RESONANCE  INTEGRAL  IS 

.Figure  l6a»  riontaska,  1  Fast  Group,  Chain  1 


9.6172E+01  CM2 
9 . 5923E+0 1  CM2 
1 . 7292E-0 1  CM2 
9 • 6096E+0 1  CM2 
. OOOOE+OO  M/SEC 
2.1557E+04  M/SEC 
2.7659E+0 2 
4.3363E+02  BARNS 
4.6972E+01  BARNS 
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JOB  OLD  9ARN YARD  CHAIN  21 


THE  TOTAL  NUME 2R  OF  BROAD  GROUPS  IS  1 
OUTPUT  WILL  BE  FOR  1  BROAD  GROUPS 
THE  30UNDAR Y  FINt  GROUPS  ARt  11  12 


NUCLIDE  IS  CARBON 

ITS  NUMBER  DENSITY  IS  6.99AOE-02  PER  BARN-CM 


MICROSCOPIC  CROSS  SECTIONS*  IN  BARNS*  FOLLOW. 


SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N 

(PO) 

(PI ) 

TRANSFER 

1 

1 

.00211 

.00000 

3.51316 

.28160 

3.51527 

1 

2 

.00000 

.00000 

.00000 

.00000 

.00000 

GROUP 

SIGTR  SIGR 

SI  GA  NUSIGF 

SOURCE 

1 

3.23367  .00000 

.00000  .00000 

1.00000 

NUCLIDE 

IS  URANIUM  235 

ITS 

NUMBER  DENSITY  IS  9. 

8700E-03 

PER  BARN-CM 

MICROSCOPIC  CROSS  SECTIONS* 

IN  E’.ARNS 

,  FOLLOW. 

SCATTER 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N 

(PO) 

(PI ) 

TRANSFER 

1 

1 

.94420 

.00000 

6.76122 

1.42560 

7.70542 

1 

2 

.00000 

.00000 

.00000 

.00000 

.00000 

GROUP 

SIGTR  SIGR 

SIGA  NUSIGF 

SOURCE 

1 

7.99061  1.71079 

1.71079  4.01466 

1.00000 

Figure  16c s  Romas ka,  1  Fast  Group*  Chain  2 


NUCLIDE  IS  URANIUM  238 

ITS  NUMBER  DENSITY  IS  1.1000E-03  PER  BARN -CM 


MICROSCOPIC  CROSS  SECTIONS,  IN  BARNS, 

FOLLOW. 

SCATTER 
FROM  TO 

INELASTIC 

N-2N 

ELASTIC 
t  PO! 

ELASTIC 
(PI  ) 

total 

transfer 

1  1 

1  2 

,94824 

.OOOOO 

.00285 

.00000 

8.09556 

.00000 

1.67050 

.00000 

9.04950 

.00000 

GROUP 

1 

SIGT7 

7.57219 

SIGR  SIGA  NUSIGF 

19318  .19603  .26603 

SOURCE 

1.00000 

JOB  OLD  BARNYARD.  CHAIN  1 


ROMASKA  CORE  (A  DIREC  ENERGY  CONVERTER  REACTOR)  3  FAST  GROUP 


THE  SOURCE  USED  IS  U235  FISSION  (CRANBERG  SPECTRUM). 
SYSTEM  TEMPERATURE  IS  1273.0  KELVIN 
ATOM  FRACTIONS  ARE 


CARBON 

.86548694 

URANIUM  235 

.12213834 

URANIUM  238 

.01361218 

GROUP 

LETHARGY 

ENERGY,  EV 

SOURCE 

RELATIVE  GROUP 
FLUX,  /UNIT  U 

1 

.50 

6. 0653E+06 

2.3023E -02 

4.4849E+00 

2 

1.00 

3. 6788E+06 

1.0824E-01 

2.3895E+01 

3 

1.50 

2.2313E+06 

2.1044E-01 

5.8781E+01 

4 

2.00 

1.3534E+06 

2.3139E-01 

8.9394E+01 

5 

2.50 

8.2085E+05 

1.8048E-01 

9.1312E+01 

6 

3.00 

4.9787E+05 

1.1483E-01 

9.0867E+01 

7 

5.00 

6.7379E+04 

1.2439E-01 

1.0000E+02 

8 

8.00 

3.3546E+03 

7.2094E-03 

3.4544E+C-1 

9 

12.00 

6.1442E+01 

.OOOCE+OO 

6.9175E-02 

10 

1  S.00 

1. 1254E+00 

.OOOOE+OO 

9.7373E-07 

11 

i  r.-oo 

4. 1399E-01 

•OCOOE+OO 

1.0213E-U 

12 

17.50 

2.5110E-01 

•OOOOE+OO 

4.7077E-14 

13 

18.00 

1.5230E-01 

•OOOOE+OO 

3.0930E-16 

14 

13.50 

9.2374E-02 

•OOOOE+OO 

2.9839E-18 

15 

19-00 

5.6028 E-02 

•OOOOE+OO 

4.2305E-20 

16 

1  9.50 

3.3983E-02 

.OOOOE+OO 

8.6734E-22 

17 

20.00 

2.0612E-02 

•OOOOE+OO 

2.5074E-23 

18 

20.50 

1.2502E-02 

.OOOOE+OO 

9.9264E-25 

19 

21.00 

7.5826E-03 

.OOOOE+OO 

5.2209E-26 

20 

71.50 

4.5991E-03 

•OOOOE+OO 

3.5424E-27 

21 

22.00 

2.7895E-03 

•OOOOE+OO 

3.0158E-28 

22 

22.50 

1..6919E-03 

.OOOOE+OO 

3.1401E-29 

23 

23.00 

1.0262E-03 

.OOOOE+OO 

3.9066E-30 

AGE  TO  INDIUM  RESONANCE  C1.46EV)  13 
AGE  TO  ARBITRARY  THERMAL  (1.12EV)  IS 
THERMAL  DIFFUSION  LENGTH  SQUARED  IS 
TOTAL  MIGRATION  APIA  IS 

MOST  PROBABLE  THERMAL  NEUTRON  VELOCITY  IS 
AVERAGE  THERMAL  NEUTRON  VELOCITY  IS 
ABSORPTION  PARAMETER  IS 
SCATTERING  PER  RESONANCE  ATOM  IS 
EFFECTIVE  RESONANCE  INTEGRAL  IS 

Figure  17a;  fiomaska,  3  Fast  Group,  Chain  1 


9.6172E+0I  CM2 
9.5923E+01  CM2 
1.7292E-01  CM2 
9.6096E+01  CM2 
. OOOOE+OO  M/SEC 
2.1557E+04  M/SEC 
2.7659E+02 
4.3363E+02  BARNS 
4.6972E+01  BARNS 
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PLOT  OF  SPECTRUM,  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHARGY,  FOLLOWS 
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Figure  17b ;  Boraaaka,  5  Past  Groups,  Flux  riot 
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JOB  Old  BARNYARD  CHAIN  21 

THE  TOTAL  NUMBER  OF  BROAD  GROUPS  IS  3 
OUTPUT  WILL  BE  FOR  3  BROAD  GROUPS 
THF  BOUNDARY  FINE  GROUPS  ARE  3  6  11  12 


NUCLIDE  IS  CARBON 


ITS  NUMBER  DENSITY  IS  6.9940E-02  PER  8ARN-CM 

MICROSCOPIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW. 


scatter 

ELASTIC 

ELASTIC 

TOTAL 

FROM 

TO 

INELASTIC 

N-2N 

(PO) 

(PI ) 

TRANSFER 

1 

l 

.00808 

.00000 

1.30261 

.28939 

1.31069 

1 

2 

.01226 

.00000 

.41154 

-.1' 988 

.42379 

1 

3 

.00306 

.00000 

.00000 

.00000 

.00306 

1 

4 

.00000 

.00000 

.00000 

.00000 

.00000 

2 

2 

.00000 

.00000 

2.24744 

.38214 

2.24744 

2 

3 

.00000 

.00000 

.31718 

-.08037 

.31718 

2 

4 

.00000 

.00000 

.00000 

.00000 

.00000 

3 

3 

.00000 

.00000 

4.19494 

.29727 

4.19494 

3 

4 

.00000 

.00000 

.00000 

.00000 

.00000 

GROUP 

SIGTR  SIGR 

SIGA  NUSIGF 

SOURCE 

1 

1.62804  .42686 

. 

00000  .00000 

.34170 

2 

2.26286  .31718 

. 

00000  .00000 

.52670 

3 

3.89767  .00000 

. 

00000  .00000 

.13160 

NUCLIDE  I 
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Figure  17c:  Eomaska,  3  Fast  Groups,  Chain  2 
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Pigure  17 d*  homaaka,  g  Past  Groups,  Chain  2,  Continued 


V.  The  Calculation  of  Cross  Sections  and  Related  Constants 


In  this  section  the  equations  and  mathematical  expressions  used 
for  the  important  calculations  in  the  code  will  be  developed.  After 
a  brief  discussion  of  flux  wei Siting,  the  theoretical  methods  for  pre¬ 
dicting  the  energy  dependent  flux  will  be  described.  This  theory  con¬ 
sists  of  two  major  subdivisions,  that  for  epithermal,  or  above  thermal, 
energies,  and  that  for  thermal  energies.  Finally,  the  methods  used  to 
calculate  group  cross  sections,  both  epithermal  and  thermal,  will  be 
presented. 

Flux  Weighting 

The  microscopic  neutron  cross  section,  represented  by  the  symbol 
ct,  can  be  defined  ae:  the  probability,  measured  in  cross  sectional  area 
(barns),  that  a  neutron  will  interact  with  a  nucleus  in  a  specified 
manner.  The  macros cop  :  cross  section  E,  is  defined  as 

E  =  Na  (1) 

where  N  is  the  number  of  nuclei  per  cm3.  Thus,  the  macroscopic  cross 
section  is  a  measure  of  the  probability  that  a  neutron  will  interact 
with  a  nucleus  per  centimeter  of  path.  If  a  material  of  macroscopic 
cross  section  E  is  also  populated  by  a  neutron  density  n  (0nVcm3) 
all  with  constant  velocity  v  (cra/sec),  then  the  reaction  rate  R 
(events/cm3  -  sec)  is 

R  =  En-  (2) 

The  product  nv  in  equation  (2)  is  defined  as  the  neutron  flux, 
neutrons  per  cm2  per  sec. 

However,  the  neutron  cross  section  is  not  a  constant;  rather,  it 


73 


is  a  funct_on  of  energy.  Therefore;  the  reaction  rate  in  the  presence 
of  polyenergetic  neutrons  becomes 

R  =  \  E(E)cp(E)dE  (3) 

Jo 

where  cd(E)  is  the  energy  dependent  neutron  flux,  n(E)v(E),  and  where 
E  is  the  energy.  The  total  flux  is,  by  definition 

V  =  \  o(E)dE  (U) 

OO 

In  terms  o.  is  total  flux,  cp,  the  reaction  rate  can  be  written 

R  =  £cp  (£) 


where  £  is  some  average  macroscopic  cross  section  over  the  energy 
range  of  interest.  It  follows,  by  equating  the  right  hand  sides  of 
equations  (3)  and  (5),  that  this  average  cross  section  must  be  given 

ly 


f  £( E) cp  (E)dE 


£  = 


f  cp(E)dE  (6) 

Jo 

or,  in  other  words,  the  average  cross  section  is  a  flux  weighted 
average  of  the  energy  dependent  cross  section.  Of  course  a  similar 
expression  to  Bq  (6)  can  be  written  in  terms  of  the  microscopic  cross 
section,  i.e., 

J  J  CO 

V  0(E)cptE)dE 


a  = 


\  c,'(E)dE 
Jo 


(7) 


Equation  6  can  be  generalized  to  any  energy  range  (limits  on  the 
integral),  say  Ej  to  E?  in  order  to  produce  an  average  cross  section 
applicable  to  that  energy  range-  Such  averages  are  called  energy- 
group  cross  sections,  and  may  be  expressed  as 
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(8) 


where  £„  is  the  group  cross  section  of  the  nth  energy  group 
E“  is  the  lower  energy  boundary  of  the  nth  group 
E„  is  the  upper  energy  boundary  of  the  nth  group. 

Thai,  in  onden  to  detenmine  g\oup  cAotb  i ectionb  the  variation 
neutnon  fitux  uiith  eneAgy  mmt  be  knorn.  The  calculation  neceuaAy 
to  determine  the  energy  dependent  fitux  {an  exceeds  the  calculation 
neceMafiy  to  determine  gtioap  cao&6  6ectcon&  from  equation  [8]  above. 

In  this  code  all  of  Chain  I  is  devoted  to  the  calculation  of  energy 
dependent  flux.  The  shorter  Chain  II  utilizes  th^s  flux  to  calculate 
group  cross  sections. 

Ihe  flux  calculations  performed  in  Chain  I  are  in  solution  to  the 
energy  dependent,  Boltzmann  transport  equation  under  very  limiting 
assumptions  and  highly  idealized  boundary  conditions.  These  assump¬ 
tions  and  limitations  are  imposed  for  the  sole  purpose  of  obtaining 
a  usable  solution.  To  attempt  a  solution  for  energy  dependent  flux 
without  such  limitations  would  involve  (for  our  purposes)  a  prohibi¬ 
tive  amount  of  calculation.  In  addition  the  Boltzmann  equation  is  to 
be  applied,  in  this  code,  to  two  very  different  situations,  the  above 
thermal  energy  region  where  target  nuclei  kinetic  energies  or  veloci¬ 
ties  are  negligible  with  respect  to  the  neutron  velocities,  and  the 
thermal  region  where  the  velocities  of  the  target  nuclei  are  not  negli¬ 
gible  with  respect  to  the  neutron  velocities.  These  two  situations 
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require  different  assumptions  and  boundary  conditions  and  hence  dif¬ 
ferent  solutions  to  the  Boltzmann  equation.  Ad,  it  will  be  ieen  the 
above  thermal,  energy  dependent  $lux  i*  determined  by  the  method  o£ 
moment*  while  the  thermal  jlux  i *  determined  by  the  method  o£  Wigneri 
and  Wilkin*. 

In  either  case  we  start  with  the  steady  state  Boltzman. .  equation 
written  for  an  infinite,  non-multiplying  homogeneous  medium*  which 
scatters  and  absorbs  neutrons. 

u  cp(x,u,  u)  +  Et(u)  (C  (x,u,  u)  =  S(x,u,  a, 
dx 

+  ^dQ'^  du'  E,(u');p  (x,u'  ,u,)f(u'-ujue)  (9) 

where  x  is  spatial  position  of  the  neut) ons, 

u  is  the  neutron's  lethargy  or  logarithm  of  the  energy,  (Ref 
7:116), 

u  oru'  is  the  cosine  of  the  scalar  angle,  cosine  9,  between  the 
neutron's  direction  and  the  x  axis, 

U0  is  the  cosine  of  the  scalar  angle,  cosine  9,  through  which 
a  neutron  is  scattered, 

co(x,u,u)  is  neutron  flux  in  0nVcm2  -  sec-steradian-unit 
lethargy 

£t(u)  is  the  lethargy  dependent  macroscopic  total  cross  section, 

-l 

cm 

£«(u)  is  the  lethargy  dependent  macroscopic  scattering  cross 
section,  cm"1, 

dO'  is  the  differential  solid  angle  =  sin  9d9d0  =  -d(cos9)d0  =  -dud 9 

fCu'-ujUg)  is  the  scattering  frequency  function  or  probability 


*The  cross  sections  be  generated  by  this  code  will  frequently  be 
applied  to  media  which  are  either  multiply' ng,  hetrogeneous  or  both 
and  always  to  media  which  arc  finite.  Without  suitable  corrections 
this  can  obviously  lead  to  error. 
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that  a  neutron  with  initial  lethargy  u'  and  direction  u' 
will  when  scattered,  emerge  within  a  unit  lethargy  interval 
about  u  and  within  a  unit  solid  angle  about  u, 

S(x,u,u)  is  the  source  of  neutrons  per  uni  t  volume,  per  unit 
lethargy,  per  unit  solid  angle. 

The  scattering  frequency  function  and  the  macroscopic  scattering 
cross  section  under  the  integral  in  equation  ( 9 )  are  frequently  taken 
as  a  product 

£*(u/,  Jo)  =  SsCu')  f(u'  -•  u,  u0)  (10) 

w/iete  £s(u',u0)  16  denned  at,  the  mc/ioscopic  di^eAentiat  tautening 
CA066  Aectian  -in  cm  /t>teAaxLLan*  The  Boltzmann  equation  then  becomes, 

^  op 

u  "dx  (x'u,u)  +  Zt(u)qp  (x,u,n)  =  S(x,u,u) 

+  C  dfi'^»(u/,u0)cp(  x,u',u')du'  (11) 

The  two  directional  coordinates,  U,u'  and  the  cosine  scalar  angle  be¬ 
tween  them,  u0 ,  are  shown  in  Figure  18.  The  cosine  angle  Uq  is 
relateu  to  u  and  u'  by  the  law  of  cosines  for  a  spherical  triangle. 

Fast  Flux  Spectrum  and  Age 

In  the  energy  region  where  target  nuclei  motion  may  be  neglected 
(the  epithermal,  above -thermal  or  fast  region)  the  Boltzmann  equation 
can  be  most  easily  solved  by  assuming  that  the  source  is  isotropic  and 
consists  of  a  plane  of  infinite  area  located  at  the  coordinate  posi¬ 
tion  x  =  0,  i.e., 

S(x,u,u)  =  6(x)  (12) 


where  6  (x)  is  the  Dirac  delta  function  at  x  •*  0. 


y>'4. 


t 


Figure  18:  Neutron  Scattering  Diagram 

With  this  assumption  the  method  of  moments  may  be  applied  to 
find  the  fast  flux  as  a  function  of  energy.  However,  as  it  will  be 
seen,  the  minimum  calculation  necessary  for  the  flux  simultaneously 
produces  the  neutron  age  as  a  function  of  energy.  Since  this 
parameter  is  of  some  importance  in  elementary  calculations  it  is  also 
tabulated  by  the  code. 

The  appearance  of  the  cosine  of  the  scalar  scattering  angle  a, 
in  the  Boltzmann  equation,  equation  (9),  immediately  suggests  the  use 
of  Legendre  polynomials.  The  Legendre  polynomials  are 

F0(u)  =  1 
Pi(u)  =  u 

p8(u)  =  h  (3a2  -  1) 


.  i 


uim  in  general. 


Pa(u) 


m  . 

=  1_V  (-l)J(2n-  j)« 

'  2"/r0  j.'(n-2j)J(n-  j)» 


u 


n-a  i 


(13) 


where  m  =  ^  if  n  is  even  or  zero  or  ra  =  l|(n  - 1)  if  n  is  odd. 
These  polynomials  obey  the  orthonality  condition: 


$__p.(u)p,(U)dia-^6; 

where  6t  =  Kronecker  Delta  =  ^  n^. 

Other  properties  which  will  be  useful  are  (Ref  12:115) 

r2rr  r2TT 

\  Pn(uo)d0  =  \  Pa(u)Pa(u')d0 
uo  Jo 

where  0  is  the  plane  angle  shown  in  Figure  18,  and 

For  further  information  on  Legendre  polynomials,  the  reader  is  ref 
erred  to  Churchill  (Ref  3:200-203). 

The  flux  in  equation  (ll)  can  be  expanded  in  terms  of  these 
Legendre  polynomials  as 

CO 

^(x,u,u)  -V  — - -  Cpn(x,u’)Pa(u) 

n  L71 

n=0 

Similarily,  the  scattering  term  may  be  expressed  as 


(Hi) 


(15) 


(16) 


(i?) 


E.«u0)  -7  2m  +  l_s*(u/)Pm(Uo) 


m=0 


Utt 


(18) 


Note  that  both  the  indexes  m  and  n  are  required  since  the  integral 
term  in  equation  (11)  involves  both  the  flux  and  scattering  terms. 


i 
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Wj? 


By  using  the  orthoganality  condition,  equation  (lU),  it  is  also  easy 
to  show  that  co0(x,u)  is  the  all  angle  flux  and  cp,(x,u)  is  the  net  neu¬ 
tron  current  of  diffusion  theory.  Similarly  applying  the  orthogo¬ 
nality  condition  to  the  scattering  term  results  in 


S?(u')  »  E.(u') 

(19) 

Si(u')  =ul,(u') 

(20) 

where  Ef(u')  is  all-angle  total  scattering  cross  section  at  lethargy 
u',  and  u  is  the  average  cosine  of  the  angle  of  scatter  in  the  labora¬ 
tory  system. 

The  substitution  of  the  source  given  by  equation  (12)  and  of  the 
expansions  (17)  and  ( 18)  into  equation  (11)  yields 


n=0 


l 


2n+l  d  cp  n  (x,u) 

1|TT  ?jX 


UP„(u) 


2n+1  Et(u)  cpn(x,u)  P n 0-1 )  =  ^y^-6(x) 


Utt 


(  ao'  C  e.V>p.<U0)  j 

J  J°  m=0  w  J  L  Un  JJ 


(21) 

If  the  recurrence  relation  ( 16)  is  used  in  the  first  term  of  (21)  and 
if  the  integral  relationship  (15)  as  well  as  the  orthogonality  property 
( lit)  are  applied  to  the  integral  term  of  (21),  we  obtain 


Y  p..,*.;  ♦  „  p..t  w]* 

n=0 

(2n+l)  rt(u)cpa  (x,u)  P„  (u)  =  S(u)  6 (x) 
(2n+l)  ?n(u)  ^  du'  Ej(u')  con(x,u')  ]■ 


+ 


(22) 


Equation  (22), 


a  single  equation  with  an  infinite  number  of  terms,  may 


be  transformed  into  an  infinite  number  of  coupled  equations,  each  with 
a  finite  number  of  terms  by  employing  the  orthogonality  condition 
given  by  equation  (ll*).  This  transformation  is  carried  out  by  opera¬ 
ting  on  equation  (22),  term  by  term,  with 


C  +  1 

p^(u)(eq. (22)) du,  i  =  0,1,2,.... 


(23) 


The  results  are: 
SL  =  0: 


-~1/X;U)  +  co0(x,u)  =1  (u')  cp0(x,u,)du/  +  S0(u)  6(x) 

Ay  ’•'0  u 


dx 

1  =  1: 

2  ,  1  *<f>o  (x,u) 

3  dx  3  dx 

l  =  2: 


+  St  cPi(x,u)  =  V  li(u')  cp1(x,u/ )du7 

•Jo 


3^s(x,u)  ^2  9cp1(x,u)  .  „  ,  %  .. 

^  9x  Z  “  St  cps(x,u)  S,(u  )  cp^xju^du' 


(2U) 


If  equations  (21;)  are  truncated  by  assuming  that  cp.+1(x,u),  cp.  (x,u), 

...,  <Pm(x,u)  are  zero,  then  a  solution  of  the  remaining  finite  set  is 

possible  and  the  result  is  known  as  the  approximation.  The  index 

1  may  have  any  value,  e.  g.  P5 .  The  flux  in  such  an  approximation 
would  be  found  by  summing 

CD  (x,u,u)=£  2££Lc^(x,u)  P£(u)  (25) 

£=0 

However  in  this  case  we  are  not  so  much  interested  in  the  flux. 
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'*"*t  •  *  'fv  “* 


$(x,u,u),  as  we  are  in  the  volume,  angle  integral  of  the  flux,  i.e., 

■>'0/  (u)  -  \  dx  \  da-a(x,u,u)  (26) 

which  ij>  the  basic  energy  dependence  0($  the  ihxx,  which  we.  seek.  Equa¬ 
tion  (26)  is  merely  the  Zero  moment  of  the  spatial  distribution  of  the 
total  (all-angle)  flux.  The  higher  moments  would  be 

^  ^  (’-1)  =  ^  dx  ^  df)  x:(x,u,  u)  !  ^  dx  dQi)(x,u,  u) 

9  (u)  =  j  dx  j  dlx2  d(x,u,u)  /  j  dx  ^  dO^(x,u,  u) 


For  the  plane  source  which  we  have  assumed,  it  is  obvious  that  the  odd 
moments  S^1  \  9^,  ...  must  be  zero  in  value. 

Although  we  cannot  calculate  rc(x,u,  u)  exactly  (a  ?<a  calculation), 

we  can  calculate  the  moments  exactly  using  some  properties  of  the  Four- 
-ier  integral  transform  on  the  spatial  variable,  i.e.,  (Ref  18) 


°F  (x,u))  = 

JU 


p  +  OO 

^  ®  (x,u)  e“lpX  dx  =  9^(p,u) 


It  can  be  seen  that  if  equation  (28)  is  written  for  £  =  0  and  the 
special  case  p  =  0,  we  have  exactly  the  zero  moment  of  equation  (26), 


i.e. , 


9°  (u)  =  9  0(p,u)  J 

(Recall  that  v0(x,u)  =  ^  c(x,u,0)  cH.  ) 
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Thus  we  Fourier  transform  equations  (22*)  using  the  properties 
(Ref  18) 


r  3cp.(x,u) 

V  {  - - - }  =  ip  9(P»u) 


and 


Vf6(x)}  =  1 


f  ip  9f'^p>u)  +  4  ip  es(p,u)  +  Et(u)  e^p^) 


\  ^(u')  91(p,u/)  du' 


(30) 


to  obtain 

ip9i(p»u;  +  Et(u)  G0(p,u)  =  S(u)  Et(u' )0 0(p,u' )du/  (31) 


(32) 


§  ip  Mp,u)  +  I  ip  0  .i(pju)  +  Et(u)  0  a(Pju)  = 


-u 


rA  (u')  G^u'Jdu' 


(33) 


We  note  t.-nt  t.:e  transformed  Pa  equations  are  valid  for  any  ar¬ 
bitrary  p.  We  are,  however,  interested  in  the  specific  point  p  -  0. 
This  suggests  the  use  of  the  Maclaurin  expansion.  The  Maclaurin  ex¬ 
pansions  for  90  (p,u)  and  0!  (p,u)  are  chosen  as 


90(P>u)=  9°(u)  +  (-ip)  9q  (u)  *i±pr,  eg(n)  ♦  .... 


2J 


9i(p,u)=  9  °(u)  +  (-ip)  9  \(u)  +  —  0f(u)  +  .... 

2! 


(3k) 


(35) 


wnere  the  argument  of  the  Maclaurin  expan  si  on  is  (-ip)  rather  than 
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for  reasons  which  will  become  apparent  in  the  discussion  on  age  below. 
The  primes  denoting  differentiation  in  the  expa>  >n  have  been  re- 


placed  by  supercripts. 

That  is 

9o'(«)  5  9®(u) 

(36) 

Under  this  notation  we 

note  that  the  zero  moment  is  symbolized 

9°  =9°(u) 

(37) 

Davison  and  Sykes  (Ref  50^3)  show  that  the  n-th  moment  of  the 
neutron  flux  can  only  involve  spherical  harmonics  of  order  n  or  less. 
Further,  they  show  that  due  to  the  odd-even  nature  of  the  functions 
involved 

^  Xa0g(x,u)dx  -  0  (38) 

only  when  both  n  and  jf  are  both  even  or  both  odd.  Therefore,  in  the 
equations  for  the  second  moment  of  the  flux  (to  be  derived)  only 
,  6q,  ...  and  9* ,  9f,  ...  will  be  non  zero.  Thus  the  Maclaurin  ex¬ 
pansions  become 

90  (p,u)  =  9°  (u)  +  (u)  +  ...  (39) 

81  (p,u)  =  (-ip)  9j  (u)  +  — 9?  (u)  +  ...  (UO) 

By  substituting  equations  (39)  and  (UO)  into  the  transformed  Pn 
equations,  equations  Ol),  (32)  and  (33)  ana  equating  like  powers  of 
-ip, we  finally  obtain  the  moments  equations  which  are: 

(“ip)0  {  ^t(u)  9°(u)  =  S(u)  S?t(u')  ?°(u')  du'  )  (Ul) 

(-io)1  f  Zt(u)  9i(u)  -  +  \U  Ej( u')  8?(u')  du'f  u2) 

L  3  <3o  J 
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(-ip)2  {  Et(u)  e|(u)  =  2  0i  (u)  +  ^  Z°.(u')  6q  (u' )  du'  } 

(-ip)3  {  Et(u)  e?(u)  =  0Q  (u)  +  2  effu)  +  £*  E;(u')  0?(u')  du'  }  (U») 


I-t  muA-t  be  emphasized  that  these  moments  equations  oJie  fiigoKous.  The 
spheniaal  harmonic  expansion  was  quite  general  and  not  terminated  at 
some  finite  l. 

Equation  (1*1)  alone  is  sufficient  to  calculate  the  zero  moment, 

0q  (u)  (which,  again,  is  the  volume  angle  integral  of  the  energy  depen¬ 
dent  flux)j  equations  (1*1  and  (1*2)  alone  are  sufficient  to  calculate 
0^(u);  and  equations  (1*1),  (1*2)  and  (1*3)  alone  are  sufficient  to  cal¬ 
culate  0q  (u) .  This  permits  a  direct  determination  of  neutron  age. 
Neutron  age  in  a  slab  geometry  is  defined  as  one  half  of  the  second 
spatial  moment  of  the  flux, 

t  =  (1*5) 

But  the  second  spatial  moment  is  given  by  equation  (27)  which  can  also 
be  expressed  as 


0(s)(u) 


/%'  /-» 

\  dx  \  dfl  x3cp(x,u,u)  \  x2  cp0(x,u)  dx 

W  “00  —  CD 

CO  O  p  +  00 

\  dx  \  dflq)(x,u,u)  \  <p0(x,u)  dx 

O  -  oo  _  oo 


S -  {  X3  Cp0(x,u)  }p  =  o 
^  l  cp0(x,u)}PM0 
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PIW^  JPWI",P.WP» 


jefapv^aiU  m  g « .11  * wmwf  lliw piu  mjwj . 


Since 


Vu*  j£(x,u))  ■  i>  £!!*&>?) 

V 

(another  Fourier  transform  property,  Ref  18),  it  follows  that 


<U7) 


T  =  3g  0  (2)  =  3g 


d2  90  (P»u) 

bp2 


■1. 


P  =  0 


0  ,  i2[(-i)a  9|(u)]  ,  e*(u) 

«  -  =  ^ 


[0O  (p.u)]p  x„  [e"(u)]  9q  (u) 

So  that  a  simultaneous  solution  of  (Ul)>  (1*2)  and  (1*3)  yields  both 
the  desired  energy  dependent  flux  [9£(u)]  as  well  as  neutron  age 
[%0§(u)/es(u)].  it  is  important  to  note  that  these  quantities  were 
determined  without  truncating  equations  (31),  (32).... 

Recall  that  the  purpose  here  is  to  develop  a  few-group  machine 
program  to  generate  group  cross  sections  and  other  related  reactor 
constants.  Therefore,  tne  moments  equations  must  be  put  into  multi¬ 
group  notation,  and  integrated  over  a  specific  lethargy  range  u^  to 
Consider  the  zero  moment  equation,  Eq  (Ul)  In  group  form  it  is 


+ 

Ua 


(*U  n  (»u  n 

Ju-  Et(u)9o(u)  du  =  V  -  S(u)  du  + 

un  Un 


\  .  du\  r°  (u')  9°  (u7 )  du' 
Jun  Jo 


(U9) 


From  the  definition  of  a  flux  weignted  group  cross  section,  Eq  (8),  it 
is  easy  to  see  that  the  first  term  becomes 


fill  t 


\  Et(u)  9°(u)  du  =  Stn  9°  n  (50) 

**n 

where  Etn  is  the  macroscopic  total  cross  section  of  the  n-th  group  and 
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is  tiie  total  group  flux.  (Remember  that  in  this  case  tne  trans¬ 
formed  flux  has  the  property  of  being  equal  to  the  spatially  integra¬ 
ted  all  angle  flux  in  real  space.)  Similarly  the  second  term  is 


+ 

Un 


S(u)  du  =  X, 


(5D 


where  y ,,  is  the  fraction  of  the  source  emitted  in  the  n-th  group  pro¬ 
viding  S(u)  is  a  normalized  source.  The  last  term 


,-.un 

\  _  du  \  Z°,  (u')  flo(u')  du' 
Ju,  Oo 


(52) 


can  be  considered  as  a  group  of  n  double  integrals  where  the  integral 
from  o  to  u  on  du'  has  been  broken  into  n  interval,  Au(,  Au^,  Au'-,, ... 
It  is  important  that  the  integration  be  carried  out  on  the  variable  u 
first  since  tne  limits  of  integration  on  u,  i.e.,  u„  to  u(( ,  are  ac¬ 
tually  functions  of  u'  and  the  maximum  Au  of  the  scattering  nucleus. 
Thus  the  last  term  becomes 


L;  duL;  r°(u,)  du'  *  1;  d“L^  °s  fu')du' 

...+  ^  du  ^  E°(u')  du'  + 

-'u;  ^Au'n_! 


f:  du\A  ,  S?(u')  00  (O  du'  (53) 

dun  JAun 

Now,  each  double  integral,  say  the  j-th,  is  a  measure  of  the  neutron 
scatter  transport  from  the  j-th  to  the  n-th  group.  Thus  this  last 
term,  B?  (53),  can  be  written  as 


+ 


O  aO 

-ll-n  01 


v-o 

—  i  2-  a 


eo  a 


+ 


v 


o 

s 


:-n 


9 


o 

0  1 


j=l 


(5u) 


Substituting  Bq  (50), 
terms  yields 

(St,  ~  - 


(5i)  and  (5U)  into  Eq  (1*9)  and  rearranging 


n-1 

.,-.K.  =  Vn  +  Y  E?j-a  9°  j 

Lj 

j=l 


or 


n-l 

Sr.C  ~  x  n  +/  sSi-.  e°; 
j=I 


(55) 


where  Era  is  the  macroscopic  total  removal  cross  section  from  the  n-th 
group,  that  is  scatter  removal  plus  absorption  removal.  Note  that 
this  equation  is  essentially  a  neutron  balance.  Treating  the  first 
roment  equation  similarly  yields 


(St,  -  s1.,-,)  e,\ 


n-l 


J-  n 


(56) 


Equation  (56)  is  the  neart  of  the  difi  rence  between  the  multigroup 
moments  method  and  tre  age  diffusion  methods.  Note  cnat  the  dif— 
feren.  c  term  multiplying  8^  a  on  the  left  is  not  tue  common  transport 
cross  section,  in  that  it  only  subtracts  away  the  within  group  scat¬ 
tering  term.  Rather  than  assume  that  collision  densities,  (£*%),  a r 
slowly  varying,  this  method  explicitly  accounts  for  the  linearly 
anisotropic  term  u  I,  .  Using  the  term 

^Se'.n  =  (S  t  n  ~  £*a-a)  (5?) 


Eq  (56)  becomes 


y  p3 

“!{tn  oa 


P° 

■'On 


n-l 

)  *—  5  J—  a  e,1 . 

j=l 


(58) 
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Finally  the  second  moment  equation  can  be  v/rit/ben 


n-1 


§cn  2  9j  a  +  Y  £sj-n  9^, 


(59) 


j=l 


The  group  age,  t„,  may  be  defined  in  a  similar  manner, 


Tn  = 


Jln  9°^  du 


4nGo(u)  du 


(60) 


However,  group  age  must  be  more  carefully  interpreted  than  group 
cross  sections.  The  only  average  age  whicn  has  any  meaning  is  the 
average  age  of  neutrons. which  leave  tne  group.  Therefore,  t8  ,  is 
always  to  be  associated  with  the  lower  energy  boundary  of  the  group. 
By  equation  (U8)  group  age  may  be  written  as 


l9|n 


*  2  9° 


(61) 


OR 


This  may  be  rewritten  by  dividing  equation  (59)  by  2Ern9°n  to  yield 

n-1 


zrs  9, 


On 


j=l 


(62) 


But,  9^  can  be  expressed  as 

-  2’,  OS  j 

Therefore,  the  group  age  expression  becomes 

n-1 


(63) 


T„  - 


I _  fgl  +  V  S0  .1 

r  po  |blE  /  9°-'  ’  -1  I 

t,rn  ‘'on 

Finally  by  defining  ar  as 


(620 


.1=1 


^n  “  N  arn  0°n 


(65) 


89 


4 


a 


n 


=  Xn 


n-1 


j=i 


qi  J~  n 
°r  1 


a  i 


(66) 


01  n 


a, 


n-1 

—  y  rt,..  e;, 


y*  ictn  ^rn  ^  ^  i|tD  j =1 

,  n-l  n 

9in  1  V  a'J-“ 

T  =  - +  —  ) - a4  t j 

On  1»  <7rn 


(67) 


(68) 


where  N  =  total  nuclei  density,  nuc]ei/cm3  x  10~24  ,  and  the  o's  are 
the  respective  microscopic  cross  sections  in  barns  with  the  subscript 
notation  the  same  as  in  previous  equations.  The  group  fluxes,  the 
and  the  group  ages  can  be  determined  by  solving  tne  above  equations. 
Cross  section  data  is  read  in  from  a  cross  section  library,  which  will 
be  discussed  later,  and  calculation  is  begun  with  group  one  where  the 
scatter  in  term  is  zero.  Calculation  proceeds  consecutively  through 
all  remaining  groups  down  to  lethargy  17  (O.hUi  ev).  All  group  fluxes 
are  then  normalized  to  0(16)  =  1  in  order  to  be  normalized  to  the 
thermal  flux  at  the  same  value.  It  should  be  noted  however  that  the 
input  cross  sections  are  generated  from  a  modified  version  of  GAM 
(Ref  13)  which  does  not  include  the  effects  of  resonance  ibsorption. 
Therefore,  resonance  absorption  must  be  accounted  for  in  the  code. 

This  is  done  by  calculating  the  resonance  escape  probability. 


Resonance 

The  resonance  escape  probability  is  the  probability  that  a  neu¬ 
tron  will  escape  resonance  capture  during  slowing  down  and  is  given  by 
Murray  as  (Ref  16:62) 
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(69) 


p  =  exp 


B, 


ttVS 


dE] 
•»  eJ 


>E  E- 

where  p  is  tee  resonance  escape  probability 

E0  is  the  reference  energy,  usually  1C  Mev 
§  is  the  average  logarithmic  r.iergy  decrement  per  collision 
Z,  is  the  macroscopic  absorption  cross  section,  cm*1 
Z t  is  the  macroscopic  total  cross  section,  cm-1 
Murray  also  points  out  that  the  exponent  in  Eq  (69)  may  be  written  as 


h.  r»^o 

\  a  • 


x  +  N 


±2jL 


dE 

E 


(70) 


where  If  is  the  nuclide  density  of  the  resonance  absorber  x  10~24  oer 


cm 


a*"  is  the  microscopic  absorption  cross  section  of  the  resonant 
absorber,  barns 

and  Ss  is  the  macroscopic  scattering  cross  section,  cm-1 
Note  that  the  term  in  the  integral  (70)  depends  on  the  total  scatter¬ 
ing  cro'-s  section  per  absorber  atom  Z «/N \  which  is  a  function  of  the 
moderator- fuel  ratio.  The  integral  term  in  (70)  above,  defined  as  the 
resonance  integral,  is  usually  approximated  oy  empirical  expressions. 
Isbun  \fief  ll:i;f>9)  lists,  some  of  these  expressions: 

(for  U83e) 

n 

171 


(PJ).ff  -  2.69  (— )°'4' 
V  JJS36  J 


0  s  “Si 5  k000 

n238 


(71) 


/n(RI).ff  =  5.6U  - 


163 


(Z,/N338)0  *65  N 


>  UCOO 


338 


(72) 
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and  (Rl)eff  =  260  for  E,/H?3e  =  », 
and  for  Th232 

(RI).„  -  6.33  W'  -50°  '  m 

(Rl).ff  =  70  J±2-  '■  U50C  (710 

where  (RI)«ff  is  the  effective  resonance  integral  in  barns. 

E,  is  the  macroscopic  scattering  cross  section  of  tne  mixture 
over  the  resonance  region,  cm-1. 

The  resonance  escape  probaoility  given  in  equation  (69)  then  becomes 

P  =  exp  [-  (Rl)ef ;  N A/?  S,]  (75) 

where  the  is  for  the  mixture. 

The  effect  of  resonance  absorption  is  then  introduced  Into  equa¬ 
tions  (55),  (56)  and  (59),  the  moments  equations,  by  increasing  the 
ar9  and  orlo  (resonance  takes  place  w:  groups  9  and  10).  tnru  the 
equations  (Ref  '  :l/9) 

Or,  =  ara  -  l[V  (76) 

where  a”*  s  A,(fil).f  f /iu  (7?) 

where  ora  =  microscopic  removal  cross  section  of  croup  n,  bams 

A s  =  an  empirical  factor  to  apportion  the  resonance  absorption 
between  groups  9  and  Lu. 

Au  =  the  lethargy  interval  over  which  resonance  is  effective. 
Murray  gives  a  value  of  5.5. 

This  has  the  effect  of  reducing  the  downscatter  to  groups  10  and  11  by 
appropriate  factors  without  violating  neutron  conservation  within  the 
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groups.  The  code  allows  the  user  the  option  of  providing  his  own  value 
of  (Ri)eff  should  he  so  desire. 

The  equations  for  determining  the  fast  flux  spectrum  and  age  are 
now  complete.  The  equations  specifying  the  thermal  flux  spectrum  are 
developed  in  the  next  section. 

Thermal  Flux  Spectrum 

In  elementary  reactor  physics  the  thermal  flux  spectrum  is  usu¬ 
ally  assumed  to  be  in  thermal  equilibrium  with  the  moderator  producing 
a  Maxwellian  distribution.  This  is  not  true  however  since  absorption 
removes  neutrons  and  thus  disturbs  the  equilibrium.  On  the  other  hand 
the  moments  equation  used  for  epithermal  energies  can  no  longer  be 
used  since  they  take  into  account  only  downscatter,  whereas  the  possi¬ 
bility  of  upscatter  must  be  considered  when  the  energy  of  the  modera- 
ter  nuclei  is  comparable  to  the  neutron  energy. 

Wilkins  (Ref  20)  derived  a  second  order  differential  equation 
which  describes  the  thermal  neutron  spectrum  in  an  infinitely  heavy 
Maxwellian  gas  and  allows  absorption.  The  thermal  spectrum  generated 
by  this  code  is  calculated  by  solving  the  Wilkins  equation.  This  solu¬ 
tion  is  expressed  in  terms  of  a  cross  section  determined  parametex-, 
the  absorption  parameter,  which  in  turn  (in  the  case  of  a  mixture  of 
nuclides)  is  expressed  in  terms  of  the  thermal  utilization,  another 
cress  section  oarameter.  These  two  parameters  are  developed  next, 
followed  by  a  development  of  the  Wilkin's  equation. 

Thermal  Utilization 

The  thermal  utilization  f ,  is  defined  as  the  ratio  of  thermal  neu¬ 
trons  absorbed  in  fuel  to  all  thermal  neutrons  absorbed  or, 
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Cv f  £.],a, 


L  vf 


V  1 
-*-‘2  0  d 


[vf  ' ll »  «  ! 


[vf  Sj,od 
[vf  e, 


t  7R^ 

\  >  '■i 


wiiere  vf  is  the  volume  fraction.  However,  the  thermal  neutron  flux  in 
the  fuel  elements  of  a  heterogenous  reactor  will  be  depressed  due  to 
the  strong  absorption  in  the  fuel.  This  effect  complicates  the  cal¬ 
culation  of  the  thermal  utilization.  The  proper  expression  in  the 
case  of  heterogeneous  reactors  is  (Ref  l6:6?-8u) 


1 

J. 


£ 


1 


+  Lz£.  d  f  + 

L  vf  T.  Jf  B  e  i 


(2-1) 


where  F  =  K0w  coth  (K0w)  (in  the  case  of  slab  geometry) 

E  =  Kj(W-w)  coth  lKj(W-w)j  (For  cylindrical  geometry  see  Ref 

16:87-88) 

Ko  =  £</D0 
Kf  =  Ii/Da 

W  =  half  the  distance  between  fuel  element  midpoints 
w  -  the  distance  from  the  midpoint  of  the  fuel  element  to  tne 
boundary  of  the  fuel  element 

and  the  subscripts  0  and  1  refer  tc  fuel  and  moderator  respectively. 

In  effect  F  is  actually  the  ratio  of  the  flux  at  tne  fuel  surface  to 
the  average  flux  in  the  fuel  while  (E-l)  accounts  for  the  fact  that 
tne  flux  in  the  moderator  rises  from  a  minimum  at  the  fuel  element  to 
a  maximum  midway  between  elements. 

In  the  code  tne  F  and  E  factors  must  be  read  in  for  heterogeneous 
cell  calculations  but,  they  are  assumed  equal  to  unity  for  homogeneous 
systems.  In  the  case  of  heterogeneous  reactors  in  which  thin  ~iosely 
spaced  fuel  plates  are  used  tne  system  can  be  assumed  homogeneous  from 
a  calculational  point  of  view. 
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The  AuSOiuti Oi"i  Pai'diutjtei*  is  aeiineu  **c 


2Aa,(E)/E/kT  (80) 

os(E) 

where  A  =  nuclide  mass,  amu 

k  =  Boltzmann's  constant  in  MeV/(°K)4 
T  =  temperature,  deg.  Kelvin 

Since  A  is  used  only  for  thermal  calculations,  <J  ,(E)  is  assumed  con¬ 
stant  and  a»(E)  is  assumed  to  vary  as  l/v.  Then 

2Ao  »0 


where  o40  =  microscopic  absorption  cross  section  at  energy  kT  in  barns. 
Weinberg  and  Wigner  (Ref  19:337)  point  out  that  for  a  mixture  of  nuc¬ 
lides  the  absorption  parameter  should  be  divided  by  the  factor  (1-f). 
Therefore,  to  find  the  absorption  parameter  for  a  mixture  of  nuclides 
the  following  equation  is  used: 


,  I 

1  i _ 

(1-f) 


f 


where  2^  =  total  macroscopic  cross  section  of  the  i-th  nuclide,  and 
tlie  summing  variable  i,  refers  to  all  moderator  nuclides. 
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The  Wilkins  Equation 


The  Boltzmann  transport  equation  for  steady  state  conditions  in 
an  infinite  homogeneous  medium  was  given  in  equation  (8).  It  is  re¬ 
peated  here  for  convenience. 


u 


-^—  ro(x,u,u)  +  ''t  (u)  i'{x,u,u)  S(x,u,u) 
3  X 


d.V  \  F.,( u' )  ?(x,u' ,u')  f(u'  -  u;  u0)du' 


(“) 


However  in  the  thermal  case  we  shall  taxe  the  source,  S(x,u,u)  as 
zero,  and  introduce  the  neutrons  as  a  boundary  condition  by  specify¬ 
ing  some  arbitrary  number  of  neutrons  slowing  down  past  the  upper 
energy  limit  of  the  tiiermal  group.  In  the  absence  of  a  source  and  be¬ 
cause  the  media  is  infinite  the  variation  of  flux  with  fotn  spatial 
position  and  direction  disappears.  That  is  0(x,u,u)  =  0(u)  and 

3  0 (x,u,u).  _  ^  Thus  the  Boltzmann  equation  becomes 
3  x 


[£»(u)  +  £s(u)]  P(u)  =  C  F,(u'~  u)  co(u')du 

v)  —  00 

or  in  terms  of  energy 

(X(E)  +  l.i(E)]  s(E)  =  \  E.(E'  -  E)  =>(E')dE' 

«Jo 


(83) 


(8U) 


where  the  limits  on  the  integral  are  correctly  o  to  ®  because  of  the 
the  possibility  of  upscatter  as  well  as  downscatter. 


where 


E.(E-  E')dE' 


r>  /  /  \  ;  !  \  %  f 

^,vu  )  i(u  -  u)au 


(63) 


and  E,(E -  E' )dE'/Es(E)  is  tne  probability  that  a  neutron  at  energy  E, 
will  be  scattered  into  the  energy  range  E'  to  E'  +  dE' . 

Both  Hurwitz  (Ref  10)  and  Wilkins  (Ref  20)  solve  this  equation 
for  a  heavy  gaseous  moderator.  Hurwitz1 s  solution  is  more  easily 
followed  and  his  method  will  be  used  to  develop  the  Wilkins  equation. 

According  to  the  principle  of  detailed  balance  for  neutrons,  in 
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(86) 


the  case  of  zero  absorption 

S(E#  -  E)  M(E')  =  >1(E-  E' )  M(E) 
where  M(E)  =  tne  Maxwellian  flux  at  energy  E.  As  Hurwitz  states, 

"Tuis  [is]  a  property  of  tne  kernel  which  is  independent  of  the  flux". 
Using  this,  equation  ( 8U)  can  be  rewritten  in  terms  of  the  new  variable 

Y(E)  =co(E)/M(E)  (87) 

where  y(E)  is  tne  deviation  from  a  Maxwellian  flux.  Thus  equation 
(8U)  becomes 

S.(E)  V(E)=C“z(E-E/)  [Y(E')  -  Y (E)jdE'  (88) 

Jo 

Hurtitz  then  uses  a  Taylor  series  expansion  of  T(E)  about  the  point 
E'  =  E  to  give 


[-rrr]T(E)  - 

L  r..  e  J  At?  2 


dE  2 
(Ref.  10:283) 


AE 


,2  d3y 


dEJ 


Ss(E) 

where  AE  =  2(2kT-  E)/A 
Ap  =  )4EkT/A 
and  A  -  mass  number  of  the  moderator  nuclei. 
Substituting  for  AE  and  AE2  in  equation  (89)  yields 


(89) 


EkT  +  (2kT -  E)  —  -  -y 
dE2  dE  k 


0 


(90) 


where  A  is  the  absorption  parameter  defined  by  equation  (8l). 

It  is  convenient  at  this  point  to  introduce  a  new  variable  x 
which  is  a  dimensionless  velocity  defined  by 


x  =  J  E/kT 

In  terms  of  this  new  variable 

x2N(x )  =  2  Eco(E)  =  2cd(u) 


(91) 


(92) 
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and  Y(E)  =  Cx~3  axp  (x2)  N(x)  (91 

where  N(x)  =  number  of  neutrons  at  velocity  x  f  and  C  is  a  constant. 


Now  noting  that 


Y(E)  =  f(x) 

Y'(E)  =  f'(x)  % 
ob 


y"(E)  =  f"(x)  f—  i8  f'(x)  [— ] 

L  dE  J  i  j  L  dE  J 


where 


dx  =  _1_  ^d  £*£  =  -  1 

dE  2xkT  dE3  it(kT)2x; 


and  where  the  primes  denote  differentiation  with  respect  to  tne  res¬ 
pective  independent  variables. 

We  find  that 

exp  (x3) 

*'(E)  ,  -  r  x-3  N/(x)  +  (2x-s  _  2x-4)  n(x),  (98) 

2  kT  J 

»"( E)  ■  |  x~>  N-/W  „  (kx-=  .  N<(x)  , 

U(kT) 2  1 

(8x-6  -  8x-4  +  Ux~2)  N(x)  |  (99) 

Substituting  Eq  (95)  and  (96)  into  Eq  (87)  results  in  the  well  known 


Wilkins  equation 


xN"(x)  +  (2x3  -  l)N'(x)  +  (ipc-A)N(x)  =  0  (100) 


This  equation  is  solved  by  setting 


N(x)  =  x2  exp  (-  x2)M(x) 


which  converts  Eq  (100)  into 


xM"(x)  +  (3-2x2)M'(x)  -  AM(x)  =  0 


(101) 


(102) 


This  equation  is  solved  by  the  power  series  method.  The  method  is 
straightforward  and  for  an  expansion  about  the  point  x  =  0  results  in 


,.ii.  ii  ■fiilii  lSW|PtiUlL,«'PliUU 


M(x)  =  £  r.n  xn 
n=0 


(103) 


A 

where  aj  "  a0  j 
A  2 

a3  -  a0  2if 

an  ~  n(n+2)  C2(n  -  2)  an_2  +  A  an_x] 

and  a3  is  arbitrarily  set  equal  to  1.  Equation  (102)  can  also  be  ex¬ 
panded  about  the  point  x  =  a.  In  this  case  by  letting  x  -  a  =  V,  Eq 
(102)  is  converted  to 

(V  +  a)  M"(x)  +  (3-2V2  -  l;Va-2a3)M/(x) -AM(x)  =  0  (101*) 


and  the  solution  is 

M(x)  = 


Y  bn(x  -  a)n 
n=0 


(105) 


where  ^ 

b»  =  - - —  [bn-1(n-l)(n-l-2a2)  - 

a(l  -  n)n 

ba_2  {2ja(n-2)  +  A}  -  bn_3  2(n-3)] 

In  the  code  the  values  M(x),  M'(x),  N(x)  and  x^x),  (which  is  « 
°^))  are  calculated  at  every  0.25  lethargy  intervals  beginning  at 
u  -  23  and  continuing  to  u  =  16.  The  calculation  employs  Eq  (103) 
for  the  first  point  then  uses  Eq  (105)  for  all  subsequent  points, 
the  point  x  =  a  used  in  Eq  /105)  being  the  point  x  in  the  previous 
calculation.  Terms  are  added  until  relative  contributions  of  addi¬ 
tional  terms  to  M(x)  and  V' (x)  are  less  than  10”6  relative  value.  All 
thermal  fluxes  are  normalized  to  0(16)  =  1  to  match  the  epithermal 
fluxes  from  the  moments  calculation.  To  arrange  for  the  proper  in¬ 
crementing  of  the  variable  x,  such  that  Ax  corresponds  to  0.25  leth- 


argy  units,  the  following  is  used: 

*> 

x(u=23)  =  exp  [(-  u.^do^/B^iery8] 

xs  =  Xi  +  Ax 

Ax  =  xx  [ exp  (.12$)  -  l] 

T  is  temperature  in  °  Kelvin 
Thermal  Constants 


U06) 

(107) 

(108) 


In  the  code  the  thermal  range  is  arbitrarily  defined  as  0  to  1.2$ 
av  which  is  00  to  16  on  the  lethargy  scale.  All  ihex-mal  constants  are 
calculated  for  this  interval. 

Most  Probable  Neutron  Velocity 

The  most  probable  neutron  thermal  neutron  velocity  is  found  by 
locating  the  peak  of  the  N(x)  vs  x  curve  by  the  code.  However,  since 
values  of  N(x)  are  found  for  0.25  lethargy  units,  the  value  of  vp  is 
not  exact.  The  true  probable  velocity  is 

vp  =  [X(16,500  (10?) 

(T  in  °K) 

Average  Neutron  Velocity 

The  average  neutron  velocity  is  found  by  Sq  (103)  except  that  the 
equation  is  evaluated  at  x  =  x,  which  is  the  average  x,  rather  than  at 
x  *  xp.  The  average  normalized  neutron  velocity  x  is 

,-x(u  =  iS) 


X 


x  N(x)  dx 


x  = 


rx(u-i«) 

^  N(x)  dx 


(no) 


The  integral  over  N(x)  can  be  found  by  integrati.'.j  Bq  (100) 

xN'(x)  +  2(xs  -  1)  N(x) 


X  N(t) 


dt  = 


(HI) 


100 


wutjrc 

N'(x)  =  2  [x°  exp  (-  x2)  M(x)]  {  -  -  l}  + 

L  x 

x2  exp  (-  x2)  M'(x)  (112) 

Thermal  Diffusion  Length  Squared 

The  thermal  diffusion  length  squared  is  defined  af  one  sixth  of 
the  mean  square  distance  a  neutron  travels  from  the  tine  it  becomes 

— 4 

thermal  until  capture.  The  total  distance  r,  a  neutron  travels  in 
making  n  collisions  is 

A  A  A  A 

r  »  7  i  9X  +  ra  8S  +  r3  03  +  ...  +  ra  0a  (113) 

where  rt  8 (  is  the  magnitude  and  direction  vector  of  the  neutron  be¬ 
tween  collisions.  The  total  mean  square  path  is  the  scalar,  or  dot, 
product  of  equation  (113)  with  itself.  Note  that  the  cross  product 
terms  will  involve  the  dot  products  of  the  various  neutron  direction 

A  A  A  A 

vectors  9j  ,  and  9t  •  9j  =  cos  9U  where  is  the  angle  between  9. 

8 j .  But  for  isotropic  scatter  the  angle  9U  between  two  successive 

paths  i  and  j  is  completely  random.  Thus,  there  are  as  many  positive 

cosines  as  negative  ones,  and  on  the  average, cancel  each  other  out. 

So  _  K 

r2  =  £  rj?  (lUi) 

k=l 


where  k  is  the  collision  number  and  K  is  the  total  number  of  collisions. 
It  can  be  shown  that  any  r2  is  twice  the  mean  free  path  squared  for 
that  collision  (16:276),  therefore 


where  X  is  the  average  neutron  mean  free  path  over  the  thermal 


ranee. 


is  the  total  average  thermal  macroscopic  total  cross  section 


Sicas  (Ref  17:lii-l6)  shows  that  n,  the  number  of  collisions  a  neutron 
makes  is  equal  to  tne  ratio  of  the  scattering  rate  to  the  absorption 
rate.  Hence,  equation  (115)  becomes 


pl*  125  ev 

2  \  Et0  (E)  dE 


3#  125  o  *f 


(116) 


Ea(E)  0(E)  dE 


Now,  the  thermal  diffusion  length  squared,  L3,  is  one-sixth  this  value. 
The  quantity  Et  can  be  found  by  assuming  that  E,  is  constant  and  E, 
has  a  l/v  dependence  in  the  thermal  range.  Since  x  =  /g/kT,  L3  can 


be  written  as 


E*  \  x  N(x)  dx 
o 

3  s3  s.0  y 


(117) 


N(x)  dx 


L3  = 


3  E3  Et, 


(118) 


where  xQ  =  x  corresponding  tc  1.125  ev, 

E»o  =  Cx  Ea(x)]x=5Cp,  the  most  probable  value.  Howe cer,  E  is 
itself  an  average  and  can  be  rewritten  as 


E  -  E.  +  E. 


(119) 


where 
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(120) 


r 

Et(x)  XN(x)  dx 

- 

\  x  N(x)  dx 
°o 

Hence  the  equation  for  L2  can  finally  be  written  as  (Ref  17:J2) 


L2  = 


a*  x 


3  a.0  E(a.  +  a,0(x)  N]2 


(121) 


where  N  is  again  the  nuclei  density  x  10“24  cm-3 

cr  is  the  microscopic  cross  section  with  subscripts  t&->  same  as 
before  in  bams 

x  is  the  average  normalized  neutron  velocity  and  the  energy 
range  for  the  calculations  is  understood  to  range  from  0  to 
1.125  ev  in  the  code. 


Migration  Area. 

The  total  migration  area  is 


M2  =  L2  +  t 


Th 


•'122) 


where  is  the  age  to  thermal. 


In  the  code  the  migration  area  is  found  by  adding  the  moments 
methods  age  to  lethargy  16  (1.125  ev)  to  the  value  of  L2  computed  from 
lethargy  16  to 

Thermax  Cross  Sections 

Ihe  2200  tn/sec  cross  sections  for  the  appropriate  nuclides  are 
input  data  to  the  code  from  cross  section  libraries.  However,  these 
cross  sections  must  be  adjusted  to  compensate  for  temperatures  above 
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I  braries. 

|  Actually,  an  iterative  procedure  should  be  used.  In  Eq  (123), 

j  a40  =  [x  a,(x)]x=Xp  .  However,  the  thermal,  cross  sections  read  into 

the  program  are  2200  meter/sec  cross  section  values.  They  are  then 
adjusted  for  temperature,  assuming  that  the  flux  is  Maxwellian.  But, 
when  the  flux  is  not  Maxwellian  the  most  probable  thermal  neutron 
velocity  xp  is  not  equal  to  xp  for  a  Maxwell-Boltzmann  distribution. 
Therefore,  tne  value  of  xp  used  in  calculating  o»0  may  not  be  correct, 
i  The  correct  procedure  would  be: 

I 

1  (l)  adjust  the  absorption  cross  section  using  equation  (123)  at 

i 

i 

■  the  end  of  the  thermal  calculations. 

! 

(2)  re- input  the  value  of  a,  calculated  by  equation  (123). 


O'1  recalculate  the  thermal  spectrum  using  the  new  value  of  a,, 
(it)  return  to  step  (1)  above  until  a,  remains  constant. 

To  keep  computer  run  time  to  a  minimum  the  iterative  procedure 
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shown  above  is  not  used. 

Cross  Section  Collapsing 

The  theory  presented  in  the  previous  paragraphs  is  all  incorpora¬ 
ted  in  Chain  I  of  OLD  BARNYARD  to  produce  flux  as  a  function  of 
energy,  both  epithermal  as  well  as  thermal,  in  an  infinite  arbitrary 
mixture  of  nuclides.  This  energy  dependent  flux  is  expressed  as  a 
histogram  of  11  epithermal  values  and  9  thermal  values.  These  flux 

r 

values  are  employed  in  Chain  II  to  calculate  flux  weighted  average 
cross  sections  over  any  number  of  epithermal  groups  from  1  to  11  and 
a  thermal  group.  The  exact  number  of  epithermal  groups,  called  broad 
groups  to  distinguish  them  from  the  input  fine  groups,  are  at  the 
specification  of  the  user.  The  boundaries  of  the  broad  groups  are 
restricted  to  values  which  are  boundaries  of  the  fine  groups. 

Broad  group  cross  sections  are  calculated  from  the  formula 


BG  i 


CT 


I*' 


(125) 


BG 

where  on  is  the  microscopic  cross  section  of  n-th  kind  for  a  broad 


group 

<Pi  is  the  total  flux  of  the  i-in  fine  group 

aj{  is  the  microscopic  cross  section  of  the  n-th  kind  for  the 

i-th  fine  group. 

The  interrelationships  which  exist  among  some  of  the  cross  sec¬ 
tions  output  by  the  code  is  not  always  obvious.  The  following  expla- 

!  nation,  therefore,  may  be  useful. 

\ 

t 

1 

I 

i 
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Total  Transfer 

The  group  to  group  total  transfer  cross  section  is  that  value 
which  when  multiplied  by  the  losing  group  flux  gives  the  rate  of  trans¬ 
fer  from  j  to  k  and  is  the  sum  of  inelastic,  elastic  (PO)  and  twice 
the  n  -  2n  cross  sections  from  j  to  k.  Note  that  elastic  (PI),  does 
not  contribute  to  total  transfer. 

Removal 

t 

The  removal  cross  section  is  that  value  which  when  multiplied 
by  its  broad  group  flux  gives  the  rate  of  removal  of  neutrons  from 
that  group.  Thus  the  removal  cross  section  is  the  sum  of  absorption, 
inelastic,  elastic  (PO)  and  the  n  -  2n  removal  cross  sections. 

Transport 

The  transport  cross  section  is  defined  as 

Itr  =  St  -  uZ,  =  It  -  I,1  (126) 

Thus  the  transport  cross  section  is  the  total  cross  section  minus  the 
elastic  (PI)  cross  section.  The  transport  cross  section  is  calculated 
as  the  sum  of  inelastic,  the  n-2n,  absorption,  and  elastic  (PO) 
cross  section.  The  transport  cross  section  thus  includes  both  within 
group  and  out-of-group  scatter  terms. 

Summary 

"OLD  BARNYARD"  was  written  to  assist  professors  and  students  en¬ 
gaged  in  a  reactor  physics  coarse.  It  does  this  by  providing  the 
few-group  cross  sections  and  related  constants  usually  required  to 
solve  problems  assigned  in  such  a  course. 


The  code  calculates  tue  energy  dependent  flux  necessary  to  flux 
weigh  the  group  cross  sections  in  the  first  of  two  sequential  chains. 
The  epithermal  flux  is  calculated  by  the  neutron  transport,  moments 
method,  and  the  thermal  flux  is  calculated  by  solving  the  Wilkins 
equation  by  the  power  series  method.  Resonance  effects  are  included 
by  the  use  of  er>.,irical  equations.  The  second  chain  collapses  the 
eleven  fins  epithermal  groups  to  any  number  of  broad  groups  from  one 
to  eleven. 

Finally  one  should  always  keep  in  mind  the  limitations  of  these 
calculations.  The  accuracy  of  any  output  is  dependent  on  the  input 
cross  sections  employed.  It  must  be  remembered  that  the  input,  eleven 
group,  cross  sections  are  the  result  of  collapsing  the  68  group  cross 
sections  from  GAK-1  using  two  specific  standard  problems.  This  large¬ 
ly  limits  the  accuracy  of  the  code.  Also,  the  calculation  of  the  most 
probable  thermal  neutron  velocity  is  only  accurate  within  0.25  lethargy 
units.  At  the  present  time  resonance  calculations  are  effectively 
limited  to  one  nuclide.  This  is  because  chain  two  includes  resonance 
absorption  in  calculating  the  broad  group  cross  sections.  If  two 
resonance  nuclides  are  included  in  chain  two,  the  resonance  absorp¬ 
tion  effect  is,  in  effect,  counted  twice  by  assigning  all  the  resonance 
absorption  occurring  in  the  mixture  to  each  of  the  resonance  nuclides. 
Finally,  the  equations  developed  in  this  ’•  ter  and  used  in  the  code 
inherently  apply  limitations.  Recall  that  the  moments  equations  were 
based  on  an  infinite  homogeneous  medium  and  a  plane  source  of  infinite 
dimensions  at  the  origin.  Similarly  the  Wilkins  equation  is  applicable 
to  an  infinitely  heavy  Maxwellian  gas  and  a  source  free  medium. 
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Appendix  A 

Glossary  of  Symbols  Used  in  Chain  One 

The  following  is  a  list  o*  all  variables  used  in  chain  one  of  the 
code.  The  abbreviation  "ec"  is  defined  as  atom  fraction.  See 
Appendix  F  for  an  example  of  the  scattering  matrices.  In  the  scatter¬ 


ing  matrices  the  elements  of  the  matrix,  a.,  ,  are  defined  as  follows: 
a  diagonal  element  if  j=k,  an  upper  triangular  element  if  j<k,  a  lower 
triangular  element  if  j>k. 


Name 

Location 

Meaning  or  Use 

t 

l 

AN 

0804 

a  in  the  recursion  formula  of 
n 

eq  (100) 

1 

0805 

b  in  the  recursion  formula  of 
n 

eq  (102) 

• 

m 

0503+01 

a  ,  in  the  recursion  formulas 
n-i 

eq  (106) 

of  eq  (100) 

and 

. 

AN  2 

0503+05 

a  _  it.  the  recursion  formulas 
n-z 

above 

AN3 

0805 

a  ,  in  the  recursion  formulas 
n-o 

above 

AREA 

081G+01 

The  area  under  the  N(x)  vs  X  curve  (see  eq 

108) 

AREAX 

0809-03 

Summing  variable  to  find  the  area  under  the 

xN(x)  vs  X  curve 

ATOMS 

1061-01 

The  last  five  alphameric  characters  of  the 

th 

n 

• 

nuclides  identification 

- 

1025+03 

A  summing  variable  for  use  in  the  first  moment 

equation,  eq  (64) 


k:1  2s  +  En,n*  +  2En2n  i-k 

3=1 
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Ndi.-vi 

Location 

Meaning  or  Use 

C 

0810 

N'(x)  see  eq  (109) 

DELTA 

1061+02 

Absorption  parameter  of  the  n**1  nuclide 

0503 

Absorption  parameter  of  the  mix 

DENBR 

1061+01 

Total  nuclide  density  of  the  mix  (atoms  per 
bam-cm) 

DENS 

1061-01 

A  code  number  used  in  chain  2 

1051+04 

The  quantity  fl/3(  for  1136  in  e9  (65) 

DSUM 

0806+04 

Summing  variable  for  terms  of  the  derivative  of 
M(x)  from  the  Wilkins  equation 

DTERM 

0806 

Derivative  of  the  n**1  term  in  the  power  series 
expansion  of  the  Wilkins  equation 

E 

1002+01 

The  E  factor  for  heterogeneous  cell  calculations: 
E  =  1  for  homogeneous  mixture 

ELSQRD 

0851+02 

Thermal  diffusion  length  squared 

EMSQRD 

0852+01 

The  migration  area 

F 

1029+01 

A  normalizing  factor  such  that  F[FLX(10)]  =  1 

0808+01 

The  value  of  N(x)  for  the  point  x  =  p 

0816+01 

Flux  per  unit  lethargy  x2N(x) 

0810+04 

A  normalizing  factor  for  thermal  fluxes  such 
that  *(u=16)  =  1 

0109+02 

Storage  for  the  maximum  group  flux.  Used  to 
normalize  the  flux  spectrum  to  a  maximum  value 
of  100 

FF 

1002+01 

The  F  factor  for  heterogeneous  calculations: 

F  =  1  for  homogeneous  mixture 

FIS 

1061+02 

A  code  number:  Is  this  a  fission  nuclide: 

0  =  NO,  1  =  YES 

FLX 

1027 

The  group  flux  times  the  removal  cross  section. 
It  is  o  in  the  moments  equations 

1023 

The  group  flux  per  unit  lethargy 

Ill 


Naitia  t-oca+ion 


0500+08  Internal  index  used  to  read  data  for  all 

nuclides  in  the  mix 


0503+09 


0814+01 


Index  used  to  identify  the  energy  group  in  the 
thermal  calculations 

A  gate  for  locating  the  peak  of  the  N(x)  vs  x 
curve 


1006+03 


1006+01 


INDEX  1043+01 


An  initial  value  for  the  index  J  in  SS  109+01 
if  a  thermal  calculation  was  not  performed 

The  maximum  number  of  groups ,  below  any 
scattering  group,  to  which  scatter  occurs 

The  number  of  groups  from  which  scatter  occurs 

Code  indicating  what  source  has  just  been  read 


1065+01 


1065+03 


1007-02 


3001+01 


0500+01 


1067+02 


1067+01 


Internal  index  used  in  various  DO  loops 

A  code  number  for  the  type  of  resonance 
absorber:  l=UorW,  2  =  Th 

A  gate  to  determine  if  the  desired  source  has 
been  read  in:  0  -  still  seeking  source,  1  = 
found  source.  Read  rest  of  source  deck  into 
a  dead  end 

Converts  log  (FLUX)  to  an  integer  for  flux 
spectrum  printout 

The  group  number  of  the  lowest  group  to  which 
scatter  occurs,  from  a  given  group 

A  sorting  variable  in  outputting  the  fast 
fluxes 

A  gate  to  separate  thermal  and  epithermal 
fluxes  and  a  counter  to  sort  the  epithermal 
fluxes 

A  gate  to  initialize  TSIG9,  TSIG10,  and  PSIG9 

A  gate  to  index  KZ  by  increments  of  2  after  the 
first  elastic  (P0)  and  elastic  (PI)  are  stored 

A  counter  used  in  sorting  the  XSAH's 


1 

■ - — 

i 

■ 

Name 

Location 

Meaning  or  Use 

I 

LE 

2002+01 

Index  used  to  convert  the  k  index  of  0j-*k  of 
elastic  (P0)  and  (PI)  for  proper  storage  in 
the  lower  half  of  the  scattering  matrices 

I 

• 

LGP 

1001 

Subscripted  variable  to  read  the  problem  I.  D, 
card 

1 

• 

1047 

Used  to  print  out  the  group  with  the  unit 
source 

■ 

i 

1050 

Identifies  source  as  users  own 

1 

♦ 

i 

LQ 

1061+02 

A  code  number:  indicates  type  of  elastic 
scatter:  1  =  (n,n*),  2  =  (n,2n),  3  -  both, 
h  =  neither 

1 

> 

LS 

1061+02 

The  number  of  entries  in  the  cross  section 

library  starting  on  the  second  card  (except  ; 

vof  if  present) 

m 

i 

2002 

Index  used  to  convert  the  j  index  of  aA  .  of 
elastic  (P0)  and  (PI),  see  LE  above  J 

1 

i 

i 

( 

H 

1067+01 

A  code  number:  indicates  type  of  media: 

0  =  homogeneous,  1  =  heterogeneous 

1 

* 

! 

.1 

HOD 

1067+01 

A  code  number:  Is  nuclide:  1  =  moderator, 

2  =  fuel,  3  =  other 

■ 

i  ' 

f 

i 

NA 

1061+02 

A  code  number:  Is  epithermal  absorption 
appreciable:  0  =  NO,  1  =  YES 

I 

i 

< 

i 

NEW 

1061+01 

Number  between  2  and  21  indicating  source  to 
be  used 

i 

V 

1 

NRNUc 

1060+01 

The  number  of  nuclides  in  the  problem 

i 

i 

1 

1 

P 

1061-01 

Nuclide  density  of  the  n^  nuclide  (atoms 
per  barn-cm) 

£ 

1061+01 

Atom  Fraction  cf  the  n**1  nuclide 

■ 

) 

1 

i 

1051+07 

A  summing  variable  for  the  ^0^  in  eq  (65) 

! 

0503+06 

The  point  x=a  for  use  in  eq  (102)  for  a  series 
expansion  of  the  Wilkins  equation  about  the 
point  x=a 

I 

' 

0808 

The  value  of  the  variable  x  for  the  n**1 

thermal  group  j 

I 

f  ! 

i 

;  j 
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Name 

Location 

0813 

PREVF 

0503+02 

0816 

PREFX 

0809-02 

PROB 

0728 

PSI 

1061+02 

PSIGF 

500+02 

PSIGM 

0500+03 

PSIGO 

0500+04 

PSIG9 

0705+04 

PU 

1045 

0995 

RES 

1060+01 

RHO 

1051+01 

RIEFF 

0995 

SA  1027+01 


Meaning  or  Use 

Lower  lethargy  boundary  for  the  various  thermal 
group  fluxes 

A  storage  variable  for  N(x)  from  the  previous 
calculation 

N(x)  at  the  peak  of  the  h(x)  vs  x  curve 

Temporary  storage  of  xN(x)  from  the  previous 
calculation.  Used  to  find  the  area  under  the 
xN(x)  vs  x  curve 

The  resonance  escape  probability 

The  average  logarithmic  energy  decrement  £ 

Sum  of  the  products  of  the  atom  fraction  times 
the  microscopic  thermal  absorption  cross  sec¬ 
tion  of  fuel  nuclides  only 

Same  as  PSIGF  except  for  moderator  nuclides 
only 

Same  as  PSIGF  except  for  other  nuclides  only 

Sura  of  the  products  of  £  times  average  of  the 
scattering  cross  section  of  groups  9  and  10  for 
the  mixture  (Average  scattering  cross  section 
in  the  resonance  region) 

A  dead  end  to  read  through  the  remaining 
undesired  sources  after  the  desired  one  has 
been  found 

Nuclide  density  of  the  resonance  nuclide(s)  in 
atoms  per  barn-cm.  Immediately  converted  to 
atom  fraction 

A  code  number:  Is  e  resonance  calculation  to  be 
performed:  0  =  NO,  1  =  YES 

A  holding  variable  for  the  lower  lethargy 
(upper  E)  boundary  of  the  n**1  epithermal  group 

A  code  number:  Is  the  resonance  integral  to  be 
calculated:  0  =  calculated  by  code,  non-zero  = 
value  of  resonance  integral  to  be  used. 

First  moment  of  tne  fourier  transformed  flux 


Name 

Location 

Meaning  or  Use 

SAM 

1061-01 

The  first  five  alphameric  characters  of  the 
nuclide's  identification 

1025+02 

The  summation  term  in  eq  (64).  Also,  with 
r j ,  it  is  used  in  eq  (66) 

SAMM 

1061-01 

Tht  second  five  alphameric  characters  of  the 
nuclide's  I.  D. 

0730 

Recriprocal  of  3  times  the  total  nuclide 
density 

1023+01 

The  age  to  the  Indium  resonance  (1.46ev) 

SCAT 

0601+01 

The  scattering  per  resonance  atom 

SIGA 

1061+02 

The  for  the  n**1  nuclide 

a 

SIGATH 

0851+01 

The  ox  for  the  mixture 
a 

SIGDEL 

0500+06 

The  sum  of  (o^A)  a.f.  of  the  mix  where  a.f. 
is  atom  fraction 

SIGS 

1061+02 

The  c^of  the  n^  nuclide 
s 

SIGSTH 

1066+04 

The  a  of  the  mixture 

8 

SIGT 

1066+03 

The  o*h  of  the  nth  nuclide 

SIGTT 

0991+01 

The  sum  of  the  product  of  the  atom  fraction 
tiroes  the  o*h  of  all  nuclides  except  fuels 

SIC9 

0705 

The  total  elastic  scatter  from  group  9  for  the 
nth  nuclide.  See  Appendix  F 

SIG10 

0705+02 

The  same  as  SIG9  except  the  scatter  is  from 
group  10 

SOURCE 

1043+01 

The  fraction  of  the  source  in  each  of  the  11 
epithermal  groups 

SR 

1005 

The  removal  cross  section 

SSN 

1009+01 

Scattering  matrix  for  (n,n*)  and  elastic  (P0) 
cross  sections,  (a  12  x  13  matrix):  Upper 
triangle  (J,  K+l)  for  (n,n*),  diagonal  and 
lower  triang’e  (LS,LE)  for  elastic  (P0) 

STR 

1004+02 

The  multigroup  transport  cross  section, 

VwEsj-j 

115 


Name 


Location 


Moaning  or  Use 


SUM  0806+03  Summing  variable  for  terms  in  the  power  series 

expansion  of  the  Wilkins  equation.  See  eq 
(100)  and  (102). 

S2N  1009+01  Same  as  SSN  except  it  is  for  (n,2n)  and  elastic 

(PI)  values 

T  1027+02  The  group  age  i  .  use  in  eq  (66) 

TEMPTR  1060+01  The  system  temperature  T,  in  degrees  Kelvin 

TERM  806+02  The  individual  terms,  a  Xn,  in  the  solution  to 

the  Wilkins  equation,  eq  (100)  and  (102) 

TERMU  1061+02  (1-{S)  thermal  of  the  nth  nuclide 

TOTWT  1027+04  A  summing  variable  for  the  summation  term  in 

eq  (66) 

TSIG9  705+01  The  total  scattering  cross  section  from  group 

9  for  the  mix 

TSIG10  705+03  Same  as  TSIG9  except  it  applies  to  group  10 

U  1065+02  Subscripted  variable  for  the  upper  lethargy 

(lower  E)  boundaries  of  the  11  epithermal 
groups 

UTIL  502  The  thermal  utilization 

VF  1061+02  The  vo^of  the  nth  nuclide 

VPR0I)  0815  The  most  probable  thermal  neutron  velocity 

VBAR  854+01  The  average  thermal  neutron  velocity 

't’H 

X(K)  1027+03  The  lethargy  width  of  the  n  group 

603+01  Storage  for  the  symbols  used  in  the  flux 
spectrum  printout 

XBAR  810+02  The  average  normalized  thermal  neutron  velocity 

x 

XSAM  0997  Temporary  storage  for  the  cross  section  library 

1030  Storage  for  the  fast  and  thermal  fluxes  prior 

to  their  output 

0698  Temporary  storage  for  the  grid  for  spectrum 

printout 
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Name 


XTON 


Y 


Location  Meaning  or  Use 

0120  This  prevents  computer  underflow  if  the  flux 

is  too  low.  Note  log  (0)  =  -« 

0806  xn“l  in  equation  (100)  or  (102) 

0806+01  x°  in  equation  (100)  or  (102) 

0503+08  The  value  of  x  corresponding  to  lethargy  U  -  23 

0809+02  The  point  x  for  the  next  Wilkins  calculation: 

x2  =  X1  +  X1  CexP  (*125-1)3 
The  quantity  (x-a)  in  equation  no.  (102) 


0810 


n  i*  n  n 
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Appendix  B 


Chain  One  Source  Beck 


THIS  CARD  ORIGINS  THF  PROGRAM  AT  ADDRESS  10970  IN  MEMORY.  IN  ORDER 
TO  DO  THIS  YOU  MUST  USE  A  SPFCIAL  SET  OF  SHORT  SU9ROUT l NES. 

JOS  OLD  BARNYARD.  CHAIN  1 
ORIGIN  10970 

STORING  IN  COMMON  THOSE  DATA  USED  IN  CHAIN  2  OR  IN  SUBPROGRAMS 
COMMON  LGP ( 13),KE,FLX(11>,U!11) . SOURCE (  1 1 ) .DENBR  .NRNUC  .  TEMPTR , 
lXBAR.Pf IGM.R1EFF.RES 

DIMENSION  SA( 11 ) *STR( 11 ) ,SR( 11) .XSAM1 168) »SSN( 12.13) »S2N( 1 ?  *  1 3 ) »X( 
111)  »T(12) 

C  READING  IN  THE  FIRST  TWO  DATA  CARDS 

1001  READ  1060.LGP 
PUNCH1C60.LGP 
1060  FORMAT ( 1 3A5 ) 

READ  1067 .NEW, NRNUC. DENBR .TEMPTR.  RE S.M 
PUNCH1065 
1065  FORMAT!//) 

TYPE  1060 .LGP 

C  READING  IN  THE  DECK  OF  SOURCES 

REAC  1070. U 
K  =  0 

DO  1046  J=1.10 
IF(K)iC45. 1043. 1045 

1043  RFAD1060.LGP 
READ1O70, INDEX, <  3URCF 

1070  FORMAT! 12N) 

IF! INDEX-NEW) 1046.1044,1046 

1044  K= 1 

GO  TO  1046 

1045  READ  K-60.PU 
READ.PU 

1046  CONTINUE 

IF INEW-1 0)1 050,1049, 1047 

C  SETTING  THE  UNIT  SOURCE  IN  ITS  GROUP  (IF  APPLICABLE) 

1047  SOURCE (NEW-10)=1. 

LGP! 10)=NEW-10 
PUNCH  1071.LGP 

1071  FORMAT (9A5.I3.3A5) 

GO  TO  1048 

C  READING  IN  THE  USERS  OWN  SOURCE  (IF  APPLICABLE) 

1049  READ1070. SOURCE 

1050  PUNCH  1060. LGP 

1048  TYPE  6C0 

600  FORMAT! /16HBEGIN  EPITHERMAL) 
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PUNCH  1058.TEMPTR 

1088  FORMA  T  (  /22HSYSTFM  T  EMC  ERA  TURF.  IS  ,F6.1»7ri  KtLV'.N) 

PUNCH  1059 

1059  FORMAT ( / 1 3HA TOM  FRACTIONS  ARF /  ) 

C  CLEARING  THF  MATRICES  OF  ANY  PREVIOUS  PROBLEM 

DO  1000  J= 1  *  1 1 
STR ( J ) =0 • 

1000  SR ( J i =0 • 

DO  1002  J= 1  *  1 2 
DO  1002  <=1*13 
SSN( J,K>=0. 

1002  S2N ! J  »K ) =0 • 

FF=  1  . 

E=  1  • 

C  READING  IN  THE  RESONANCE  DATA  CARD  (IF  APPLICABLE) 

IF(RES)995,994,995 
995  READ, JJJ,PU,RIEFF 
PU=PU/DENBR 

(  IF  APPLICABLE ) 

READ  THE  F  AND  E  FACTORS  FOR  HE  T EROGENOUS  CELL  CALCULATIONS 
994  IF(MJ 500*500.501 
501  RFAD.FF.E 

SETTING  UP  THE  DO  LOOP  FOR  EACH  OF  THE  NUCLIDES 
500  PR08=1. 

KKK=0 

PSIGF=Oe 

PSIGM=0. 

PSIGO=0. 

SIGTTs.O 
SIGDEL=.0 
SIGSTH1  0» 

DO  1020  1=1 »NRNUC 

READING  IN  THE  NUClIDE  DATA  CARDS 
READ  1061,  SAM,  SAMM,  ATOMS,  P,  MOD,  DENS 
1061  FORMAT (5X»3A5,10N1 
P=P/DENBR 

READ  1  "'61  ,  SAM,  SAMM,  ATOMS  ,LS»NA,LQ,PSI  ,SI3A,  DELTA, SI  GS  ,  TERMU,  VF  ,F  I S 
PUNC H 1066, SAM, SAMM, ATOMS, P 
1066  FORMAT (5X,3A5,F15«8) 

CORRECTING  THE  THERMAL  CROSS  SECTIONS  FOU  TEMPERATURE  AND 
CALCULATING  THE  THERMAL  SIGMA  TOTAL  FOR  EACH  NUCLIDE 
S I GA=SIGA*2 • 2E+( 5/( 1 • 28E+04*SQRT ( TEMPTR ) ) 

DELTA=DELTA*2*2E+05/( 1.28E+04*S0RT( TEMPTR) ) 

SIGT  =  S  I  GA+S‘I  GS 

CALCULATING  TCTAL  THERMAL  SCATTERING  SIGMA 
5IGSTH-SI GSTH+P*S I GS 

C  CALCULATING  THERMAL  UTILIZATION  AND  ABSORPTION  PARAMETER 

GO  TO  (993,992,996! , MOD 
992  PSIGF=PSIGF+P*SIG4 
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GO  TO 

997 

996 

PS  I GO= 

f  SIGC+P+S 

IGA 

GO  TO 

991 

993 

PS  I  GM= 

PSIuM+P*S 

IGA 

991 

SIGDEL 

=  sigdel*-s 

I  GT*P*DEL I  A 

5  IGTT  =  5  I GT  T  +P*  S I 6T 

C  READING  IN  ThE  WHOLE  CROSS  SECTION  LIBRARY  FOR  THE  NUCLIDE 

997  READ  1067  * ( X  SAM (J) »  J= 1 »LS ) 

1067  FORMAT (6N) 

KZ  =  0 
KKZ  =  0 

C  CHECKING  FOR  ABSORPTION 


IF(NA) 1006. 

1006*1004 

1004 

DOlf 05  J=1 ♦ 

11 

KZ=K'Z  +  1 

STR( J)=STR( 

J)+XSAM(K7)*P 

1005 

SR ( J ) -  SR ( J  ) 

+  XSAM ( KZ ) VP 

1006 

KZ=KZ+1 

IN=XSAK(KZ ) 

KZ=KZ+1 

I  IN=XSAM(KZ 

) 

DO  1017  J=1 

,  IN 

KK=J+I IN 

IF(KK-12 ) 1008*1008, 1007 

1007 

KK=12 

1008 

DO  10)7  K=J 

♦  KK 

GO  TO  (1009 

.1010.1009,1011  )  *L0 

1009 

KZ=KZ+1 

C  INELASTIC  CROSS  SECTION*  GROUP  TO  GROUP 

SSN( J.K  +  1 )=SSN( J.K+1 )+XSAM(KZ  )*P 
GO  TO  1012 

1010  KZ=KZ+1 

C  N-2N  CROSS  SECTIONS*  GROUP  TO  GROUP 
S?N( J.K+1 )=S2N( J.K+1 )+XSAM(KZ)*P 
GO  TO  1012 

1011  IF(KKZ)2001. 2000*2001 

2000  KZ=K7+3 
KKZ  =  1 
GOT02002 

2001  KZ=KZ+2 

C  PUTTING  ELASTIC  ( PO )  AND  ELASTIC  (PI)  IN  THE  BOTTOM  OF  SSN  AND  S2N 

2002  LS=13-„ 

LE=13-K 

SSN(LS.LE) =SSN(LS*LE;+XSAM(KZ )*P 
S2N(LS*LE) =S2N(LS,LE)+XSAM(KZ+1 )*P/3. 


1012 

STR ( J ) 

=  STR(  J 

)+XSAM( KZ )*P 

IFIK-J 

) 1016, 

1013.1016 

1013 

GO  TO 

( 1017, 

1014,1017,1015)  , 

1014 

SR ( J ) = 

SRI J>- 

XSAM(KZ  )*P 

GO  'O 

1017 

101S 

STR ( J ) 

=  STR  ( J 

)-XSAM( K7+1 )*P/3 

GO  TO 

1017 

120 


1  °  1  6  SKI  JJ=5,R(^)+XSAM(KZ  )*P 

mi  7  CONTINUE 

GO  TOf 1018.1018.1019,700) .LG 

1018  10=4 

GO  to  moft 

1019  LQ=. 

GO  TO  1006 

C  READING  THE  FISSION  CROSS  SECTIONS  USED  IN  CHAIN  2  INTO  DUMMY  STORAG 

700  IF(FIS)t98,698.69V 

699  READ  1067 , ( X  SAM ( J ) , J  =  4 1 » 5 1 ) 

C  READING  IN  THE  GRID  BACKGROUND  FOR  THE  SPECTRUM  PRINTOUT 
698  READ  2722, (XSAM(J) ,J=41.120> 

2222  FORMAT ( £041 > 

C  CHECKING  FOR  RESONANCE  CALCULATION 

If (KKKJ705.701 . 705 
7^1  IF(RES)704. 1020, 704 

704  KKK= 1 

C  CALCULATING  PSI*SIGMA-S  FOR  THE  MIXTURE 

T  S I G9=0 . 

TSIG10=0. 

°SIG9=0« 

705  S I G9=SSN (4.4)+SSN(4,3) +SSN (4.2) +SSN (4»l)-TSlG9 
TSIG9=SSN(4,4)+SSN(4,31 +SSN (4.2) +SSN (4.1) 

SIGl0=SSN(3,3)+SSN( 3,2>+SSN(3,l )-TSIGl0 

TS I C  10=SSN  (  3  » 3  )+SSN  (  3.2  )  4-SSN (3.1) 

PSI G9=PS I G9+PS I *P* ( SI G9+S I G1 0 ) *• 5 

1020  CONTINUE 

C  ALL  CROSS  SECT IC  NS  HAVE  BEEN  CALCULATED  AT  THIS  POINT 

IF(KKK)720.730,720 

720  TYPE  601 

601  FORMAT ( / 1 5HBEG IN  RESONANCE) 

C  CALCULATING  SCATTERING  PER  RESONANCE  ATOM  (IF  APPLICABLE) 

SC  AT=«5*(TSI  G9+T  S  I  G1>j  )  /PU 

C  SELECTING  THE  EMPIRICAL  EQUATION  FOR  RESONANCE  INTEGRAL 

IF  (RIEFFJ728.721.728 

721  GO  T0( 722,725' ,JJJ 

722  IF(SCAT-4000. >723,723,724 

727  RIFFP=2,69*SCAT**,471 

GO  TO  778 

724  RIEFF=EXPF(5.64-163./(SCAT**.65>  ) 

GO  TO  728 

725  IF(SCAT-4500. >726,726,727 

726  RIEFF=8.33*SCAT**,253 
GO  TO  728 

727  R I EFF=7  0. 

C  CALCULATING  RFSONANCE  ESCAPE  PROBABILITY 

778  PROB=FXPF(-PU*R!EFF/OSIG9> 
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C  ADJUSTING  GROUP  9  AND  10  REMOVAL  CROSS  SECTION 

c  for  resonance  escape 

SR(9)=SR(9)  fRIFFF*.  1 *DU/L0GF ' 6 1 . 442 / 1 .  1254) 

SRU0)=SRt  101+R  '7,7F*.9*PU/L0GF(61.A42/1.  1254) 

STRt  P)  =SR( 9)  •  SS,  ,4)-S2N«4,4) 

STR(10)=STR(  1  0  )  -i  jSN  (  3.31-S2N!  3,31 

C  HANSENS  FINITE  DIFFERENCE  EQUATIONS  FOR  EPITHERMAL  FLUX  AND  AGE 

730  SAMM=1./0.*DENBR*DEN8R> 

PUNCH  1051 

1051  FORMAT ( //57X « 14HRELAT I VE  GROUP/5HGRGUP * ?X , SHLEThARGY  »  7X  » 10HENFRGY , 

1  EV»8X  »6HS0URCE  *6X» 14HFLUX*  /UNIT  U/ 1 
RHC=.0 

DO  1029  K=  1  ♦  1 1 
LEa13-‘' 

DENS=SAMM/<STR(K)*SR(K)  1 
AT0MS=0. 

T0TWT=0. 

P=0. 

T (< 1 =0» 

FLX(K1=0. 

SA ( K 1 =0* 

IF(K-1)1027, 1027, 1025 

1025  DO  1026  J-1,<-1 
l.S*  13— J 

SAM=FL> ( J)*(SSN(LS*LE)+SSN( J ,K+1 ! +2 . *S2N ( J,K+1) )/SR( J) 

ATOMS=ATOMS*SAM 

TOTWT  =  TOTWT  +SAM*T ( J 1 

1026  P=P*S2N(LS«LE)*SA(J> 

1027  FLX ( K ) =  ATOMS+SOURCE ( K 1 
SA(K)=FLX{K)*DENS+P/STR<K) 

T(K)=(SA(K)+TOTWT ) /FLX(K) 

X ( K ) =U( K ) -RHO 

RHO=U(K) 

1029  CONTINUE 

F=SR( 101*4. /FLX( 10) 

C  CALCULATING  TOTAL  FLUX  PER  GROUP 

DO  1023  K= 1  ♦  1 1 

1023  FLX(<1=FLX(K)/SR(<3 
DO  1 03C  K=1 » 10 

'  030  X3AM  (  40-K  )  -=F*FLX  (  K  1  /X  {  K ) 

C  NOTICE  THAT  TOTAL  FLUX  PER  GROUP  HAS  BEEN  LEFT  IN  STORauc  IN  rLX(ll) 

C  CALCULATING  AGE  TO  INDIUM  RESONANCE 

SAMM= ( T ( 10 1 -T ( 9 1 }*0.935+T(9) 

C  CHECKING  FOR  MODERATOR  IN  THE  SYSTEM 

IF(PSIGM) 1024, 108 » 1024 
C  .  STARTING  THE  THERMAL  CALCULATION 

1024  TYPE  602 

602  FORMAT! '13HSEGIN  THERMAL) 

C  THE  ABSORPTION  PARAMETER  AND  Il/l-F) 

UTIL=0, 
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IF ( PS  I GF )502»5P3,502 

S^?  UTIL  =!./((  (P5IGM+PSIG0)  /PS  IGF  !*FF+E) 

SO*  OFL  T4  =  5  IGDEL/SIOTT* ( 1 . /( l.-UT IL ) ) 

ARtAX=»0 
PREFX=.0 
PREVF  =  0  • 

11=0 

AN2=1. 

P  =  0. 

C  ARRANGING  FOR  INCREMENTS  OF  0.25  LETHARGY  UNITS 

ETODU=FXP< .1251-1. 

C  Y  HERE  EQUALS  X.  SINCE  X  WAS  SUBSCRIPTED  ABOVE 

Y  =  EXP( - 1 1 • 5  »  *SORT ( 1.0E+12/ ( 8.616*TEMPTR ) ) 

1=0 

C  DOING  THE  SERIEI  SOLUTION  TO  WILKINS  EQUATION 

AN1=AN?*DELTA*. 333333 
DSUM=AN1 

80?  SUM=AN2+AN1*Y 
XTON=Y 
AN3=0. 

N=4 

803  N=N+ I 
XN=N-3 

IF  I P) 805.804.805 

C  THE  RECURSION  RELATION  FOR  X=0 

804  AN=(2.*(XN-2. / *AN2+DEL7A*AN1 ) / ( XN* I XN+2. ) I 
GO  TO  806 

c  THE  GENERAL  RECURSION  RELATION 

805  AN= I AN1* ( XN-1 • ) * ( XN+1 »-2»*P*P ) — AN2* ( 4.  *P* ( XN-2 • ) +DELTA ) -AN3*2»* ( X 
IN-3. 1 > / ( P*XN* ( l.-XN) ) 

AN3= AN2 

806  DTERM=XN*AN*XTON 
XTON=XTON*Y 
TERM=AN*XTON 
:.UM=SUK+TERM 
DSUM=DTERM+DSUM 
AN?=AN1 

AN1=AN 

C  CHECKING  FOR  CONVERGENCE  OF  THE  SERIES 
IF< A8S < TERM/ SUNU-1.0E-06) 307, 807, 803 
8Q7  IF (A9S(DTERK/DSUM)-1.0E-06) 808,808,803 
808  P= Y+P 

C  FINDING  N  t  X 1  FOR  THE  POINT  X  =  P 

F=EXP( -P*P ) *SUM*P**2 

C  LOCATING  THE  PEAK  OF  THE  NIX)  CURVE 

I F ( I  I  1816,814,816 
814  IFIF-PREVF 1815,815,816 

C  CALCULATING  THE  40ST  PROBABLE  VELOCITY 
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n  n 


8 1  s  VPROB  =  SOR  T  (  1 .6477C  +  0A*  YEMPTR  )  *P 
t  1  =  1 

816  PRCvPsF 

C  CALCULATING  FLUX  PER  UNIT  !J  WHICH  IS  X  *  X  x  N  (  X  ) 

F  =  F‘«P4  p 
1  =  1+1 

C  STORING  THE  VALUES  OF  THERMAL  FLUX 

X c- AM  (  n=F 

C  CALCULATING  THE  AREA  UNDER  THE  XN(X)  CURVE  BY  TRAPEZOIDAL  RULE 

F  =  F/P 

ARFAX=AREAX+( F+PREFX ) * Y*. 5 
PREF  X  =  F 

!F( 1-29) 800*810*810 

809  AN2  =  SUF. 

AN1=DSUM 

y  =  P  +  P*E  T ODU 
Y  =  Y-P 
GO  TO  802 

C  CALCULATING  THE  AKEA  UNDER  THE  FLUX  CURVE  ANALYTICALLY 

810  C=2.*F*( l./(P*P)-l. )+P*P*EXP( -P*P) *D$UM 
ARE/ =(P*C+?.*(P*P-1.1*F/P) /DELTA 

C  FINDING  the  AVERAGE  NORMALIZED  VELOCITY 

X8AR=AREAX/AREA 

NORMALIZING  THE  rLUXES  AND  MATCHING  THE  THERMAL  AND  EPITHERMAL 
FLUXES  IF  A  THERMAL  CALCULATION  WAS  PERFORMED 
F=1./XSAM(29 . 

DO  812  J= 1 « 29 
812  XSAMt J)=F*XSAM< J) 

GO  TO  109 

108  11=30 
TYPE  1028 

1028  FORMAT (36HTHERE  IS  NO  MODERATOR  IN  THE  SYSTEM. /36HN0  THERMAL  CALCU 
1LATI0N  WILL  BE  DONE. ) 

109  F=0. 

C  SETTING  THE  MAXIMUM  FLUX  TO  100 

DO  105  J= I  I  *  39 
I F ( F-XSAM ( J ) ) 106 » 106 1 105 
106  F  =  XSAMJ) 

105  CONTINUE 

DO  111  J= 1 1  *  39 

111  XSAM ( J ) =  XSAM IJ)*100»/C 

C  OUTPUTTING  THE  FLUXCS 

0=16.5 

DO  131  <=1.23 
IFIK-10) 112.112.113 

112  <<=A0— K 

1022  PUNCH  1052.<»U(<) » 1.0E+07*EXP (-U(K) ) .SOURCE l O .XSAM(KK) 

GO  TO  131 

113  IF(DSIGM) 115,131.115 
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115  P=P+ • 5 

j=47-<-K 

313  PUNCH  1052.K.P.I .0E+07*EXP(-P) *0. »XSAM{ J> 

131  CONTINUE 

1052  FORMAT ( I3.F15.2.3X.1P3E16.4) 

C  OUTPUTTING  THE  THERMAL  CONSTANTS 

0|JNCH  (-5  0  *  SAMM 
IF(PSIGM) 130*  860  *  130 
130  PUNCH  851 »T ( 10  I 

850  FORMAT! //5X.35HAGE  TO  INDIUM  RESONANCE  (1.46EV)  IS* 1PE20.4,4H  CM2) 

851  FORMAT ( /5X  *  36HAGE  TO  ARBITRARY  THERMAL  !1.12EV)  IS*  1PE19.4*4H  CM2 
1  ) 

C  CALCULATING  THE  DIFFUSION  LENGTH  SQUARED  AND  MIGRATION  AREA 

SIGA1 K=PSIGM+PSIGF+PSIGO 

ELSQRD=.33333«SIGSTH*XBAR/(SIGATH*!DENBR*'SIGSTH+SIGATH/XBAR) ; **2 ) 
PUNCH  852  *ELSQRD 

852  FORMAT  (  /  5X  *  3  5HTHEKMAL  DIFFUSION  LENGTH  SQUARED  IS*lPE20.4»4l-i  CM2) 
EMSQRD=T ( 10 ) +ELSQRD 

PUNCH  853 » EMSQRD 

853  FORMAT! /5X*23HTOTAL  MIGRATION  AREA  IS* 1PE32.4.4H  CM2) 

PUNCH  (  54  » VPROP 

854  FORMAT! /5X,41HM05T  PROBABLE  THERMAL  NEUTRON  VELOCITY  IS* 1PE14.4.6H 
l  M/SEC) 

VBAR  =SQRT < 1.6477E+04*TEMPTR )*XBAR 
PUNCH  £5 5. VBAR 

855  FORMAT (/5X.35HAVERAGE  THERMAL  NEUTRON  VELOCITY  I S» 1PE20*4*6H  M/SEC 
1) 

PUNCH  856. DELTA 

856  FORMAT! /5X.23HABSORPTION  PARAMETER  IS»1PE32.4) 

C  OUTPUTTING  THE  RESULTS  OF  THE  RESONANCE  CALCULATION 

IF;RES1859, 860*859 
859  PUNCH  857, SCAT 

857  FORMAT! /5X»  32HSCATTERING  PER  RESONANCE  ATOM  IS* 1PE23.4.6H  BARNS) 
PUNCH  858.RIEFF 

858  FORMAT! /5X *31HEFFECT I VE  RESONANCE  INTEGRAL  I S » 1PE24.4 ,6H  BARNS) 

86n  PUNCH  861 

861  FORMAT ! 78X  » 1H- ) 

C  WARNING  THE  USEF  TO  SET  SWITCH  FOR  SPECTRUM  OUTPUT 

TYPE  862 

862  FORMAT! /52HPUT  SWITCH  1  ON  FOR  SPECTRUM  PRINTOUT.  PRESS  START.) 
PAUSE 

IF! SENSE  SWITCH  1)3000.2999 

C  THE  SPECTRUM  PUNCHOUT  SUBPROGRAM 

29Q9  TYPE  863 

863  FORMAT! /28HNO  SPECTRUM  WILL  BE  PUNCHED.) 

GO  TO  1091 

300C  TYPE  603 

6N3  FORMAT! /14HBEGIN  SPECTRUM) 

C  STORING  THE  SYMBOLS  FOR  THE  GRAPH 

X { 1 ) = 1 HX 

x(2)=:hl 
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X ( 3 ) = 1 HH 
PUNCH  3l'l 

300  FORMAT ( 69HPlOT  OF  SPECTRUM,  IN  FLUX  PER  UNIT  LuTHARGY  VERSUS  LETHA 
1RGY,  FOLLOWS,//! 

PUNCH  3P 1 

ADJUSTING  MINIMUM  FLUX  TO  PREVENT  COMPUTER  UNDERFLOW  AND 
REDUCING  FLUX-S  TO  LOG  SCALE  80  UNITS  WIDE 
DO  300!  J= 1 1  *  39 

IFIXSAMI J) -A. 902E-03) 120.3001  ,3601 
120  X  SAV  (_i5=A«°82E-03 

jO'M  X  S  AM { J )  =  ( t LOG! XS AVI  J) ) /2. 302891+1. 1^20.  +  .5 
KK  =  9 
1=0 

C  ONE  PASS  THROUGH  THE  LOOP  FOR  EACH  LINE  OF  OUTPUT 

DO  3005  J=II ,93 
IF(KK)3011, 3014, 3011 

C  CHECK  I  MG  FOR  EPITHERMAL  OR'  THERMAL  FLUXES 

3011  IF! J-29) 3009,3009,3012 

3012  PU=J 

P  =23. 2-PU/A. 

C  ADJUSTING  FOR  DIFFERENT  WIDTH  GROUPS 

IF(P-U(KK)  )  301 3  »  30 1 3 , 3014 
301 ^  KK=KK-1 

3014  K=XSAN (39-KK) 

GO  T0  °oo8 

3000  K=XSAM(J) 

3008  1=1+1 

C  NUMBERS  FOR  THE  LETHARGY  SCALE 

LE=23-J/a 

C  CHECKING  FOR  VALUES  TOO  LARGE  FOR  SCALE 

IF( <-79)3015*3016  *  3006 
3006  N=  3 
<  =  79 

GO  TO  3004 

C  CHECKING  FOR  VALUES  TOO  SMALL  FOR  SCALE 

3015  IF(K-3)3C19, 3017, 3016 

3"!7  GO  TO  '3018,3016, 3016, 3016, 3018), I 

3018  N=2 
K  =  4 

GO  TO  3O0A 

3019  1F(K-2)3018, 3017, 3017 
yn 16  N=1 

3004  KZ  =  <  +  3S 
KKZ=K+4 l 

GO  T0( 3002,3003,3003,3003,3002 ) , I 

C  OUTPUTTING  THE  I  >ECTRUM  GRAPH 

3002  PUNCH  302, LE , ( XSAMI LS) »LS=43*KZ ) »X(N) , ( XSAM(LS) ,LS=KKZ,120 ) 

302  FORMAT ( I  2 ,78A1  ) 


GO  TO  3005 

3003  PUNCH  2222  *  i  XSAMi  LS )  »LS=4i  «KZ  j  »X  i  N  i  »  i  XSAM  i  l5  j  »l5=KK2  ♦  120 ) 

3005  CONTINUE 
PUNCH  30 1 

301  FORMAT ( /2H.1 , 1 7X ♦ 1H 1 » 1 8X , 2H 10 ♦ 1 8X . 3H100 » 1 5X »4H1000/ ) 

1091  TYPE  1092 

1092  FORMAT ( /14HEN0  OF  CHAIN  1//39KTO  COLLAPSE  CROSS  SECTIONS  LOAD  CHA! 
IN  2/) 

C  THIS  PUNCH  MAKES  THE  407  LISTER  ADVANCE  A  SHEET  OF  PAPER 

PUNCH  861 

C  THIS  PUNCH  AVOIDS  RUNNING  OUT  THE  LAST  PUNCHED  OUTPUT  CARD 

PUNCH  1065 
STOP 
END 

S  EOJ 
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Appendix  C 


Glossary  of  Symbols  Used  in  Chain  Two 

The  following  is  a  list  of  all  variable  names  used  in  chain  two 
of  the  program.  Their  location  in  the  code  is  referenced  by  statement 
numbers,  and  their  meaning  or  use  in  the  code  at  that  location  is 
given.  Some  variable  terms  are  carried  over  to  chain  two  from  chain 
one  through  the  use  of  the  COMMON  statement  at  the  beginning  of  the 
chains.  -  See  Appendix  F  for  a  sample  of  the  scattering  matrices. 


Name  Location 

A  0054-01 


1109+02 


ATOMS 

0701 

DENS 

0701 

0061 

DELTA 

0702 

FIS 

0702 

FLX 

0003 

Meaning  or  Use 

A  summing  variable  for  all  the  group  to  group 
elastic  (P0)  flux  weighted  cross  sections 

p* 

(La.  ,  for  all  j  in  the  "scatter  from"  broad 
i 

group,  and  all  k  in  the  "scatter  to"  broad 
group) 

A  summing  variable  for  all  the  flux  weighted 
fine  group  absorption  cross  sections  in  any 
given  broad  group 

The  last  five  alphameric  characters  of  the 
nuclides  name 

A  number  indicating  whether  microscopic  or 
macroscopic  cross  sections  are  desired: 

1  =  micro,  2  =  macro 

If  macroscopic  cross  sections  are  desired, 
this  variable  is  set  equal  to  the  nuclide 
density  (atoms  per  barn-cm) 

Thermal  absorption  parameter  of  the  n**1  nuclide. 
Not  used  in  chain  two  calculations 

A  code  number:  Are  epithermal  vo^  included  for 
this  nuclide:  0  =  NO,  1  =  YES 

The  fraction  of  the  total  broad  group  flux  that 
exists  in  each  fine  group  within  that  broad 
group 
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Name 


Location 


I 


I  IN 
IRES 

J 


K 


0310 

0307 

1108 

0205+01 

1006+01 

1115 

0601+01 

0002-03 

1112 

0049 

1004 

1006+04 

1120 

4040 

0604+03 


Meaning  or  Use 

An  integer  number  used  to  decrease  the  request¬ 
ed  number  of  broad  groups  if  the  flux  in  a 
broad  group  is  too  low  (done  to  prevent  divi¬ 
sion  by  near  zero  quanities) 

Index  of  the  DO  loop  of  SS  301-02  to  SS  305 . 

Set  equal  to  12  to  immediately  the  DO  loop  if 
NBG  is  changed 

Index  to  determine  if  thermal  broad  group  will 
be  output 

Index  of  outermost  of  four  DO  loops.  It  indexes 
through  the  epithermal  broad  groups.  For 
scatter  from  group  n  to  group  m  it  is  the  value 
of  n 

The  number  of  groups  from  which  scatter  occurs 

A  code  number:  Is  this  a  resomnce  nuclide: 

0  =  NO,  1  =  YES 

A  subscript  to  determine  the  boundaries  of  the 
broad  groups 

An  index  used  in  normalizing  the  fraction  of 
the  broad  group  flux  that  exists  in  each  fine 
group 

Index  for  clearing  matrices 

Index  to  put  the  cross  section  library  into 
XSAM 

Index  to  sort  the  aa  into  their  proper  SA 

Index  of  "he  outer  DO  loop  used  to  load  the 
cross  sections  into  the  scattering  matrices. 

It  is  the  j  index  of 

Index  of  the  third  of  four  nested  DO  loops. 

For  scatter  from  broad  group  n  to  broad  group 
m,  it  accounts  for  all  fine  groups  in  r. 

Index  of  the  outer  DO  loop  for  broad  group 
output 

Temporary  storage  for  the  requested  number  of 
broad  groups 
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Warn* 


KK 


KKZ 


KZ 

L 


LE 


LGP 

LQ 


Loc?.ti™ 

0002-01 

1008 

1122 

1109 

1070+04 

1007-02 

1120+01 

1067+02 

1067+01 

0205+02 

2002+01 

0604+01 

0702 


Meaning  or  Use 


Index  of  inner  DO  loop  to  account  for  all  fine 
groups  in  each  broad  group 

Index  of  the  inner  DO  loop  used  to  load  the 
cross  sections  into  the  scattering  matrices. 

It  is  the  k  index  of 

Index  of  the  fourth  of  four  nested  DO  loops. 

For  scatter  from  broad  group  n  to  broad  group 
m,  it  accounts  for  all  fine  groups  in  m 

Index  of  the  inner  DO  loop  for  broad  group  out¬ 
put.  It  indexes  through  the  fine  groups  within 
any  given  broad  group 

Both  a  gate  and  an  index  to  identify  the  first 
fine  group  with  a  flux  less  than  10~12  (if  any) 

The  group  number  of  the  lowest  fine  group,  to 
which  scatter  occurs,  from  a  given  broad  group 

An  index  used  with  the  variable  k  in  SS  1122. 

It  assists  in  proper  scanning  of  the  scattering 
matrices  (Example:  prevents  scanning  on  or 
below  the  diagonal  for  (n,n*)  or  (n,2n)  cross 
sections ) 

A  gate  to  cause  KZ  to  increment  properly  after 
reading  in  the  first  elastic  (P0)  and  (PI) 
values  into  the  scattering  matrices:  AKZ  =  2 
for  elastic  scatter,  AKZ  »  1  for  all  others 

A  counter  to  sort  the  cross  sections  into  the 
proper  matrices 

Index  of  the  second  of  four  nested  DO  loops. 

For  scatter  from  group  n  to  group  m,  it  is  the 
value  of  m 

Index  used  to  convert  the  k  index  of  o •  ^  for 
elastic  (P0)  and  (PI)  cross  sections.  3 These 
cross  sections  are  read  into  the  lower  half  of 
SSN  and  S2N 

A  subscripted  variable  for  the  numbers  of  the 
fine  groups  that  form  the  boundaries  of  the 
broad  groups 

A  code  number:  indicates  type  of  inelastic 
scatter.  1  =  (n,n*),  2  =  (n,2n),  3  =  both, 

4  =  neither 


130 


MMP 


Name  Location 

LS  0702 

2002 

M  1112-02 


MOD 


0701 


MXL 


0604+04 


0004+01 


NA 


0702 


NBG 


0604+01 


NRNUC  0004+01 


P 

PS] 


0701 

0702 


0054 


RES 


0704 


Meaning  or  Use 

The  total  number  of  entries  in  the  nuclide 
cross  section  library  starting  with  the  second 
library  card  (does  not  include  epithermal  \>Of 
of  fission  nuclides) 

Index  used  to  convert  the  j  index  of  o j for 
elastic  (P0)  and  (PI)  cross  sections.  These 
cress  sections  are  read  into  the  lower  half  of 
the  SSN  and  S2N  matrices 

An  address  variable  to  cause  the  program  to 
branch  to  the  proper  statement  number  after 
clearing  the  matrices .  It  is  a  statement  num¬ 
ber  that  can  be  varied.  Used  with  GO  TO  M 
statements 

Is  nuclide  moderator,  fuel,  or  other?  Not  used 
in  chain  2  however  it  must  be  read  in  so  that 
correct  number  is  read  into  DENS 

The  number  of  broad  groups  requested  plus  one. 
This  allows  for  output  of  thermal  group  values 

Index  of  a  DO  loop  (SS  0004+01  to  SS  0055)  used 
to  read  through  all  nuclides  of  the  mixture 

A  code  number:  Is  epithermal  absorption  appre¬ 
ciable:  0  =  NO,  1  =  YES 

The  total  number  of  broad  groups  for  chain  2 
output 

The  total  number  of  nuclides  in  the  mixture. 
COMMON  from  chain  1 

Nuclide  density  in  atoms  per  barn -cm 

The  average  logarithmic  energy  decrement  £ .  Not 
used  in  chain  2 

A  summing  variable  for  the  group  to  group  elas¬ 
tic  (PI)  flux  weighted  cross  sections.  See  use 
of  A  in  SS  0054-01 

A  code  indicating  whether  or  not  a  resonance 
calculation  was  performed  in  chain  1.  If  yes 
(RES  =1)  the  resonance  absorption  will  be  in¬ 
cluded  in  calculating  the  absorption  cross 
sections  of  the  resonance  nuclide.  COMMON  from 
chain  1 
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f  -  *1  <■"  n.r 


Name 

RIEFF 

SA 

SAM 


SAMM 


SIGA 

SIGS 

SIGTR 

SOURCE 

SR 

SSN 

STR 


Location 

0700 

1005 

0003-01 

0701 

1122+01 

0315 

0701 

1122+02 

1109+01 

0702 

0702 

0314+02 


Meaning  or  Use 

The  effective  resonance  integral  of  the  reso¬ 
nance  nuclide(s).  COMMON  from  chain  l 

Microscopic  epithermal  absorption  cross  section 

The  sum  of  all  fine  group  fluxes  in  a  given 
broad  group 

The  first  five  alphameric  characters  of  the 
nuclide  name 

A  summing  variable  for  group  to  group  (n,n*) 
flux  weighted  cross  sections.  See  the  use  of 
A  in  SS  0054-01 

Dummy  storage  for  th«  f,.*id  background  cards 

The  second  five  alphameric  characters  of  the 
nuclide  name 

A  summing  variable  for  group  to  group  (n,2n) 
flux  weighted  cross  sections.  See  the  use  of 
A  in  SS  0054-01 

A  summing  variable  for  the  broad  group  flux 
weighted  vof 

Microscopic  thermal  absorption  cross  section 

Microscopic  thermal  scattering  cross  section 

Thermal  transport  cross  section:  C(l-jj)oa+oa]* 
DENS 


0004  The  fraction  of  the  total  source  in  each  broad 

group 

0901  The  broad  group  flux  weighted  removal  cross 

section .  The  absorption  contribution  is  added 
in  SS  0045-01 

1009+01  Scattering  matrix  (12  x  13)  for  (n,n*)  and 

elastic  (P0)  cross  sections.  Upper  triangle 
(J,  k+1)  is  for  (n,n*).  Diagonal  and  lower 
triangle  (LS,  LE)  is  for  elastic  (P0) 

0054+02  The  broad  group  flux  weighted  common  transport 
cross  section.  The  absorption  contribution  is 
added  in  SS  0045 


Name 

Location 

Meaning  or  Use 

S2H 

1010+02 

Scattering  matrix  (12  x  13)  for  (n,2n)  and 
elastic  (PI)  cross  sections.  Upper  triangle 
(J,  k+1)  is  for  (n,2n).  Diagonal  and  lower 
triangle  is  for  elastic  (PI) 

TERMU 

0702 

— ■  th 

Thermal  (1-u)  of  the  n  nuclide 

TT 

0054+01 

A  summing  variable  for  the  flux  weighted  group 
to  group  total  transfer  cross  sections.  See 
the  use  of  A  in  SS  0054-01 

VF 

0702 

Thermal  vo^  of  the  n**1  nuclide 
f 

XBAR 

1130 

The  average  normalized  thermal  neutron  velocity 
x.  COMMON  from  chain  1 

XSAM 

0049 

Temporary  storage  for  the  entries  of  the 
nuclide  cross  section  library 

Appendix  D 


Chain  Two  Source  Deck 


C  THIS  IS  PART  2  OF  OLD  BARNYARD.  CROSS  SECTION  COLLAPSING 

t  J08  OLD  BARNVARD  CHAIN  2 

$  OR  Ci  I N  10970 

C  NOTICE  COMMON  STATEMENT  IS  SAME  AS  IN  CHAIN  1 

COMMON  LGP( 13) .KE.FLXf 11 ) »Ut 11 ) .SOURCE! 11 ) , DENfcJR ♦NRNUC , TEMP TR , XBAR 
1 , PS  I GM. R I EFF  ,RES 

DIMENSION  S A  C 1 1 »  *  STR  <  11) .SR  ill  )  »XSAM(200) .SSN( 12 . 1  3 ) , S2N (  1  2 . 1 3  )  ♦ 
1VSFI 1 1 ) 

TvPt  604 

604  FORMAT ( / 1 6HBFG I N  COLLAPSING) 

C  READING  IN  NUMBER  OF  BROAD  GROUPS  AND  THEIR  BOUNDARY  FINE  GROUPS 

READ  107C.N8G. (LGP! J) . J=2,NBG+1 ) 

1070  FORMAT (12N) 

C  CHECK  FOR  2ERO  FLUX  AND  REv'ISE  THE  GROUP  STRUCTURE  IF  NECESSARY 

K=NBG 
MXL=NBG+1 
LGPINBG+? ) =12 
KK  =  0 

DO  305  J=1 *LGP ( NSG+ 1 ) 

IF(KK)3G1.301.305 

301  IFCFLX; J1-1.0E-12 ) 303.305, 30S 
303  KK=J 

DO  305  1=2. KK 

I F ( J- ( LGP ( 11+11)310.310.305 
310  NBG= I -1 
3C7  1=12 

305  CONTINUE 

IF(KK) 308.303.309 
309  TYPE  110C.KK.N8G 

1100  FORMAT ( /37UFLUX  APPROXIMATELY  ZERO  IN  FINE  GROUP » I  3 » 2X , 2 7HT0  PREVE 
1NT  DIVISION  RY  ZERO. /39HWHEN  FLUX  WEIGHTING  CROSS  SECTIONS  ONLY. I 
23.2X.28HFPI THERMAL  BROAD  GROUPS  WILl,/,71HBE  output,  the  FOLLOWING 
3FINE  GROUP  BOUNDARIES  WILL  BE  USED  ACCORDING  TO, / . 14HSENSE  SWITCH 
41) 

TYPE  1 j  0 1 . I LGP  t J ) t  J=2  »K  +  1 ) .12 
TYPE  1102, ( LGP  t  J ) , J=2,N6G+1 ) ,  12 

1101  FORMAT (2HON.3X, 121 3) 

1102  FORMAT ( 3 HOFF . 2X , 12  I  3 ) 

TYPE  1103 

1103  FORMAT! 30HWHERE  12  IS  THE  THERMAL  GPOU°. ♦ / » 3 1HSE T  SENSE  SWITCH  l.P 
1PESS  START.) 

PAUSE 

IF (SENSE  SWITCH  1)1105,1106 
1106  LGP ( NEG+2 ) = 1 2 
K  =  N3G 
PSIGM^O. 
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r- 


M  X  L  =  H  v  G  -r  I 
GO  TO  308 
1105  LGP ( K  +  2  ) =12 
MXL=K+1 
308  LGP ( 1 ) =0 

IF(PSIGM) 1107,1108, 1107 
1108  1=0 

TYPE  1909 

1909  FORMAT ( 44HTHFRF  WILL  RF  NO  THERMAL  BROAD  GROUP  OUTPUT.) 

GO  TO  1110 
110"  1=1 

1110  PUNCH  1111 ,MXL-1 +1 *NRG+I , { LGP ( J ) »J=2»K+2) 

TYPE  111  1 ,MXL-1 *1 *NBG-i . (LGP!  J)  *J=2.K+2) 

1111  FORMAT ( 33HT HE  TOTAL  NUMBER  OF  BROAD  GROUPS  IS* 1 3»/ ♦ 1 SHOUT PUT  WILL 
1BE  FOR*  I3*2X»1 2HBR0AD  GROUPS* /.  28HTHE  BOUNDARY  FiNF.  GROUPS  ARE.12 
2  I  3  » / ♦ 30HWHESE  12  IS  THE  THERMAL  GROUP.) 

NORMALIZING  THE  FINE  GROUP  FLUXES  IN  EACH  BROAD  GROUP  AND 
ADDING  UP  THE  SOURCE  IN  EACH  BROAD  GROUP 
DO  2  J- 1 ,NBG 
SA ( J ) =0. 

DO  2  K  =  LGP ( J )  +  1  *  LGP ( J+l ) 

2  SA! J)=SA( J)+FLX(K) 

DO  A  J= 1 *NBG 
SAM=0. 

DO  3  <=LGP(J)+1,LGP( J+l) 

SAM=SAM+SOURCF(K) 

3  FLX(K)=FLX(K)/SA( J) 

4  SOURCE! J)=SAM 

C  ONE  TRIP  THROUGH  THE  DO  LOOP  FOR  EACH  NUCLIDE 

DO  55  N=1 ,NRNUC 

C  CLEARING  THE  MATRICES  AND  SETTING  SCWE  TERMS  TO  ZERO  FOR  LATER  USE 

ASSIGN  704  TO  M 
I RES=0 

1112  DO  1002  J=  1  ,  1  2 
DO  1002  K= 1 » 1 3 
SSN ( J,K ) =0 • 

1002  S2N ! J. K ) =0 • 

1113  DO  1000  J=l»ll 
STR ( J ) =0. 

SR! J)=0. 

SA ; J ) =0. 

1000  VSF ( J) =0. 

1114  A=0. 

SAMM=0. 

SAM=0* 

R  =  0. 

TT  =  0  . 

GO  TO  M 

C  READING  IN  THE  NUCLIDE  DATA  CARDS 

704  IF(RES)1115, 701*1115 

701  READ  1061  *SAM,SAMM*  ATOMS  *P»M.OD*DENS 
GO  TO  702 

1115  READ  1061  * S AM, SAMM, ATOMS  »  P , MOD*  DENS  *  I  RES 
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K6i  FOR*  .*T(5X,3A5,lON> 

70?  R  E AD  1 0  6 i »  S  A  m , S  A  MM , A  T  0M  S  »  L  S  »  N  A  »  L  Q  *  P  S  I » S I GA , DEL T A  ,5 I GS * TERMU ♦ VF * F I S 
punch  2CC » s am , samm ,  atoms  * p 

200  format ( //l 1HNUCI IDE  IS  ,3A‘>,//25h  ITS  NUMBER  DENSITY  IS.1PE11.4 
1 ,  12H  PER  6ARN-CM,//  ) 

C  CHECK.  I NG  FOR  MICRO  OR  MACROSCOPIC  CROSS  SECTIONS 

I  F { DENS-1 • 161.62 ; 6  I 

61  DENR=o 
PUNCH  2^2 

202  FORMAT ( 46HMACR0SC0P I C  CROSS  _-CTIONS.  IN  PER  CM,  FOLLOW.) 

GO  TO  A 9 

62  PUNCH  201 

201  F0RMAT(45HMICR0SC0PIC  CROSS  SECTIONS.  IN  BARNS.  FOLLOW.) 

C  READING  ’N  THE  WHOLE  CROSS  SECTION  LIBRARY  FuR  THE  NUCLIDE 

40  READ  "1067, IXSAM(j).J-l.LS) 

1067  FORMAT (6N) 

KZ  =  0 
KKZ  =  0 


C  CHECKING  FOR  ABSORPTION 

IF(NA) 1006. 1006, 100A 
1004  DC  1005  J*l.ll 
KZ  =  KZ-*-l 

)oos  SAt J)=SA( J)+XSAM(KZ )*CLX( J)*DENS 
I  F ( IRES! 700,1006,700 

700  SA(9)=SA(9)+( ( R I EFF /LOGF ( 61 .442 / 1 . 1 254 ) ) * . 1*FLX ( 9 } *DENS ) 

.  SAt  10)=SA(  1 G  )  -*-  (  (RIEFF  /  LOGF  (61.A42/1.1254  )  )  *.  9*FLX  (  10  )  *DENS  ) 

1006  KZ=KZ+1 
IN=XSAMIKZ) 

KZ=KZ+1 

IIN=XSAV(KZ) 

DO  1017  J= 1 » IN 
KK=J+I IN 

IF (KK-12) 1008, 1008, 1007 

1007  KK= 1 2 

1008  DO  1017  K=J.KK 

30  TO  ( 1009,1010,1009,1011  )  ,LO 

1009  KZ=KZ+1 


r  INFLASTIC  CROSS  SECT  I  ON .GROUP  TO  GROUP 

SSN ( J »K  +  1 ) =S$N ( J  » K+ 1 ) +XSAM ( KZ ) 

GO  TO  1017 

1010  KZ=KZ+1 

C  N-2N  CROSS  SECTIONS,  GROUP  TO  GROUP 

S2N(  J,K+1  )  =S2N  (  J  ,  <-:■  1  )  +XSAM  {  KZ  ) 

GO  TO  1017 

1011  IF(KKZ)2001, 2000, 2001 

2000  KZ=KZ+1 
KKZ=  1 

GO  TO  2002 

2001  KZ  =  KZ  +2 


C  PUTTING  ELASTIC  (P 0)  AND  ELASTIC  (PI)  IN  THE  BOTTOM  OF  SSN  AND  S2N 

2002  LS=13-J 
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LE= 13-K 

SSN!LS*LE) =SSN!LS.L£)+XSAM<KZ ) 

~  52N ( LS»LE  )"=S2N  ( LS~»LE  )+XSAK"(  KZ+1 )  /  3* 

1017  CONTINUE 

GO  TO  ( 1018.1018.1019.1020) »LQ 

1018  LQ=4 

GO  TO  1006 

1019  LQ=2 

,  GO  TO  1006 

1020  PUNCH  205 

205  FORMAT! //7HSCATTER . 34X .  7HELAST I C* 6X . 7HELAST I C . 7X . 5HTOTAL/8HFROM  T 
10*6X*9HINELASTIC*7X*  4HN- 2N *  9X  *4H ( PO ) * 9X  *4H ( PI ) ,6X » 8HTRANSFER/ ) 

C _ OUTPUTT  I NG  _GROUP_TO  GROUP  VALUES _ _ 

DO  5 3  1  =  1  *  NBG  '  '''  . .  '  '  ~ 

DO  53  L= I *MXL 
ASSIGN  1120  TO  M 
GO  TO  1114 

1120  DO  54  J=LGP{ I ) +1  *  LGP ! 1+1) 

KKaLGP ( L ) +1  __  _  _  _ 

. IF(KK-J)Tl21. 1122.1122 

1121  KK= J 

1122  DO  54  K=KK »LGP ( L+l ) 

SAM=SAM+SSN ( J»K+1 )*FLX! J }*DENS 
SAMM3SAMM+S2N ( J»K+1 ) *FLX ( J ) *DENS 

_ LS- 13-J _ 

LE*=  13-K 

A*=A+SSN(LS,LE)*FLX(  J)*DENS 
54  R=R+S2N<LS.LE)*FLX( J)*DENS 
TT-TT+SAM+2«*SAMM+A 
STR ( I ) =STR ( I J+SAM+SAMM+A-R 

.  _ IF(I-L  >900*901  ♦900  _ _  _ 

901  SR < I) *SR ( I ) -SAMM 
GO  TO  53 

900  SRI  I )=SR( I J+SAM+SAMM+A 
53  PUNCH  103.1 *L»SAM,SAMM*A*R*TT 

C  THERMAL  GROUP  TO  GROUP  VALUES _ . _ 

SIGS*SIGS*DENS 
IF! PS  I GM >56*903*56 

56  PUNCH  103»MXL.MXL*0.*0.»SIGS» I1.-TERMU>*SIGS.SI6S 
103  FORMAT! 13. 15. 5F13-5> 

C  CALCULATING  AND  CJTPUTT1NG  BROAD  GROUP  VALUES 
903  PUNCH  220 

220  FORMAT ( //5HGROUP .5X.5HSIGTR.9X.4HSI GR»  9X .4HSIGA  »  5X . 6HNUSIGF  »6X. 
16HSOURCE ) 

421  IF(FIS)4040. 4040*4041 

4041  READ  1067»!VSF!I), 1=1*11) 

4Q4Q  DO  46  J= 1 «NBG _ 

ASSIGN  1109  TO  M 
*  GO  TO  1114 

1109  DO  45  K=LGP( Jl+l.LGP! J+l ) 

SAMM=SAMM+VSF ( < ) *FLX ( K ) 

A=A+SA!K) 

SR ( J  >  =SR ( J ) +SA ( K ) 

45  STR! J)=STR( JJ+SAIKT 


T 


46  PUNCH  102.J.STRI  J)  ,SR<  J)  ,  A  .DENS*  SA.-1M  ,  SGUKCl  <  o  ) 

C  OUTPUTTING  THERMAL  BROAD  GROUP  VALUES 

I F ( PS  I GM  )  3 1 4  » 3 1 5  » 3 1 4 

314  VF= (VF/XRAR ) *DENS 

S I GA= ( S I GA /XPAR ! *DENS 
SIGTR=(SI6S*TF.RMU+SIGA) 

PUNCH  102.MXL .SIGTR.SIGA.SIGA.VF 
102  FORMAT (  I4.5F12.5) 

C  READING  AND  IGNORING  THE  GRID  BACKGROUND  CARD  IN  EACH  LIBRARY  DECK 

315  READ  1061,  SAM 
55  CONTINUE 

IF(PSIGM)  1130. 1093,1130 
1130  PUNCH  lll.Xi.AR 

111  FORMA  T  ( //48HMAXV/ELL-BOLTZMAN  FACTOR  =  1.126*  AVERAGE  X  =  »N) 
1093  TYPE  1063 

1063  FORMAT  (/ 15HE.ND  OF  PROGRAM.  /  ) 

PUNCH  1065 
1065  FORMAT (//IX) 

STOP 

END 
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List  of  Neutron  Source  Deck 


.5  1.  1.5  ?.  2.5  3.  5.  S.  12.  16.  17. 

THE  SOURCE  USED  IS  'J235  FISSION  (CRAN6ERG  SPECTRUM). 

2  .023023  .10824  .21044  .23139  .18048  .11483  .12439  .0072094  .0  .0  .0 
THE  SOURCE  USED  IS  U233  FISSION. 

3  .021495  .10472  .20811  .23202  .18252  .11674  .1267  .0073845  .0  .0  .0 

THE  SOURCE  USED  IS  PU239  FISSION. 

4  .024567  .11047  .21062  .22996  .17917  .11412  .1239  .0071986  .0  .0  .0 
THE  SOURCE  USED  IS  PU241  FISSION. 

5  .026852  .11615  .215  .22949  .11077  .1190  .0068574  .0  .0  .0 
THE  SOURCE  USED  IS  CF252  FISSION. 

6  .046142  .14786  .2292  .21936  .15472  .09531  .099137  .0055827  .0  .0  .0 
THE  SOURCE  USED  IS  PO-RE  ( WHI TMORE-BAKER > . 

7  .2342  .357  .261  .1025  .036  .0085  .0  .0  .0  .0  .0 
THE  SOURCE  USED  IS  PO-BE  ( COCHRANE-HENR Y > . 

8  .257  .369  .251*4  .079  .028  .010  .0  .0  .0  .0  .0 
THE  SOURCE  USED  IS  RA-BE  (HILL). 

9  .239  .365  .189  .109  .065  .033  .0  .0  .0  .0  .0 
THE  SOURCE  UCED  IS  YOUR  OWN  AS  SHOWN  BELOW. 

10.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 

THE  SOURCE  USED  IS  A  UNIT  SOURCE  IN  GROUP  NO. 

10  .0  .0  .0  ,0  .0  .0  .0  .0  .0  .0  .0 


l 
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where  m-n  is  the 
scattering  cross 
in  the  matrix* 


location  of  the  inelastic  (Po) 
section  from  group  m  to  group  n 


Figure  19»  Inelastic  and  silastic  (Po)  Scattering  Matrix 
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